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PREFACE 


Active  control  technology  has  steadily  emerged  and  offers  a means  for  improving  perfor- 
mance and  operational  flexibility  and  of  changing  air  vehicle  design  concepts,  procedures  and 
methods.  The  use  of  active  control  devices  on  aircraft  appears  beneficial  in  at  least  six  different 
areas:  ( 1 ) quasi-steady  load  reductions,  (2)  flutter  suppression,  (3)  reduced  fatigue  loading, 

(4)  improved  stability  and  control,  (5)  improved  ride  qualities,  and  (6)  reduced  margins.  The 
Flight  Dynamics  Laboratory  of  Wright  Patterson  Air  Force  Base  played  a leading  role  in 
initiating  the  study  effort  on  active  controls,  and  since  then  groups  in  Europe  as  well  as  the 
United  States  have  made  marked  advances  in  developing  the  concepts,  in  building  hardware, 
and  in  testing,  both  in  windtunnels  and  in  full-scale  flight. 

The  problem  area  is  appropriate  for  joint  consideration  by  the  NATO  nations,  and  work 
appears  to  follow  logically  under  the  cognizance  of  the  SMP,  FMP  and  GCP.  In  the  October 
1 973  meeting  (the  Hague),  the  SMP  held  a plenary  session  on  Active  Control  Devices,  in  which 
five  papers  were  presented.  At  their  April  1975  meeting  in  Brussels;  the  SMP  held  a specialist 
meeting  on  Flutter  Suppression  and  Structural  Load  Alleviation,  at  which  eight  papers  were 
presented.  In  October  1974,  the  FMP  and  GCP  held  a four-day  symposium  on  the  subject 
entitled  “Impact  of  Active  Control  Technology  on  Airplane  Design”. 

At  the  April  1977  meeting  (Lisbon)  the  SMP  held  a specialist  meeting  on  Structural 
Aspects  of  Active  Controls,  at  which  seven  papers  were  presented.  These  Proceedings  are  a 
compilation  of  the  papers  presented.  The  specialist  meeting  covered  here  was  organized  and 
conducted  by  the  Ad  hoc  Group  on  Structural  Aspects  of  Active  Controls. 

It  is  hoped  that  effort  within  the  SMP  has  helped  in  the  dissemination  of  knowledge  on 
the  subject  and  that  these  Proceedings  will  help  promote  useful  thoughts  and  contribute  to  the 
development  of  means  for  increasing  the  safety,  performance,  and  utility  of  aircraft. 


John  C.HOUBOLT 
Chairman,  Ad  hoc  Group 
on  Structural  Aspects  of 
Active  Controls 


ili 


CONTENTS 


Page 

PREFACE  Wi 

Reference 

A PRACTICAL  OPTIMUM  SELECTION  PROCEDURE  FOR  A MOTIVATOR  IN  ACTIVE 
FLUTTER  SUPPRESSION  SYTEM  DESIGN  ON  AN  AIRCRAFT  WITH  UNDERWING  STORES 

by  M.R.Turner  and  C.G.  Lodge  1 

IMPACT  OF  A COMMAND  AND  STABILITY  AUGMENTATION  SYSTEM  ON  GUST 
RESPONSE  OF  A COMBAT  AIRCRAFT 

by  K.D.Collmann  and  O.Sensburg  2 

ACTIVE  FLUTTER  SUPPRESSION  ON  AN  AIRPLANE  WITH  WING  MOUNTED  EXTERNAL 
STORES 

by  H.Honlinger  3 

AIRPLANE  MATH  MODELING  METHODS  FOR  ACTIVE  CONTROL  DESIGN 

by  K.L.Roger  4 

CONSISTENCY  IN  AIRCRAFT  STRUCTURAL  AND  FLIGHT  CONTROL  ANALYSIS 

by  R.W.Schwanz  5 

YC-I4  CONTROL  SYSTEM  REDUNDANCY 

by  W.T.Hamilton  . 6 

ETUDE  EN  SOUFFLERIE  D’UN  SUPPRESSEUR  DE  FLOTTEMENT 

par  R.Destuynder  7 


t PBAfyrTRAT.  OfTimm  smanrioH  fBOCHlDBB  TOE  A 


MOTIYATOE  IH  ACTIVE 


yUTHTHl  SDPFiaSSIOH  SYSTEM  IBSIC8 


OH  AM  AIBCTArr  WITH  nwiWBWTIT.  STOBBR 


N.  R.  Tuzner 

British  Alxorsft  Corporation, 
CooDsrclal  Alroraft  Siylslon, 
Hilton,  Bristol,  Boglsnd. 


C.  G.  hodgs 

British  Alroraft  Corporation, 
Militaxr  Aircraft  Siylsion, 
Warton,  Lanos,  ftigland 


A^stQot 

Tbsorotloal  aotlTo  flutter  control  of  a Tsrlahls  swoop  wing  with  oxtomal  atoros  with 
four  ooahlnationa  of  storo  oonfiguxatlon/wing  swoop/Haoh  Aiahor  was  studied.  Elootrloally 
aodifiod  outputs  of  a struotnro  aountod  txansduoor  wore  usod  to  driwo  an  auxlliaxjr  oontrol 
surfaoo  on  tho  wii«  or  storo.  The  host  tranaduoor/foroo  positions  on  tho  wing  sad  storos 
wore  found  using  Hjnpilst  Plots  and  roproaontlng  ths  oontrol  surfaoo  loads  hp  point  forces. 

Tho  object  was  to  see  if  a oosMon  aetlTO  flutter  oontrol  systea  using  a control 
surface  on  ths  wing  oould  be  found  for  a range  of  stores,  Maoh  Buabers  and  wing  sweep  angles. 

Blffloultles  were  due  to  two  instabilities  with  close  Ar^uencles  in  two  of  tbs  configurations 
and  Tory  low  darings  in  sons  of  ths  stable  nodes. 

1.  IHTROllOCTIOH  J 

The  use  of  aotire  flutter  oontrol  could  hare  algnlfloant  adrantages  oyer  ths  astabllahad  I 

asthods  of  structural  nodlfioatlon  and  adding  concentrated  naaees.  Boweyer,  if  no  suitable  control  j 

surface  is  already  ayallable  to  act  as  a notlyator  for  an  aotlye  flutter  oontrol  systea,  an  anzlliaxy  1 

oontrol  surface  must  be  introduced.  j 

This  paper  describes  a practical  optlnnn  selection  procedure  for  the  position  of  ths  auxiliary 
oontrol  surfaoe  and  ths  feedback  transducer.  The  procedure  is  Illustrated  by  applloatlon  to  the 
oontrol  of  flutter  of  a yarlable  sweep  aircraft  with  underwlng  stores. 

2.  OBJECTIVE  I 

A theoretical  study  was  aade  of  the  use  of  aotlye  controls  to  suppress  ths  flutter  of  a | 

yarlable  sweep  wing  (Figure  1)  with  the  four  oases  of  store  configuration,  Mach  Nuabsr  and  wing  sweep  j 

angle  shown  in  Figure  2.  Ths  wing  had  no  suitable  oontrol  surfaoe  already  ayallable  which  oould  j 

be  used  for  actlye  flutter  oontrol.  Therefore  the  addition  of  an  auxiliary  control  surface  on  the  j 

wing  or  store  was  necessary.  j 

The  object  was  to  find  the  aost  effeotlye  control  surfaoo  position,  transducer  position  and  j 

feedback  law  to  suppress  flutter  for  each  of  the  foiar  oases.  It  was  also  intended  to  find  the  best  j 

w,,,  system  using  a control  surface  on  the  wing.  This  la  logistloally  more  appealing  than  to  use 
oontrol  surfaces  on  each  store.  j 

3.  mathematicat.  wnmn.  j 

3.1  Basic  Wing-Plus-Store  Characteristics 

Sixteen  branch  nodes  for  each  wing  and  storo  combination  were  used.  These  were  zioiDSllsed  and 
the  first  el^t  wore  used  for  the  aotlye  flutter  oontrol  study.  The  fuselage  was  aseumsd  fixed 
Inboard  of  tbs  wing  plyot.  Two  dlmsnslonal  aerodynamic  dariyatlyes  wore  used  in  a strip  theory  sense 
for  ths  wing  flutter  aerodynasdos.  These  aerodynamlo  derlyatlyes  were  ohosen  to  giye  the  same 
flutter  speeds  and  similar  suborltloal  oharaoterlstlos  to  those  giyen  by  three  dimensional  lifting 
surface  theories. 

Because  (l)  sons  stable  roots  had  axtremsly  low  aerodynamio  damping 

(li)  the  introduction  of  1%  structural  damping  into  the  first  three  modes  excluded  most 
of  ths  basic  store  flutters 

ths  study  was  dons  with  no  structural  damping  in  the  first  three  normal  modes  and  1%  structural 
dam>ing  in  modes  U to  8. 

Figure  3 shows  the  basic  flutter  oharaoterlstlos  for  the  four  cases.  Beoause  the  instabilities 
occur  in  the  first  three  roots,  which  are  close  in  frequency,  only  these  roots  are  shown  in  tho 
flutter  plots.  Hotios  that  for  oases  1 and  3 there  are  two  unstable  roots  close  in  frequency. 
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3.2 


togtili«£j_Con^2l_^22lS£ 


To  otaooM  the  heet  position  for  the  auxiliary  control  aurfaoe,  the  control  aurfaoe  foroee 
were  represented  by  a point  force  P 


where 


P = ipV^KofWW 

's^V  ° etreaa  dynamic  preeaure 

K,  = overall  feedback  gain 
= feedback  filter 
W s transducer  output 


This  assumes: 

(l)  tbs  actuator  transfer  function  is  unity 
(ll)  the  control  surface  rotational  stiffness  is  infinite 
(ill)  the  frequency  dependant  part  of  the  control  surface  aerodynamics  can  be  Ignored 


and 

(iv)  the  inertia  foroee  due  to  control  surface  acceleration  are  negligible. 

Therefore,  the  generalised  excitation  forces  due  to  the  control  surface  force  acting 
on  the  modal  system  can  be  sliqply  represented  by: 

lA^P 

Where  is  the  vector  of  modal  deflections  at  the  point  force  position. 

The  Justification  for  using  a sli^le  point  force  to  represent  the  aerodynemic  loading  due  to 
control  surface  rotation  is  as  follows.  We  consider  strip  theory  representation  of  the  control 
surface  aerodynaialcs  to  bo  adequate  for  preliminary  design  calculations  on  a high  aspect  ratio 
wing.  If  strip  theory  was  used,  the  aerodynamic  stiffness  forces  would  be  represented  by  a 
point  force  on  eaoh  chordwlse  strip.  Therefore,  if  the  strips  were  the  width  of  the  control 
surface  we  would  have  the  single  point  force  representation  used  in  this  study. 


Aesuaqjtions  (iii)  and  (iv)  are  Justified  in  the  interest  of  providing  an  economic  preliminary 
design  procedure,  while  assumption  (i)  can  be  cosq>ensated  for  by  a simple  filter  if  nooeeeaiy. 

Generally,  the  error  duo  to  the  point  force  aasuiiq>tlon  will  only  be  significant  at  relatively 
high  frequencies  where  it  is  difficult  to  accurately  model  the  aeroelastlc  data  and  actuator  anyway. 

Obviously  the  most  favourable  positions  given  by  the  point  force  study  should  be  checked  with 
a full  representation  of  the  control  surface  forces.  The  actuator  transfer  function  and  impedance 
will  emerge  as  requirements  from  the  design  study. 


1,.  yUTTER  COHTROL  SYSTEM  BESIGH 


U.1  Design  Reoulrements 

(1)  A design  speed  is  chosen  below  whloh  the  aircraft  must  iu>t  be  unstable.  This  is  typically 
above  the  maxtmin  fll^t  speed  of  the  aircraft. 

(ii)  Stability  nsist  be  maintained  in  a 2$  ft. /sec.  tuned  (1 -cosine)  gust  at  the  design  speed. 

(iii)  Overall  aircraft  stability  and  control  oheraoterietios  must  not  be  affected  by  the  flutter 
control  system. 

U.2  Design 

To  achieve  these  design  requirements  with  a sliq>le  flutter  control  system  with  maxisium 
reliability,  low  sensitivity  and  Tnlniimm  weight  penalty,  the  following  are  aimed  for: 

(l)  Only  one  transducer  and  one  control  surface  will  be  used. 

(il)  Adequate  closed  loop  gain  and  phase  stability  margins  are  to  be  available  at  all  frequencies 
and  speeds.  This  is  to  cover  manufacturing  and  wear  induced  tolerances  in  the  feedback 
hardware  - particularly  in  the  actuator.  It  may  also  be  iwcessary  to  cover  variations  in 
the  basic  aircraft  transfer  function  either  from  the  predicted  value  or  from  variations  in 
flight  conditions  such  as  fuel  distribution  or  altitude. 

In  this  study  the  requirements  will  bet6db(i.e.  «2  and  +2)  and  160°  because  these 
are  the  margins  in  comsion  use  in  control  systems  eiigineering. 
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(ill)  Only  flrat  ordar  flltara  ar«  to  be  ua«d  In  ttaa  fMdbnok  law  to  obtain  tho  naoaaaaxy  ptaaaa  j 

ohangaa  and  attanuationa.  Notoh  flltars  will  not  ba  oaad  to  raaowa  tha  raaponaa  of  an  ! 

undaalrabla  atxuotuxal  noda  aa  thla  oxaataa  a aanaltlwlty  pxoblan  dua  to  altbar  Inaoourata  ] 

pzadlotlon  of  tha  nodal  fraquanoy  or  ita  variation  with  fU^t  oondltlona  or  foal  atata.  1 

(iv)  To  avoid  aanaitivlty  to  unoartalntlaa  In  tha  control  aurfaoa  aarodynaaloa  tba  actuator  ’ 

tranafar  function,  raaponaa  of  tha  hi|^r  fraquanoy  a»daa  in  tha  faadback  aignal  will  ba  j 

filtarad  out.  j 

(v)  No  faadbaok  gain  and  phaaa  aohaduling  with  apaad  will  ba  oonaldarad.  | 

(vl)  A oomion  faadbaok  ayatam  whloh  can  ba  uaed  for  all  four  caaaa  will  bo  aaarohad  for.  ] 

U>3  Bealan  Prooadura  | 

For  aaoh  oaaoi  | 

(l)  Araduoa  If/qulat  Plota  at  tho  daalgn  apaad  for  a rango  of  ooabinatlona  of  point  load  hA 
tranaduoer  poaltlona. 

(ll)  Chooaa,for  each  Nyqulat  Plot  from  (i),  a faadbaok  gain  and  a filter  whloh  naka  tho  unatabla 
modaa  atabla,  and  tha  atabla  modaa  ramain  atabla,  with  tha  boat  poaalbla  gain  and  phaaa 
margina. 

(ill)  Choose  tho  bast  force/ transducer  position  combinations  talcing  into  account  the  Sealgn  Alms.  j 

(iv)  For  each  of  these  chosen  oomblnatlons  produce  I^quist  Plots  for  speeds  below  tha  design  | 

spaed  to  sea  If  a simple  constant  control  law  would  satisfy  tha  Design  Alms  at  all  speeds.  \ 

(v)  Find  the  marlmum  force  required  (from  which  the  slaa  and  rotation  of  the  control  surface,  j 

and  tha  actuator  oharacterlstlos  such  as  rata  limit,  are  obtainable)  by  calculating  tha 
response  of  the  controlled  structure  to  a tuned  25  ft./seo.  (l-ooslna)  guat  at  tha  design  speed. 

U.U  Point  Force  and  Transducer  Positions  Inveatlaatad  I 

The  point  force  positions  considered  on  tha  wing  (Figure  U)  ware  five  spanwlse  stations  at 
5$%  chord  (representing  forces  due  to  2$?^  chord  trailing  edge  control  surfaces)  together  with  twanty 
chordwlse  positions  near  the  wing  tip.  The  wing  transducer  positions  (Figure  U)  were  fifteen 
vertically  and  five  meaaurlng  streamwise  Inoldanoa. 

Tha  force  and  tranaduoer  positions  on  tha  stores  (Figure  5)  ware  vertical  and  lateral  at 
each  end  of  the  store  with  two  rotation  transducers  msasuring  store  pitch  and  yaw.  System  designs 
using  store  excitation  have  been  produced  purely  for  ooiq>arleon  purposes,  since  It  is  well  Imown 
that  design  of  store  excitation  systems  la  relatively  straightforward. 

From  operational  considerations  It  la  best  to  have  both  the  control  force  and  the  transducer 
on  the  wing.  Next  best  Is  the  force  on  tha  wing  and  the  tranaduoer  on  the  store.  However,  this 
combination  Is  probably  \maooeptable  because  of  tha  possible  unreliability  of  any  transfer  function 
across  the  store/wlng  Junction.  Therefore  the  study  has  been  with  the  force  and  transducer  on  the 
wing  or  both  of  them  on  the  store. 

5.  RESPLTS  OF  FIOOTEB  COmHOL  SYSTEM  DESIGN  CAICDLATIONS 

5.1  Basic  Aircraft  Characteristics 

Figures  6 and  7 show  ths  basic  aircraft  flutter  oharaoterlstlcs  and  tha  daalgn  speed  frequanoles  | 

and  dampings.  j 

At  tha  design  speeds  | 

I 

Case  1 has  two  Inatabllitlea,  at  4.U2  hs.  (-1.56%)  and  U.73  h».  (-1.08%).  There  Is  low 
damping  in  the  stable  roots  at  10.47  hs.  (1%)  and  35.04  hs.  (1.48%). 

Case  2 has  an  instability  at  6.27  hs.  (-1.87%)  with  very  low  damping  at  8.09  hs.  (0.26%)  and 
low  damping  at  11.8  hs.  (0.98%)  and  24.8  hs.  (1.31%).  { 

Caaa  3 has  two  Instabilities,  at  6.28  hs.  (-1.93%)  and  6.73  (-6.12%)  with  low  A.«r<"g  I 

In  the  stable  root  at  9.57  hs.  (0.88%). 

Case  4 has  an  Instability  at  7.04  hs.  with  a simII  negative  dMipSng  of  -0.42%.  However  the  | 

second  root  at  6.21  hs.  has  an  extremely  small  positive  damping  of  0.06%.  There  is  also  low  daq>lng  I 

at  9.69  hs.  (1.09%),  15.16  hs.  (1.33%),  36.04  hs.  (1.49%)  and  56.31  hs.  (2.09%).  | 

5.2  Chosen  feedback  eystems 

Ths  Hyqulst  Plots  for  all  combinations  of  force  and  transducer  position,  for  all  four  oases 
at  their  design  speeds  were  produced.  Then  the  Hyqulst  Plots  for  the  most  promising  force  and  I 

transducer  positions  were  produced  for  speeds  below  tha  daalgn  apeada.  From  these  ths  chosen  best  | 

nine  combinations  of  force  and  tranaduoer  on  the  wing  are  shown  In  Figure  8.  j 
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5.3  C—  1 rorot  and  lYnaitnor  on  th«  wing  (Fl«ur»  9) 

FiroB  thB  point  of  viow  of  doublo  ■ntloloofcwioo  onolrolonont  of  th*  (-1,0)  point  to  •tnblXiB# 
tho  two  inotnbiUtloB  at  J».U2  and  U.73  tha  only  foroaa  at  55%  chord  to  giwa  raaaonabla  raaulta 
axa  P12  and  F9^. 

For  thaaa  force  poaltlona,  tha  uaalhl  tranaduoor  poaltlona  are  W1,  W12,  W21  andodIO  toocl1l*. 

Howewer,  highar  nodae  (in  particular,  16. U ha.  and  30.1*  ha.)  cannot  bo  filtarad  out 
flrat  ordar  filtara  (althar  alnglo  or  in  oaooada)  ao  thay  can  only  bo  ram>wad  by  doublo  intagration 
and  c*^*^1*^  tha  sign  of  tha  faadbaok  gain# 

The  hast  ramalnlng  of  thaaa  ia  P96/  « 110. 

For  «110,  loada  ware  oonaidarad  at  a range  of  ohordwlaa  poaitlona  along  tho  wing  tip.  Ptob 
thaaa  tha  ohoaan  faadbaoka  are  P90  and  P100.  Aa  tha  gain  and  phaaa  m^glna  are  alnl^  for  thaaa 
two  load  pooitiona  any  force  poaition  batwoan  thaae  two  (including  P96)  can  ba  expected  to  bo  aa 
aooaptablo.  Thle  haa  tha  advantage  of  inaoneitlvity  of  tho  feedback  ayaton  to  centre  of  preeaure 
of  tha  aerodynaalca. 


Whexeaa  tha  gain  nargina  are  adaquata,  tho  phaaa  aarglna  at  tha  dealgn  apood  are  aa  low  aa 
± 30°.  Thoea  phaaa  nargina  Inprova  aa  the  apaad  deoreaaaa.  ij 

5.1j  Caae  2 Force  and  'Pranadnoar  on  the  Wing  (Figure  IO)  j 

P12  axoitoa  tha  unatablo  node  well  relative  to  tha  other  oodea.  For  P12/W85,  tho  velocity  i 

feedback  will  have  to  ba  Intagratod  to  avoid  large  reaponae  at  hi^  frequency  and  a high  paaa 

filter  la  then  required  to  remove  tha  low  frequency  reaponae  and  balance  the  phaaa  marglna.  | 

Aa  the  force  movea  outboard  along  tho  55%  chord  line,  the  l^quiat  Plota  are  apollt  by  higher  j 

n»de  reaponae,  particularly  the  8.O9  hz.  mode,  which  tend  to  diatort  the  flutter  mode  o^le^ 

Thla  ia  notiooable  with  P100/WU3  where  tha  8.09  hz.  roaponao  ia  tha  raaeon  for  the  low  ( t 25  } pnaae 

margina.  Note  that  thie  L.09  hz.  node  haa  only  0.26%  critical  damping  and  Ita  effect  would  be  ^ 

conaiderably  reduced  if  it  wore  given  1%  atruotaral  damping!  j 

Similarly,  but  more  reallatloally,  the  other  choaan  feedback,  P90/  oC  110  loaee  phaae  margin  j 

because  the  laorge  response  In  the  6a22  h»a  stable  mode  distorts  tha  flutter  mode  clrclea  3 

] 

5.5  Caae  3 Force  and  Tranaducar  on  tha  Wing  (Figure  11)  j 

Again,  doublo  encircling  of  the  (-1.0)  point  ia  required.  For  tranaduoara  at  55%  chord  along 
the  wing  and  wing  tip  forcoa  between  25%  and  75%  chord  : moat  plota  are  apollt  althar  by  large 
reaponaea  in  hi^r  modaa,  particular  in  the  ll^tly  daaq^d  9.6  hz.  mode,  or  by  poor  phaae  margina 
due  to  dlatortion  of  tho  6.73  “c8e  reaponae. 

Of  the  ohoaen,  P90/V1  haa  amall  phaae  margina  at  tho  dealgn  apeed  (1:32°)  and  would  benefit  i 

from  phaae  acheduling  with  apeed  becauae  of  the  large  phaae  lag  of  tho  unatablo  modea  going  from  , 

the  design  speed  to  the  flutter  speeda 

For  P100/W85  the  phaae  margina  are  very  amall  ( < * 15°)  ! 

5.6  Caae  li  Force  and  Tranaducer  on  the  Wing  (Figure  12)  ! 

Vhareaa  we  would  expect  the  7.0U  hz.  mode  to  be  eaay  to  atabillee  becauae  it  ia  cn^  Juat 
unatablo  (-0.1*3%  critical  damping),  the  proaonoe  of  the  extremely  lightly  damped  (+0.07%)  mode 
nearby  at  6.22  hz.  oreatea  probleme. 

Becauae  of  this,  and  tha  large  reaponae  in  tho  36  hz,  mode  whan  tho  loading  la  eft  on  the 
wing  tip,  tho  beat  force  poaition  la  forward  on  the  wing  tip. 

Hence,  P90/W1  and  P90/(X.11l*  ware  aolaoted.  Both  give  good  gain  and  phaaa  margina. 

5.7  Effeot  of  Flutter  Control  Svatem  at  Off-Beeign  Conditione 

Figure  13  ahowe  a typical  variation  of  I^quiot  Plot  with  apeed.  The  exaaqple  given  ia  for 
Inboard  elope  ( « 110)  feedback  to  a wing  tip  trailing  edge  control  aurface  (PlOO)  for  caae  1 at 
apeeda  below  tho  dealgn  apeed.  Tho  oooond  and  third  Hyquiet  Plota  ahow  the  aooond  root  lobe 
expanding  aa  the  apeed  reduoea  and  tha  root  beoomea  leee  unatable  (0.1*73  TR)  until  it  haa  a 180  phaae 
ahlft  and  rotataa  olookwiaa  whan  it  beoomea  atable  (0.1*11*  TO).  The  next  two  Myqulot  Plota  ahow  tho 

happening  to  tha  firot  unatablo  root  going  from  0,355  VR  to  O.296  TO. 

Aa  tho  apeed  roduoeo  further,  tho  lobea  contract  (0.237  TO)  aa  damping  Inoreaaaa,  than  expand 
(0.II8TO)  aa  the  aerodynamic  damplnge  reduce.  With  the  reduction  of  tho  aerodynamic  forcee,  the 
lobea  rotate  to%*arde  tho  zero  apeed  phaaing. 

At  very  low  apeed  there  la  generally  a danger  of  feedback  inatability  duo  to  very  low  open 
loop  damping  (otruotural  alone)  and  a I80  phaaa  ahlft  beoauae  the  inertia  force  now  dominatea 
the  control  aurface  force.  Thla  ia  particularly  ao  when  velocity  feedback  ia  uaed. 


L. 


J 


Flgur*  1U  shows  ths  oorxssponding  olossd  loop  frsquenoles  and  dasipings  against  spaad.  Also 
shown  are  the  gain  and  phaaa  margins. 

Figure  15  shows,  for  the  same  faedbaok  law,  ths  olosed  loop  frequencies  and  daagtlngs  agalnat 
feedback  gain  at  the  design  speed.  Notice  there  la  a marl mum  possible  olosed  loop  damping  of  about 
for  the  first  unstable  root  and  about  19^  for  the  second  unstable  root. 

5.8  Summary  of  Hesults  for  Force  and  Transducer  on  the  W^iy 

Figure  8 shows  the  chosen  nine  feedback  systems  with  force  and  transducer  on  the  wing. 

For  six  of  these,  the  force  or  transducer  Is  at  the  Inboard  end  of  the  wing  where  the  modal 
deflections  are  small.  Because  of  this.  It  Is  necessary  to  Include  the  aircraft  rigid  body  and 
fuselage  modes  in  the  analysis. 

Also,  the  rigid  body  modes  need  to  be  Included  because  the  chosen  control  systems  use  dis- 
placement, or  Integrated  dlsplaoemsnt.  faedbaok  and  it  Is  naoessary  to  obaok  that  high  pass  filters  can 
be  used  to  preyent  the  flutter  control  system  from  affeotlng  the  oyerall  aircraft  stability  and  control 
oharaotarlstios . 

No  comnon  system  inside  the  Design  Aims  was  found,  nor  did  it  appear  likely  that  one  would  be 
found  even  If  gain  and  phase  scheduling,  and  notch  filters,  were  allowed. 

However  the  best  oonaton  system  used  a force  at  the  wing  tip  at  2$%  chord  (P90)  together  with 
Inboard  slope  ( 0OIO)  feedback  for  Inboard  stores  and  Inboard  leading  edge  displacement  (V1)  feedback 
for  outboard  stores.  For  this  comaon  system,  the  phase  margin  requirement  would  need  to  be  relaxed 
to  about  1 30  at  5 hz. 

Therefore  an  all  moving  tip  looks  the  most  promising  motivator  but  its  aerodynamics  forces  are 
uncertain  and  need  to  be  Investigated. 

5.9  Control  Forces  on  the  Store  (Figure  16) 

The  feedback  systems  using  a control  force  and  transducer  acting  vertically  on  the  store  front 
are  excellent  because  : 

(1)  they  provide  6db  and  '^60°  margins  at  all  speeds, 

(11)  a slniple  feedback  law  Invariant  with  speed  can  be  used, 

(ill)  a comnon  force/ transducer  position  can  be  used  for  all  four  store/wlng  sweep/Mach  Number 
combinations. 

For  each  case,  the  behaviour  of  the  velocity  transfer  function  Is  like  that  of  a system  of  un- 
coupled, lightly  darned  single  degrees  of  freedom,  where  the  unstable  modes  are  separated  from  the 
stable  modes  by  180°  giving  perfect  stability  margins. 

Also  there  Is  no  response  In  higher  frequency  modes  because  each  store  on  Its  pylon  aots  as  a 
mechanical  filter. 

This  mechanical  filter  and  the  velocity  feedback  remove  ths  possibility  of  closed  loop 
Instability  of  hl^  frequency  modes.  The  velocity  feedback  reduces  Interference  with  the  overall 
aircraft  stability  and  control. 

This  Is  merely  a manifestation  of  the  Importance  of  struotural  damping  on  store  flutter  which 
has  provided,  and  continues  to  provide,  effective  passive  flutter  suppression. 

6.  r.nsT  bbspohSB  CADCUIATIONS 


The  response  to  a 1-ooalne  gust  of  aaq>lltude  25ft./seo.  and  tuned  to  give  earl mum  response  in 
the  flutter  modes  was  oouiputed  at  the  design  speed  at  sea  level  with  the  feedback  system  operating. 

This  showed  that  the  force  that  the  control  surface  must  be  oapable  of  producing  to  make  this 
response  stable  depends  upon  its  position. 


Control  Surface  Position 

Approximate  Bequired  Force 

Wing  Tip 

3,000  lb 

Wing  Boot 

7,000  lb 

Store  Front 

UOO  lb 

Tbsse  results  show  a surprisingly  small  dlfferenoe  between  the  required  forces  at  the  wing  tip 
and  wing  root,  and  an  even  more  surprisingly  small  force  required  on  the  store  oompared  with  the 
force  required  on  ths  wing  tip. 
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This  •oooaaioal  ato^r  baa  raaaalad  tliat  In  oadar  to  paooaad  ftirtliar  wa  aaad  a battar  rap- 
raaaatatloD  of  a «ii«  tip  oontaol  aorfaoa,  oootrol  aaocodpnanloa,  iaartia  and  aetoator. 

Vltbont  thla,  no  ooaiB  faadbaok  w/mUm  ooold  ba  foond  for  a foroa  and  taanaAioar  on  ttaa  wli« 
whloh  aatiafiad  tha  Saai^  Aina  for  all  four  oaaaa.  Ria  dlffleoltlaa  wava  aaialjr  dna  to 

(l)  tvD  Inatabilltiaa  at  Tazj  eloaa  fMqaaoelaa 

and  (il)  wny  low  implagB  In  aoaa  of  tha  atabla  nodaa. 

In  addition,  owarall  alxoraft  nodallinK  ia  aaaantial  ainoa  pieoblana  with  rigid  body  djmaniea 
axa  likaly. 

It  ia  poaaibla  that  tha  aaaxoh  for  a ooaon  qratan  would  ba  aora  aoooaaaftL.  If  thara  waa  a 
battar  rapraaantatloo  of  tha  atmotiixal  ila^li^  and  tha  Saai^i  iiaa  warn  ralazad.  In  partioalar,  it 
would  ba  aaaiar  if  tha  atabilitj  aatgia  raquixananta  wara  radnoad.  Alao,  oaii^  nana  than  on*  trana- 
dnoar  nag  halp. 

Iha  axoallant  raaolta  for  tha  fo*oa  and  tranadnoar  on  tha  atora  aza  intaraatiiw  bat  tha  t if- 
fioalty  in  oaing  aaoh  a ayatan  probably  nokaa  it  l^zaotioal.  Howawar  tha  Mall  foroaa  ra<|alrad  on 
tha  atora  probably  naana  that  it  ia  aaay  to  atvpraaa  tha  flnttor  paaaiwaly  oaii^  a naohanioal  daM>er. 

HPcatwmtTTnp  aw>  amnm  worn 

(1 ) Imraatigata  tha  aarodymnioa  of  an  all  noTliv  wlj^g  tip 

(2)  ttaa  flnttor  oontzol  affioianoiao  of  ttaa  otaoaan  ayatMa  wIim  oootrol  anrfaoa  aarodynaaioa 
and  inHrtia,  and  total  aireraft  aodaa  aza  inolndad.  Safina  ttaa  aotaator  1 nniili  Manta  for  aotiwa 
flnttar  oootrol. 

(3)  Invaatigata  ttaa  offaot  of  ttaa  aotiwa  flnttar  oontzol  ayatan  on  ttaa  cloan  wing  otability.  Safina 
ttaa  aotaator  ragnlrananta  to  avoid  introdneing  oootrol  aorfaoa  flnttar. 

(U)  Inaaatigata  ttaa  affoot  on  aotiwa  flnttar  oontrol  zaqnirananta  oft 

(i)  vazyiiw  ftwl  diatribation  in  tha  wing  and  atora, 

(ii)  vazyiiv  atora  attaotaoant  atifftaaaa  and  inartia  propartiaa 

(ill)  naing  iMdIng  adga  oontrol  aarfaoaa. 

($)  Inaaatigata  tha  naa  of  ooza  than  ona  tranadnaar  in  a flnttor  oontzol  ayatan. 

(6)  Bariaw  ttaa  gain  and  phaoa  nazgin  ragnlrananta  for  a flnttar  oontrol  ayatan. 

COiaJSDC  muvrft 

(1)  Tha  prooadnra  daoorlbad  taara  io  baaad  upon  ttaa  naa  of  l^yqniot  Flota.  Daii^  it  with  a point 
foroa  to  rapzaaant  ttaa  loading  dna  to  a oontzol  oorfaea  on  tho  wing  faoilitataa  an  aoononio 
firot  oaaaaanont  of  tho  optinnn  poaitioo  for  on  aotiva  flattar  control  ayatan  notivator. 

(2)  for  a variabla  awaap  wing  with  tmdazwii^  atorao,  no  ooMon  octiwa  flatter  oontrol  ayatan  naing 
a oootrol  aorfaoa  on  the  wii«  waa  found  irtiioh  gaaa  good  atobilitr  ""’V****  for  foor  oonblnationa 
of  Haoh  aunbar,  wing  awaap  and  atora  ooofigaration.  Tba  baot  cannon  ayatan  found  noed  on  all 
noring  wing  tip. 

Tha  onthora  a^haaiM  that  ttaa  viowa  oa^zaaaad  in  ttaia  p^ior  arc  tbair  own. 

Thla  papor  o<mtalna  raaolta  whloh  warn  aohlaaad  ondor  contract  for  tba  Brltioh  Hinlatzy  of  Safcooa. 
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Fig.2  Design  cases 


CASE  DESIGN  SPEED 
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FREQUENCY 


DESIGN  SPEED 


FLUTTER  SPEED 


0 575VR 


0-296VR 
0 459VR 


4 36  Hi 
4-73 


0 S75VR  0 2S2VR 


6-30  Hi 


0 7 57VR 


0 314  VR 
0 -3UVR 


6-20  Hi 
7 15 


0S05VR 


• 621VR 


7-13  Hi 


VR=  reference  speed 

Fig.6  Basic  aircraft  characteristics 
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Fig.7  Design  speed  frequencies  and  damping 
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FEEDBACK  LAW:  P=Kf(s)W 

P = POINT  FORCE 
K = FEEDBACK  GAIN 
f(s)=  FEEDBACK  FILTER 
W = STRUCTURAL  VELOCITY 


Fig.8  Chosen  best  flutter  control  systems  with  control  surface  on  the  wing 
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INTRODUCTION 

Calculations  of  transfer  functions  to  predict  the  structural  response  of  elastic 
aircraft  to  gust,  manoeuvres,  01  landing  impact  excitation  are  becoming  increasingly 
elaborate,  and  yet  produce  results  which  show  only  limited  agreement  with  experimen- 
tal results  in  return  for  extensive  and  time-consuming  computation.  The  theoretical 
principles  and  their  linearized  versions  required  for  these  computations  are  well 
known.  In  order  to  reduce  the  discrepancies  between  theoretical  and  experimental  re- 
sults, input  data.  Interim  results,  and  part  of  the  final  results  of  the  calcula- 
tions are  compared  with  corresponding  experimental  results  and  corrected  to  agree 
with  those.  Ref.  [11  , [21  , [31  • 


To  get  reasonable  results  for  gust  response  calculations  it  is  necessary  to  intro- 
duce the  elastic  aircraft  behaviour  as  well  as  the  command  and  stability  augmentation 
system  (CSAS)  into  the  mathematical  model.  Once  these  techniques  are  established  they 
can  be  used  for  any  kind  of  control  law  including  CCV  technologies. 


In  this  paper  it  is  demonstrated  how  calculation  results  are  influenced  by  using 
aerodynamic  interference  air  forces.  In  the  second  part  the  influence  of  the  CSAS  is 
presented.  Considering  that  all  important  elastic  degrees  of  freedom  are  used  to- 
gether with  a CSAS-model  and  unsteady  air  forces  are  calculated  with  interference  in- 
fluence the  analysis  becomes  very  expensive  and  computer-time-consuming.  To  reduce 
this  heavy  work  load  it  is  necessary  to  find  out  negligible  parameters  and  also  high- 
light the  influence  of  important  ones. 


The  figures  in  this  paper  show  that  the  influence  of  the  CSAS  on  the  dynamic  re- 
sponse is  of  prime  Interest  and  often  exceeds  the  Influence  of  the  elastic  structure 
by  far. 


It  is  also  shown  that  the  unsteady  aerodynamic  forces  should  be  determined  with 
three-dimensional  theories  including  interference  and  that  corrections  to  match  the 
steady  derivatives  measured  in  the  wind  tunnel  should  be  made. 


The  impedance  function,  control  loop  transfer  functions  are  highly  nonlinear  due 
to  the  nonlinearities  of  the  hydraulic  actuators.  All  these  functions  must  be  deter- 
mined experimentally  and  Introduced  into  the  elastic  aircraft  equation.  Response 
plots  of  the  total  system  should  be  calculated  and  compared  with  results  of  so-called 
"structural  mode  coupling  tests".  If  correlation  is  good  a major  part  for  the  inves- 
tigation of  structural  response  of  the  aircraft  due  to  various  input  functions  is 
verified. 


SYMBOLS 


A(cjl 


defined  in  equation  (3) 


13 

H.(cd) 

Hg(cj) 


matrix  of  generalized  structural  damping 


vector  of  frequency  response  function  relating 
the  response  of  generalized  coordinates  to  a 
sinusoidal  turbulence 


vector  of  frequency  response  function  relating 
the  aircraft  response  to  a sinusoidal  turbulence 


I 

k 

K 

M 


nx 

Jpx 

Ipp 

hq 


x(t) 


xe 

Mco) 

y(t) 

Y(co) 

P 

0x1“  1 

^loil 

2m 

2m 


VTT 

reduced  frequency 

generalized  structural  stiffness  matrix 

generalized  mass  matrix 

number  of  zero  crossings 

vector  of  generalized  coordinates 

variance  function  or  mean  square 

matrix  of  actuator  transfer  functions 

matrix  of  control  surface  impedance 

matrix  of  frequency  response  functions  of  CSAS 
relating  the  actuator  input  to  a sinusoidal  dis- 
placement of  the  generalized  coordinates 

matrix  of  frequency  response  functions  of  CSAS 
relating  the  actuator  input  to  a sinusoidal  dis- 
placement of  the  generalized  coordinates 

matrix  of  frequency  response  functions  relating 
the  required  aircraft  response  to  sinusoidal  gen- 
eralized coordinates 

vector  of  output  functions 

actuator  input  displacement 

vector  of  Fourier  transformations  for  x 

vector  of  input  functions 

vector  of  Fourier  transformations  for  ^ 

vector  of  control  surface  deflections 

input  power  spectral  density 

output  power  spectral  density 

vector  of  generalized  forces  due  to  motion 

matrix  of  generalized  forces  due  to  sinusoidal 
displacements  of  generalized  coordinates 


2-3 


vector  of  generalized 

forces 

due 

to 

excitation 

Qr 

vector  of  generalized 

forces 

due 

to 

CSAS 

CJ 

circular  frequency 

2.  THE  EQUATION  OF  MOTIONS  OF  A FREE  ELASTIC  STRUCTURE  WITH  CSAS 


To  describe  the  dynamical  behaviour  of  an  elastic  structure  that  can  move  freely 
in  the  space  the  following  assumption  must  be  made: 

a)  The  elastomechanical  behaviour  can  be  described  by  linear  equations.  This 
assumes  that  the  elastic  deformations  are  small  compared  with  the  dimensions 
and  the  motions  of  a structure  and  it  requires  the  validity  of  the  Hooke's 
Law. 

b)  The  dynamical  behaviour  can  be  described  by  linear  relations.  This  implies 
the  applicability  of  the  small  perturbations  method. 

c)  The  specific  weight  of  a structure  is  invariant  (for  instance  fuel  consump- 
tion is  not  considered). 

d)  The  elastic  structure  having  an  infinite  number  of  degrees  of  freedom  will 
be  replaced  by  a system  with  a limited  number  of  degrees  of  freedom.  For 
this  idealization  it  does  not  matter  whether  the  deformations  are  a finite 
sum  of  assumed  modes  or  whether  they  are  described  by  the  motions  of  a 
finite  number  of  discrete  elements. 


A system  of  linear  equations  can  be  written  in  a matrix  form. 


Mq  + Dg^+Kq= 


(1) 


The  order  of  the  quadratic  matrices  and  the  number  of  column  of  the  vectors  is 
determined  by  the  number  of  rigid  body  and  elastic  modes  which  are  considered.  The 
elastic  modes  used  are  the  primary  vibration  modes  of  the  free-free  elastic  structure. 
This  is  a common  but  not  necessary  procedure.  If  the  orthogonal  coordinate  system 
having  its  origin  in  the  center  of  gravity  is  referred  to  the  principal  axes  of  the 
structure  then  the  matrix  of  generalized  masses  ^ and  generalized  stiffnesses  iL  is 
a diagonal  matrix  because  the  rigid  body  degrees  of  freedom  as  well  as  the  elastic 
degrees  of  freedom  fulfil  the  known  orthogonality  condition.  For  the  elastic  de- 
grees of  freedom  the  elements  of  matrix^  ore  the  generalized  structural  dampings. 

The  generalized  forces  on  the  right  side  of  the  equation  are  arising  from  motions 
from  excitations I^q). 


This  paper  deals  only  with  the  solutions  of  the  equations  in  the  frequency  domain 
Ref.  . Choosing  this  method  has  the  advantage  that  the  important  functions  - actua- 
tor impedance,  CSAS  transfer  functions  and  unsteady  aerodynamic  airforces  - are  avail- 
able in  the  frequency  domain.  The  PSD-Analysis  to  determine  structure  fatique  life  is 
a direct  result  of  this  method.  The  Fourier  transformation  of  equation  (l)  leads  to 

[-0)2  M + iu)Dt  K- Q (0))1  q = Q (to)  (2) 


A((o)£=  Qq(w) 


A(U))  = -0j2M  + i0)D  + K- 0|^{0))  (3) 


or 


with 


The  transfer  function  of  the  elastic  system  can  be  determined  from 


H(0))=  [A(U))]-’  Qp(U)) 


(4) 


1^ 


The  frequency  responses  of  other  parameters  like  the  total  motions,  the  velocities, 
accelerations  and  of  the  forces,  moments  and  stresses  can  be  described  by  the  trans- 
formation 

Ho(w)  = ToH(w)  (5) 


This  transformation  is  usually  replaced  by  other  analysis  methods  for  the  sake 
of  numerical  accuracy |Ref*  [^1 • With  any  kind  of  deterministic  excitation  the  response 
in  the  frequency  domain  can  be  calculated  from 

W = Ho(w)Y((D)  (6) 


where  Yfcj)  is  the  Fourier  transform  of  the  excitation.  The  response  in  the  time 
domain  is  defined  by  the  transformation 

x(t)=  [ Mod)  e'Wt  do)  (7) 


For  stochastic  excitation  the  power  spectral  density  of  the  response  is  obtained 

from 

^x(w)=  $y(U))  (8) 

where  is  the  power  spectral  density  of  the  excitation.  The  variance  or  mean 

square  values  can  be  determined  from 


5x  = |0x(w)da)  (9) 


and  the  number  of  positive  zero  crossings  from 


Noi  = i (a)2^i(a))doj  (10) 

^xi  0 


In  Fig.  1 the  control  system  is  shown.  The  aircraft  is  the  system  to  be  controlled 
Responses  of  this  system  are  the  generalized  coordinates 


The  system  is  excited  by  gusts  or  buffeting  (changes  of  the  airplane  angle  of 
attack)  or  by  raanouvres  (signals  to  the  actuators  of  control  surfaces).  In  addition 
to  the  external  excitations  the  control  system  generates  input  signals  to  the  actua- 
tors. The  generalized  excitation  force  of  equation  (2)  can  therefore  be  partitioned 
into  the  external  excitation  force  and  into  the  force  produced  by  the  control  system. 


A((i))q  = Qq((1))+  Qp  (OJ) 


(11) 


The  Impedances  of  the  hydraulic  actuators  are  strongly  frequency  dependent.  Cal- 
culation of  frequencies  and  primary  modes  considering  the  actuator  stiffnesses  leads 
to  an  Eigenvalue  problem  with  complex  frequency  dependent  stiffnesses.  It  is  not  nec- 
essary to  solve  this  problem  if  the  normal  modes  are  determined  with  infinitely  stiff 
control  surface  attachments.  The  elastic  deformation  of  the  control  surfaces  must  be 
contained  in  these  modes.  The  degrees  of  freedom  missing  in  these  modes  can  be  intro- 
duced later  by  additional  attachment  modes  (secondary  modes).  It  should  be  borne  in 
mind  that  the  orthogonality  condition  is  no  more  fullfllled  using  primary  modes  and 
secondary  modes.  For  this  reason  the  mass  matrix  will  have  off  diagonal  elements. 

The  amplitudes  of  generalized  coordinates  for  the  primary  modes  will  be  described 
by  the  vector  ^ and  for  the  control  surface  rotations  (secondary  modes)  by  the  vector|3. 
Equation  (ll)  can  now  be  written  down.  ~ 

-qq  ~q6  S ®qD  o 

- + (12) 

Apq  AppJ 

From  the  principle  of  the  virtual-work  it  follows  that  the  components  of  the  am- 
plitude vectors  0.^0  are  the  moments  of  the  control  surfaces  and  these  must  be  zero 
because  the  primary  modes  contain  no  control  surface  motions.  The  vector  of  the  ac- 
tuator moments  shall  be  described  by  the  following  relation 


~PR  “ (13) 

It  is  assumed  that  the  dependence  of  the  control  surface  moments  from  the  actuator 
inputs  can  be  described  by  linear  transfer  functions.  The  coefficients  of  the  matrix 
T(j|j  are  the  control  surface  impedance  which  are  the  reaction  moments  of  the  actuators 
due  to  control  surface  motion.  The  second  product  of  equation  (13)  is  considered  to  be 
a part  o{  Qnn  • According  to  equation  (ll)  T(j/j  is  the  matrix  of  control  surface  stiff- 
nesses and  therefore  an  additive  term  of  matrix  • 

Duo  to  the  motion  of  the  airplane, input  signals  to  the  actuators  of  the  control 
surfaces  are  generated.  Since  the  motions  of  the  elastic  structure  were  developed  in  a 
finite  number  of  rigid  body  modes  and  primary  modes  the  input  signals  can  be  described 
as  a function  of  generalized  coordinates 

Jie=lRq(‘*>)a+5qa 

The  coefficients  of  the  matrix  Tf,^ are  the  transfer  functions  of  the  control  system 
including  the  sensors  multiplied  with  the  amplitudes  of  the  primary  modes  at  the  sensor 
stations.  The  coefficients  of  the  matrix  Tg^  are  the  input  sigrals  which  for  instance 
arise  from  mechanical  Inputs  to  the  actuators  due  to  different  deformations  of  the 
structure  and  the  control  rods  and  are  therefore  determined  by  the  vibration  modes. 

x^  therefore  only  describes  the  actuator  input  signals  introduced  by  the  control 
system.  Reactions  of  the  actuator  Inputs  due  to  control  surface  moments  are  not  con- 
tained in  Xg.  How  to  determine  the  matrix  coefficients  of  equation  (13)  and  equation 
(ik)  experimentally  is  described  in  Ref.  (5)  , [61,(71  or  the  respectxve  literature.  If 
the  control  surfaces  are  not  dependent  upon  each  other  then  only  the  main  diagonal 
elements  of  the  matrixes  of  equation  (13)  and  equation  (l4)  are  not  zero.  Equation  (13) 
and  equation  (l4)  inserted  in  equation  (12)  give 

-qq 

Apq-]px<jRq^Jsq) 
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3.  INTERFERING  UNSTEADY  AERODYNAMIC  FORCES 

The  effect  of  interference  between  control  surfaces  of  an  aircraft  in  steady  flow 
is  well  known  and  is  introduced  in  flight  mechanics  by  theoretical  estimates  and  meas- 
urements. Neglecting  these  terms  would  lead  to  wrong  flight  mechanical  results. 

The  aerodynamic  interference  is  also  important  for  unsteady  flow.  To  estimate  its 
impact  is  much  more  difficult  considering  that  no  corresponding  unsteady  measurements 
are  available  to  aid  these  est imates . After  a few  flutter  accidents  especially  on  air- 
planes with  T-tails  the  interference  between  harmonically  oscillating  surfaces  was 
analytically  investigated.  The  first  topic  which  was  investigated  was  the  interference 
between  tailplane  and  fin  especially  for  the  T-tail  , Ref . C8]  , (9^.  After  instabilities 
created  between  tail  and  wing  interference  became  known  the  theoretical  investigation 
of  Interference  between  all  surfaces  was  considered, Ref . [lOl  . A surface  method  was 
developed  by  B.Laschka  which  was  used  to  create  the  results  shown  here, Ref.  Ill]  . Ex- 
perimentally the  method  was  substantiated  by  J. Becker  in  Ref.  112]  . The  unsteady  rudder 
air  forces  were  calculated  with  the  so-called  equivalent  slope  method, Ref . I 13J • Good 
correlation  of  these  analytical  air  forces  with  experiment  was  presented  in  Ref.  [l41 . 

In  earlier  gust  calculations  without  the  complete  consideration  of  the  theoretical 
interference  the  downwash  of  the  wing  on  the  tailplane  was  estimated  and  considered  by 
introducing  a dowiiwash  factor  for  harmonical  motion.  This  downwash  (due  to  motion) 
depends  upon  the  angle  of  attack  of  the  wing  and  can  be  easily  introduced  for  the  heave 
motion  (01  a constant)  but  for  an  arbitrary  oc-distribution  on  the  wing  a useful  estimate 
is  extremely  difficult. 

Fig.  2 shows  analytical  air  force  distributions  in  spanwise  direction  with  and 
without  interference  for  the  symmetrical  case  (heave-motion)  on  the  wing  and  on  the 
tailplane.  For  the  wing  there  is  practically  no  interference  effect  whereas  on  the 
tailplane  there  is  a considerable  reduction  of  the  air  force  and  also  a reversed  phase. 
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AIR  FORCE  DISTRI0U1  ION  ON  THE  WING  DUE  TO  HEAVE  MOTION 


AIR  FORCE  DISTRIBUTION  ON  THE  TAILPLANE  DUE  TO  HE  AVE  MOTION 


WITH  AERODYNAMIC  INTERFERENCE 
WITHOUT  AERODYNAMIC  INTERFERENCE 


FIG.  2 AIR  FORCE  DISTRIBUTIONS  ALONG  THE  SPAN  DUE  TO 
HEAVE  MOTION  (k  = 0.05) 


Fig.  3 ahowa  correaponding  reaulta  for  antiaymmetrical  motion  of  fin  and  tail- 
plane.  The  air  force  diatributiona  on  the  fin  ahow  in  additon  to  the  reaulta  with  and 
without  Interference  the  reaulta  of  the  calculation  with  a fin  which  ia  aymmetrically 
extended.  It  ia  ahown  that  thia  method  almcat  givea  the  same  results  for  the  fin  as 
the  consideration  of  interference.  For  the  tailplane  there  are  considerable  differences 

The  last  section  of  this  chapter  deals  with  the  correction  of  analytically 
determined  air  forces  by  measured  values.  Between  calculated  and  experimental  air 
force  distributions  and  their  corresponding  derivatives  there  are  differences  which 
cannot  be  neglected.  There  are  various  reasons  of  this  differences  like  boundary 
layer  thickness  effects,  fuselage  Interference  etc.  In  general  there  are  only  test 
results  for  steady  flow  available,  mainly  derivatives  for  certain  parts  or  for  the 
total  airplane. 


AIR  FORCE  DISTRIBUTION  ON  THE  FIN  DUE  TO  LATERAL  TRANSLATION 


AIR  FORCE  DISTRIBUTION  ON  THE  FIN  DUE  TO  YAW  MOTION 


AIR  FORCE  DISTRIBUTION  ON  THE  TAILPLANE  DUE  TO  ROLL  MOTION 
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FIG.  3 AIR  FORCE  DISTRIBUTIONS  ON  FIN  AND  TAILPLANE 


Unsteady  pressure  distribution  measurements  especially  for  elastic  motions  are 
very  seldomly  available  to  that  a correction  of  theoretical  air  forces  due  to  elastic 
motions  is  practically  impossible. 

For  the  rigid  body  modes  a correction  of  the  theoretical  air  force  distribution 
and/or  derivatives  is  possible  and  necessary  because  usually  test  results  are  avail- 
able and  the  loads  are  mainly  determined  by  rigid  body  modes.  To  adjust  the  theoreti- 
cal derivatives  to  test  results  the  etnalytically  determined  downwash  is  modified  lo- 
cally or  globally  depending  whether  a distribution  measurement  or  measured  derivatives 
for  the  rigid  body  modes  are  available.  For  that  reason  it  is  necessary  to  find  a 
correspondance  between  the  stationary  flight  mechanical  and  the  unsteady  aeroelastical 
method.  For  the  stationary  case  (k — 0)  these  relations  can  be  found.  It  should  be 
noted  here  that  the  two-dimensional  unsteady  theory  is  not  applicable  for  some  deriv- 
atives because  there  stationary  values (k— ^ 0)  become  singular  whereas  the  three-dimen- 
sional theory  leads  to  useful  results  which  cem  be  used  to  complete  the  stationary 
derivatives. 

In  Fig.  4 a comparison  of  unsteady  gust  air  forces  on  the  tailplane  and  on  the 
total  aircraft  with  and  without  interference  is  given.  It  is  also  shown  how  a matching 
of  the  stationary  values  (k  « 0)  influences  the  results.  From  this  picture  one  can 
deduce  that  the  calculation  with  interference  already  gives  a good  match  of  the  station- 
ary values.  Not  considering  the  interference  aerodynamics  produces  much  bigger  gust 
forces  on  the  tailplane. 
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GENERALIZED  GUST  FORCE  ON  THE  TAILELANE 
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GENERALIZED  CUST  FORCE  ON  THE  TOTAL  AIRCR  AFT 
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FIG.  4 GENERALIZED  GUST  FORCE 
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4.  IMPACT  OF  AERODYNAMIC  INTERFERENCE  ON  THE  DYNAMIC  RESPONSE 

The  dynamic  response  calculations  were  made  using  the  unsteady  airforces  described 
as  follows: 


a) 

not  interfering  (corrected) 

b) 

interfering 

(corrected) 

c) 

interfering 

(not  corrected) 

In  Fig.  5 the  root  shear  force  on  tailplane  and  wing  is  presented  for  discrete 
gusts  with  various  length.  Here  it  is  shown  that  the  influence  on  the  wing  is  small 
compared  with  that  on  the  tailplane. 


WITH  INTERFERENCE  I AIRFORCESCORRirTEK  WITH 

f STEADY  WIND  TUNNEL  Dl  RIVATIVI  S 
WITHOUT  INTEREERENCtJ 


FIG.  5 SHEAR  FORCE  AND  BENDING  MOMENT  AS  A FUNCTION 
OF  VERTICAL  GUST  LENGTH  (ELASTIC  AIRCRAFT) 
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DYNAMIC  RESPONSE  CALCULATIONS  WITH  AND  WITHOUT  CSAS 


This  chapter  deals  with  the  impact  of  the  Comnand  and  Stability  Augmentation  System 
(CSAS)  on  the  gust  response. 

The  airplane  we  are  discussing  here  is  the  first  operational  sweepable  wing  aircraft 
featuring  a triplex  analogue  fly-by-wire  control  system,  mechanical  emergency  control 
and  automatic  stabilization.  The  primary  flight  control  system  provides  pitch,  roll  and 
yaw  control  by  means  of  an  all  moving  tailplane  (taileron),  a conventional  rudder,  and 
wing  mounted  spoilers.  The  tailerons  operate  in  phase  for  pitch  and  differentially  for 
roll  control.  The  spoilers  give  augmented  roll  control  at  unswept  and  Intermediate  wing 
positions  at  low  speeds,  and  also  act  as  lift  dumper  after  touch  down. 

The  flying  control  surfaces  are  actuated  by  tandem  hydraulic  Jacks.  Two  completely 
separate  and  independent  hydraulic  systems  provide  fully  duplicated  power  for  the  pri- 
mary and  secondary  flying  controls.  The  control  stick  and  rudder  pedal  movement  is 
picked  off  by  triplex  electrical  position  sensors  which  generate  the  command  signals  to 
the  CSAS. 

The  main  sensors  for  feeding  back  the  aircraft  motion  are  rate  gyros.  Both,  the 
command  signals  and  the  feedback  signals  are  passed  through  appropriate  gain  schedulers 
and  filters  before  they  are  fed  to  the  control  surface  actuators.  As  main  scheduling 
parameters  dynamic  pressure,  wing  sweep  and  a flap  switch  signal  are  used.  For  im- 
proving the  rolling  characteristics,  especially  turn  entry  and  turn  exit,  a roll  to 
rudder  Interconnect  is  implemented  and  scheduled  with  a pitch  stick  increment  signal 
in  addition  to  the  forementioned  scheduling  parameters. 

In  pitch  axis  signals  are  added,  which  compensate  for  the  pitching  moments  gen- 
erated by  airbrake  and  flap  deployment. 

A block  diagram  of  the  main  elements  of  the  CSAS  and  the  CAS  is  shown  in  Fig,  6. 


FJTCH  STICK 


FIG.  6 CSAS-CONTROL  SYSTEM 


In  Fig.  7 the  wing  root  bending  moment  for  vortical  ( 1-cos)  gust  is  depicted. 
These  bending  moments  are  presented  for  the  rigid  and  elastic  mathematical  model  with 
and  without  CSAS.  For  the  elastic  mathematical  model  two  rigid  body  modes  - heave  and 
pitch  - together  with  15  elastic  modes  of  the  free-free  structure  were  used.  The  CSAS 
was  considered  by  three  degrees  of  freedom  of  the  all  moveable  tail. 
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WING  ROOT  BENDING  MOMENT  WimC  root  BENDING  MOMENT  **1NG  ROOT  BENDING  MOMENT  FOR  VARIOUS 

GUST  LENGTH 

O RIGID  a/C  without  C SAS  A ELASTIC  A)C  WITHOUT  CSAS 

• RIGID  A/C  WITH  CSAS  A ELASTIC  A/C  WITH  CSAS 


FIG.  7 WING  ROOT  BENDING  MOMENT  AS  A FUNCTION 
OF  VERTICAL  GUST  LENGTH 


One  can  see  that  the  considerable  influence  of  the  CSAS  on  the  wing  root  bending 
moment  is  only  apparent  for  the  low  frequency  regime.  This  means  that  in  the  frequency 
regime  of  the  elastic  modes  the  influence  of  the  CSAS  is  negligible  whereas  it  must  be 
introduced  in  the  low  frequency  regime  near  the  short  period  mode  frequency  to  get 
realistic  structural  loads  due  to  gusts. 

Fig.  8 shows  the  bending  moment,  the  shear  force  and  the  torsional  moment  on  the 
fin  root  due  to  a lateral  (l-cos)  gust.  These  results  were  also  derived  with  a rigid 
and  an  elastic  mathematical  model  with  and  without  CSAS.  For  the  euitisymmetrical 
mathematical  model  three  rigid  body  degrees  of  freedom  - side  translation,  roll  and 
yaw  - together  with  l8  elastic  modes  were  used.  The  CSAS  was  considered  with  four 
modes.  Fig.  8 shows  similar  results  for  the  fin  root  as  Fig,  7 for  the  wing  as  far 
as  the  influence  of  the  CSAS  is  concerned. 


FIN  ROOT  SHEAR  FORCE  FIN  ROOT  BFNDING  MOMENT  HN  ROOT  TORSIONAL  MOMENT 
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FIG.  8 FIN  ROOT  LOADS  AS  A FUNCTION  OF  LATERAL  GUST  LENGTH 


TRANSFER  FUNCTION  OF  THE  ROOT  FIN  SHEAR  FORCE  TRANSFER  FUNCTION  OF  THE  LEADING  EDGE  TIP  ACCELERATION 


FIG.  9 FIN  TRANSFER  FUNCTIONS  FOR  LATERAL  GUST 


It  can  be  seen  from  this  picture,  that  there  is  a big  Influence  of  CSAS  on  Fin 
Shear  Force  is  considerable.  This  can  be  explained  by  the  fact  that  the  gust  spectrum 
(integral  of  equation  9)  is  mainly  determined  by  the  influence  of  the  low  frequency 
amplitudes . 

In  Fig.  10  the  frequency  of  exceedance  of  the  wing  root  bending  moment  due  to 
vertical  stochastic  gust  excitation  and  the  shear  force  on  the  fin  root  due  to  lateral 
stochastic  gust  excitation  is  shown.  These  results  were  calculated  with  the  mathemati- 
cal model  described  for  Fig.  7 and  8.  The  constants  used  for  this  calculation  were 
taken  from  the  US-Airforce  Military  Specification  886l  A (May  I960).  The  influence  of 
the  CSAS  on  the  fin  root  shear  force  is  considerable. 


FIG.  10  FREQUENCY  OF  EXCEEDANCE  OF  WING  ROOT  BENDING  MOMENT 
AND  FIN  ROOT  SHEAR  FORCE  FOR  GUST  EXCITATION 


The  last  part  of  this  paper  deals  with  the  Influence  of  CSAS  on  response  calcula- 
tions due  to  manoeuvre. In  order  to  sort  out  the  different  Influences  three  different 
kinds  of  tailplane  movements  were  considered: 

a)  the  theoretical  trapezoidal  tailplane  movement 

b)  the  trapezoidal  tailplane  movement  multiplied  with  the  actuator  functions 

c)  the  trapezoidal  tailplane  movement  multiplied  with  the  transfer  function  of  the 
command  augmentation  system  (CAS)  and  the  actuator. 

Those  three  tailplane  motions  are  presented  in  Fig.  11. 


FIG.  11  TAILPLANE  ANGLE  INPUT  SIGNAL 


These  three  motions  were  introduced  as  tailplane  manoeuvre  input  into  the  rigid 
and  elastic  mathematical  model.  In  Fig.  12  the  vertical  acceleration  of  the  center 
of  gravity  is  depicted. 
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FIG.  12  A/C  CENTER  OF  GRAVITY  VERTICAL  ACCELERATION 
DUE  TO  MANOEUVRE  (ELASTIC  A/C) 
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TAILPLANE  ANGLE 


The  corresponding  vertical  accelerations  on  the  wing  tip  are  presented  in  Fig.  13 
and  the  tailplane  root  shear  forces  in  Fig. 
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FIG.  13  WING  TIP  VERTICAL  ACCELERATION  DUE  TO 

MANOEUVRE  (ELASTIC  A/C) 
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TAILPLANE  ROOT  SHEAR  FORCE  FOR  RIGIO  A/C  TAILPLANE  ROOT  SHEAR  FORCE  FOR  ELASTIC  A/C 


FIG.  14  TAILPLANE  ROOT  SHEAR  FORCE  DUE  TO  MANOEUVRE 


Together  with  the  responses  of  the  uncontrolled  aircraft  for  the  three  different 
tailplane  modes  in  Fig.  12  to  I'l  the  responses  with  the  full  CSAS  are  also  plotted. 

The  mathematical  model  used,  was  described  for  Fig,  7.  Figs. 12  -Ik  show  that  the 
consideration  of  the  actuator  transfer  function  reduces  the  vertical  acceleration  as 
well  as  the  root  shear  force  on  the  tailplane  more  than  10  %.  The  introduction  of 
the  CAS  gives  only  a time  shift  of  the  maxima  of  the  plotted  accelerations  whereas 
the  tailplane  root  shear  force  is  reduced  more  than  20  The  mathematical  model  with 
the  full  CSAS  reduces  the  structural  loads  about  25  % - 35  % compared  with  the  re- 
sponses without  CSAS  and  trapezoidal  excitation. 

A comparison  of  the  responses  for  the  elastic  and  rigid  mathematical  model  shows 
that  the  influence  of  structural  elasticity  is  small.  A reasonable  explanation  for 
this  small  Influence  is  the  low  excitation  frequency  having  a big  frequency  distance 
.to  the  first  elastic  mode.  The  responses  for  the  trapezoidal  excitation  show  for  the 
elastic  model  that  higher  modes  of  the  elastic  structure  are  excited.  This  can  be 
explained  by  the  fact  that  this  excitation  cannot  be  steadily  differentiated  which 
does  not  occur  for  the  other  forms  of  excitation.  Such  an  excitation  of  higher  modes 
is  therefore  a result  of  a theoretical  model  of  a tailplane  motion  which  does  not 
occur  in  flight. 

An  Intensive  study  of  all  responses  show  that  the  CSAS  must  be  very  carefully 
introduced  into  the  mathematical  model  to  determine  structural  loads  due  to  manoeuvres 
whereas  the  elastic  behaviour  can  be  neglected  if  there  is  enough  frequency  space 
between  the  main  frequency  of  the  excisting  manoeuvre  and  the  frequency  of  the  first 
elastic  mode. 
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ACTIVE  FLUTTER  SUPPRESSION  ON  AN  AIRPLANE  WITH 
WING  MOUNTED  EXTERNAL  STORES 


H.  HONLINGER 

MESSERSCHMITT-BOLKOW-BLOHM  GmbH. 
Unternehmensbereich  Flugzeuge 
8 Miinchen  80,  Postfach  801160, -W. Germany 


1.  INTRODUCTION 

After  a successful  application  of  active  flutter  suppression  on  wing  store  and 
empennage  flutter  problems  [ l]  ,[2],[3^  in  wind  tunnels  an  extension  of  this  technology 
to  a full  scale  airplane  was  considered  to  be  rewarding.  The  effort  was  focused  on  the 
flight  test  of  a wing/store  flutter  suppression  system  (FSS)  with  store  mounted  vanes. 

This  program  started  in  1975  with  the  design  of  the  system  and  the  instrumentation 
of  a FIAT  G 91/T3  as  flying  test  bed.  Flight  test  was  finished  in  Febr.  1977. 

Test  objectives  of  this  study  were: 

. Provide  first  flight  experience  with  FSS  on  external  stores 

. Substantiate  and  demonstrate  a new  method  for  flight  flutter  testing  wing  mounted 
external  stores  by  use  of  this  FSS  (Automatic  Mode  Excitation  System  - AMES) 

This  work  was  carried  out  under  the  ZTL  research  contract  for  the  German  Ministry 
of  Defense  by  the  MBB  structural  dynamics  group.  The  flight  test  was  performed  by 
BWB  - LG  IV  8 - and  E 6l  at  the  German  Airforce  Test  Centre  at  Manching. 
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Two  modes,  shown  in  Fig.  2,  wing  bending  at  11.39  Hz  and  store  pitch  wxth  wing 
torsion  at  l4.l5  Hz  are  causing  the  flutter  case  which  we  intended  to  suppress. 


MODE  2 1139  Hz  MODE  3 K.15  Hz 


FIG.  2 VIBRATION  MODE  SHAPES 


2.2  Control  System 

The  damping  of  vibration  or  flutter  modes  of  wing  mounted  stores  by  oscillating 
vanes  creating  aerodynamic  forces  was  tested  in  the  wind  tunnel  to  be  a very  effective 
way  and  should  now  be  applied  on  an  airplane.  The  block  diagram  in  Fig.  3 shows  the 
control  loop  of  the  airplane  with  the  FSS. 


The  control  signal  is  produced  by  two  accelerometers  located  in  the  forward  and 
rear  part  of  the  flutter  tank.  This  signals  were  added  and  integrated  to  give  the 
pitch  angle  of  velocity  of  the  flutter  store  which  is  u;-ed  as  input  signal  to  the 
hydraulic  actuators  driving  small  vanes  attached  at  the  forward  part  of  the  store. 

The  vanes  oscillate  in  such  a way  that  the  generated  airforces  counteract  the  pitch 
motion  of  the  store  like  a velocity  proportional  damper.  The  steady  aerodynamic 
derivatives  of  the  A/C  were  not  influenced  by  these  vanes  and  therefore  no  change 
in  the  flight  mechanical  behaviour  was  expected.  For  optimization  of  the  phase  between 
control  signal  and  angle  of  attack  of  the  vanes  a manual  phase  shift  device  was  pro- 
vided in  the  control  loop. 


2.3  Analysis  of  the  Control  System 

An  analysis  of  the  airplane  with  control  system  was  carried  out  with  a computer 
program  deyeloped  at  MBB[i33*  This  analysis  predicted  that  the  system  could  be  suc- 
cessfully applied.  An  explanation  of  the  applied  Nyquist  criteria  is  shown  in  Fig.  k. 
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FIG.  4 NYQUIST  STABILITY  CRITERIA 


In  Fig.  5 the  Nyquist  diagram  for  v ■ 5OO  kts  for  the  airplane  with  FSS  is  de- 
picted. It  shows  that  the  mode  to  be  controlled  is  stable.  With  the  automatic  mode 
excitation  system  (AMES)  this  mode  becomes  unstable.  For  y • 6OO  kts  the  system  is 
fluttering  and  becomes  stable  only  with  the  FSS. 


FIG.  S NYQUIST  DIAGRAM  FOR  THE  FLUTTER  SUPPRESSION  (FSS)  AND 
MODE  EXCITATION  SYSTEM  (AMES) 


2. ^ Technical  Realization  of  the  FSS 


As  mentioned  before  the  FSS  had  to  be  implemented  into  the  two  ballasted  520  Itr. 
tanks  of  the  G 91.  For  the  sake  of  redundancy  one  tank  on  each  wing  was  installed. 

Each  tank  worked  Independently  and  only  needed  electrical  power  from  the  aircraft. 

Fig.  6 shows  a sketch  of  the  installation  of  the  FSS  in  the  520  Itr.  tank.  A 
frame  was  Installed  in  the  forward  part  of  the  tank  to  carry  the  vanes,  the  electro- 
hydraulic  power  package  and  control  electronics.  The  ballast  weight  was  clamped  in 
the  centre  part  of  the  tank. 


ACTUATOR  CONTROL 


FIG.  6 EXTERNAL  TANK  WITH  FSS  (FLUTTER  TANK) 


The  vanes  were  designed  for  transsonic  flow  (angle  of  leading  edge  30°,  axis  of 
rotation  at  35  % root  chord). 

Special  fast  vane  actuators  with  an  mechanical  fixing  device  were  developed.  In 
Fig.  7 the  transfer  function  of  the  actuator  with  vanes  is  given.  The  Bode  plot  shows 
ideal  behaviour  in  the  required  frequency  range  (up  to  Ik  Hz).  There  is  only  a small 
phase  shift  and  the  amplitudes  remain  constant  up  to  15  Hz.  The  max.  torque  the  actu- 
ators produce  is  2700  Ncm.  The  max.  possible  angle  of  attack  of  the  vanes  is  +_  10°. 
The  weight  of  the  whole  FSS  in  this  test  status  is  approx.  21  kg. 


FIG.  7 ACTUATOR  TRANSFER  FUNCTION 


FIG.  8 FLUTTER  TANK  ON  THE  AIRPLANE 


3.  FLIGHT  TEST  INSTRUMENTATION 

To  perform  the  flight  test  of  the  FSS  and  to  evaluate  new  methods  for  flight 
flutter  tests  special  flight  test  equipment  was  installed  in  the  airplane. 


A frequency  sweep  generator  was  fitted  which  gives  specified  inputs  into  the 
control  system  (variable  frequency  sine  wave).  This  input  can  be  given  in  the  open 
and  closed  loop  mode  as  well. 
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The  acceleration  pickup  locations  are  shovm  in  Fig.  9.  17  parameters  could  be 
recorded  on  tape  and  also  telemetered  for  ground  quick  look  Inspection.  The  pilot 
was  able  to  manually  control  the  FSS  with  two  control  panels. 


Fig.  10  shows  the  main  control  panel  in  front  of  the  pilot.  The  second  control 
panel  was  used  to  switch  the  hydraulic  power  on  and  off.  6 modes  of  operation  of  the 
flutter  tanks  could  be  selected  by  the  pilot.  Also  operation  of  one  side  FSS  was 
possible . 


FIG.  10  MAIN  CONTROL  PANEL 


In  the  open  loop  mode  four  different  angles  of  attack  of  the  vanes  ranging  from 
1°  to  <1°  could  be  choosen  by  the  pilot.  A warning  light  indicated  failure  or  automatic 
cutoff  of  the  system. 

All  operation  modes  could  be  startet  and  cut  off  by  the  pilot  pressing  the  trigger 
button  on  the  stick.  An  emergency  switch  was  also  installed  in  the  stick. 


The  operation  modes  of  the  flutter  tanks  were: 


1. 

2. 

3- 

4. 

5. 

6. 


Open  loop  mode 

(mode  excitation  by  frequency  sweep  inputs) 
Automatic  flutter  suppression  (FSS) 
Automatic  mode  excitation  (AMES) 


FSS  with  additional  frequency  sweep  input 

L/H  flutter  tank  mode  1 (open  loop  excitation) 
R/H  flutter  tank  mode  2 (FSS) 

L/H  flutter  tank  mode  3 (AMES) 

H/H  flutter  tank  mode  2 (FSS) 


Safety  Installations 

In  order  to  avoid  hazardous  flight  conditions  special  safety  provisions  were  made. 
Both  systems  work  independently  and  each  of  them  is  able  to  suppress  flutter  of  the 
airplane  up  to  a defined  velocity. 

If  a hydraulic  or  electronic  failure  occurs  the  vane  actuator  is  blocked  mechani- 
cally and  the  angle  of  attack  of  the  vanes  becomes  zero.  An  automatic  cutoff  device 
was  Installed  which  switches  off  the  system  in  the  open  loop  and  automatic  flutter 
excitation  mode  when  a preselected  acceleration  limit  on  the  flutter  teuiks  or  on  the 
wing  tips  is  exceeded.  In  order  to  be  shock  insensitive  this  system  was  integrating 
amplitudes.  In  the  utmost  emergency  case  the  flutter  tanks  could  be  jettisoned. 


k.  GROUND  TESTS  ON  THE  AIRPLANE 

To  check  the  stability  of  the  control  loop  of  airplane  and  FSS  a structural 
mode  coupling  test  on  ground  was  performed.  The  test  showed  no  detrimental  coupling 
of  the  FSS  with  the  vibration  modes  of  the  airframe  on  ground. 

5.  FLIGHT  TEST 

5.  1 Performance  of  Flight  Test 

From  the  beginning  of  Nov.  1976  till  end  of  Febr.  1977  l8  flights  were  performed. 
17  parameters  were  monitored  at  the  telemetry  station  during  the  flights. 

After  each  flight  the  test  data  were  analyzed  with  an  HP  5^51  Fourier  Analyzer. 
Special  programs  were  developed  to  calculate  Nyquist  diagrams  and  evaluate  damping  from 
the  open  loop  tests  with  frequency  sweep  excitation. 
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5.2  Open  Loop  Teats  in  Flight 

In  order  to  check  the  stability  and  to  optimize  the  phase  of  the  FSS  open  loop 
tests  in  flight  at  various  airspeeds  were  performed  at  the  beginning  of  the  flight 
test.  Frequency  sweeps  were  fed  into  the  control  system.  The  frequency  increased 
according  to  a logarithmic  law  from  5 Hz  to  25  Hz  and  then  decreased  to  5 Hz  within 
120  seconds. 

Fig.  11  shows  the  response  of  the  flutter  tank  to  the  frequency  sweep  excitation 
through  the  vanes.  A good  excitation  of  the  flutter  mode  can  be  seen. 
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FIG.  11  TIME  HISTORY  OF  FREQUENCY  SWEEP  EXCITATION  (OPEN  LOOP) 


In  Fig.  12  two  Nyquist  diagrams  calculated  from  these  frequency  sweeps  at  v » 350 
KIAS  and  v » 4»50  KIAS  are  presented.  The  phase  of  the  flutter  mode  is  not  yet  optimized. 
Comparing  both  plots  the  phase  shift  due  to  the  unsteady  air  forces  can  be  seen  too. 


V=  350  KIAS 
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FIG.  12  NYQUIST  PLOTS  AT  v = 350  KIAS  AND  v = 450  KIAS  (OPEN  LOOP) 
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5.3  Closed  Loop  Testa  of  the  FSS 

After  having  optimized  the  phase  of  the  flutter  mode  and  having  checked  the 
stability  of  the  system  within  the  whole  flight  range  closed  loop  tests  were  per- 
formed. No  valuable  test  data  could  be  gained  by  excitation  through  stick  jerks 
because  the  frequency  content  of  those  excitation  signals  is  too  low.  The  only  meth- 
od to  check  the  closed  loop  condition  is  to  feed  an  additional  sinusoidal  signal 
(frequency  sweep)  into  the  closed  loop  and  analyze  the  response. 

Fig.  13  shows  the  time  history  of  a closed  loop  test  with  frequency  sweep  input. 
Comparing  Fig.  12  with  Fig.  11  one  can  see  that  the  FSS  reduces  the  vane  angle  at 
the  frequency  of  about  ik  Hz.  The  damping  coefficient  was  also  evaluated  from  these 
tests  and  is  presented  later  in  a summary  plot. 
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FIG.  13  TIME  HISTORY  OF  FSS  RESPONSE 
DUE  TO  FREQUENCY  SWEEP  INPUT 


5. Teats  to  Substantiate  a New  Flight  Flutter  Test  Method  for  Stores 

The  aim  of  this  tests  was  to  demonstrate  that  the  modified  FSS  can  be  used  for 
a new  flight  flutter  test  method  as  already  described  ln[l]  . 

Fig.  5 shows  the  principle  of  the  method.  If  the  phase  of  control  signal  is 
shifted  180°  the  FSS  becomes  unstable  and  excites  the  airplane.  This  method  called 
Automatic  Mode  Excitation  (AME)  has  the  following  outstanding  advantages: 

. It  is  automatically  tuning  the  frequency  into  the  store  flutter  mode  and  excites 
it  harmonically  (provided  a suited  sensor  signal  is  choosen). 

. Switching  off  the  AMES  one  can  easyly  evaluate  the  aircraft  damping  from  the 
logarithmic  decrement  of  the  response. 

. Combined  with  the  FSS  it  allows  to  measure  directly  the  damping  trend  of  the  flutter 
mode  above  the  critical  speed. 

In  comparison  with  the  damping  evaluation  from  frequency  sweep  response  which 
cannot  be  done  without  the  help  of  computers  the  signals  produced  by  AMES  can  be 
analyzed  directly  by  the  engineer. 

Fig.  Ik  shows  the  time  history  of  AMES  at  various  airspeeds.  This  signals  were 
produced  during  the  flight  test.  As  can  be  seen  high  signal  to  noise  ratio  allows 
an  evaluation  of  the  damping  from  the  decay  at  once. 
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FIQ.  14  TIME  HISTORY  OF  AMES  v = 450  KIAS  AND  v = 470  KIAS 


The  logarithmic  increment  of  the  increasing  amplitudes  shows  the  excitation  of  the 
stable  system  which  can  be  turned  into  damping  by  shifting  the  phase  l80  . The  dampxng 
evaluated  this  way  has  to  bo  added  to  the  damping  trend  of  the  flutter  mode  to  get  the 
damping  of  the  stabilized  system.  This  is  shown  in  Fig.  15* 
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FIQ.  15  RELATIONSHIP  BETWEEN  FSS  AND  AMES 


5.5  Pemonatration  of  Active  Flutter  Suppreaaion 


During  the  flight  teat  it  vraa  found  that  the  atructural  damping  of  the  A/C  was 
higher  than  expected.  Thla  increased  the  flutter  speed  to  a speed  which  could  not  be 
reached  anymore.  Therefore  an  artificial  flutter  case  was  produced.  It  was  found  that 
the  A/C  could  be  driven  into  flutter  at  any  lower  speed  using  only  one  system  in  the 
AMES  mode. 
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FIG.  16  TIME  HISTORY  OF  THE  UNSTABLE  SYSTEM  AT  v = 300  KIAS 


Fig*  l6  shows  the  time  history  of  an  artiTicial  flutter  case  at  JOO  kts*  Using 
the  other  system  in  FSS»mode  the  instable  system  became  stable. 
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FIG.  17  TIME  HISTORY  OF  STABILIZED  SYSTEM  AT  v = 300  KIAS 


Thla  is  demonstrated  in  Fig.  17.  This  figure  shows  the  unstable  and  stabilized 
system.  As  can  be  seen,  even  additional  excitation  by  stick  Jerks  could  not  destabi- 
lize the  system. 
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6.  FLIGHT  TEST  RESULTS 

All  tests  were  performed  in  the  high  subsonic  regime.  Pronounced  transsonic 

effects  could  not  be  evaluated. 

6.1  Automatic  Mode  Excitation 

The  method  described  was  found  to  be  very  powerful  and  promising.  The  main  advan- 
tages of  the  AMES  are: 

. The  possibility  to  measure  the  flutter  speed  directly  if  it  occurs  within  the 
flight  envelope  of  the  A/C 

. Improved  extrapolation  of  the  d^lmping  trend  is  possible 

. The  test  results  can  be  directly  used  by  the  engineer  without  the  help  of  com- 
puters 

. The  flight  flutter  test  of  stores  can  be  done  with  less  flight  hour?  and  less 
risks  for  the  A/C  (one  store  configuration  in  one  flighti) 


In  Fig.  l8  damping  trends  of  the  aircraft  with  the  flutter  stores  are  presented. 

It  can  be  seen  that  the  damping  evaluated  from  the  frequency  sweep  response  is  slightly 
higher  than  the  values  found  with  AMES. 


FIG.  18  MEASURED  DAMPING  AND  FREQUENCY  TREND 


The  most  important  feature  for  the  system  is  a reliable  automatic  switch  off  of 
the  AMES.  Especially  at  higher  frequencies  and  higher  speeds  it  was  found  that  in  some 
cases  the  pilot  will  not  be  able  to  switch  off  manually  the  system  quick  enough. 


6.2  Discussion  of  Results  with  the  Flutter  Suppression  System 

. Influence  of  the  phase  on  efficiency  of  FSS: 

The  simple  control  system  used  for  this  test  was  optimized  for  one  mode,  the  flutter 
mode . 

As  can  be  seen  from  the  open  loop  tests  non  optimized  phase  reduces  the  system 
efficiency  with  the  cosine  of  the  phase.  Therefore  the  phase  optimization  needs  no 
high  precision.  The  control  system  however  has  to  compensate  the  phase  shift  caused 
by  the  air  forces  on  the  flutter  store  (approx.  20°  from  0 - 500  K1A5). 

. Influence  of  the  gain  on  the  efficiency  of  the  FSS: 

The  damping  force  produced  by  the  FSS  at  constant  airspeeds  is  proportional  to  the 
gain  of  the  system. 
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Testing  of  failure  cases 

In  a failure  case  if  only  one  flutter  tank  works  as  a damper  the  efficiency 
is  reduced.  In  Fig.  19  this  effect  is  shown.  If  the  behaviour  of  the  A/C  is 
symmetrical  this  effect  can  be  regarded  as  a gain  reduction. 


FIG.  19  DAMPING  BEHAVIOUR  OF  FSS  IN  FAILURE  CASE 


If  one  accelerometer  of  the  system  fails  the  control  system  uses  the  vertical 
acceleration  of  one  sensor  as  feed  back  signal.  This  means  theoretically 
a reduction  of  the  gain  of  50  % and  the  control  system  is  a true  ILAF  System  [17]. 

Fig.  20  shows  this  effect  on  the  Nyquist  plots. 


FvMLURE  CASE  NORMAL  SVSTEM 


-0,S  -0,3  -0,1  0,1  0,3  0,5  HO  -0,6  -0,2  0l2  0,6  1,0 

V = 300  KIAS 


FIG.  20  NYQUIST  PLOTS  FOR  FAILURE  OF  ONE  ACCELEROMETER 


These  two  tests  have  shown  that  assuming  the  A/C ' s behaviour  is  sysimetrical 
this  kind  of  FSS  can  be  easy  designed  with  high  redundancies. 


Estimated  Increase  of  the  Flutter  Speed  due  to  FSS 

The  increase  of  the  real  flutter  speed  could  not  be  demonstrated  in  flight 
In  Fig.  21  the  increase  of  the  flutter  speed  is  given  by  extrapolating  the 
measured  damping  trends.  In  this  picture  a medium  gain  was  used.  One  must  keep 
in  mind  that  the  extrapolation  given  in  Fig.  20  does  not  include  the  effects 
which  may  occur  at  high  transsonic  speeds.  Stall  of  the  vanes  and  saturation 
conditions  of  the  actuator  have  to  be  avoided. 


T ‘ T 


STABILIZED 
FLUTTER  MODE 


FLUTTER  MODE 


FIG.  21  EFFICIENCY  OF  FSS 


CONCLUSIONS  AND  RECOMMENDATIONS 

It  was  shown  that  a relative  simple  control  system  with  store  mounted  vanes 
can  be  used  for  store  flutter  suppression.  This  FSS  is  very  effective  with  rela 
tive  small  vanes  which  do  not  change  the  flight  mechanical  characteristics  of 
the  aircraft. 


First  flight  experience  with  an  AMES  was  made  and  the  method  was  found  to 
be  promising. 
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\ AIRPLANE  MATH  MODELING  METHODS 

r FOR  ACTIVE  CONTROL  DESIGN 
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The  Boeing  Company 
3801  S.  Oliver 
Wichita,  Kansas  67210 
U.S.A. 

SUMMARY 

Selected  Analytical  Methods  are  described  which  are  useful  and  practical  in  math  modeling  for  airplane  active  control 
system  design.  A technique  for  writing  state  equations  is  presented  which  is  suitable  for  incorporating  lifting  surface 
aerodynamic  solutions.  An  economical  method  of  computing  unsteady  aerodynamic  influence  matrices  is  presented  for  line 
doublets  and  plate  doublets,  the  latter  usable  at  any  Mach  number.  An  economical  way  to  analyze  three-dimensional 
turbulence  and  a convenient  way  of  using  design  criteria  in  n-dimensions  are  presented  to  aid  in  designing  for  statistical 
performance.  Recommendations  include  the  use  of  a single  airplane  math  model  for  analysis  of  multiple  performance 
parameters  and  the  use  of  control  hardware  math  modeling  during  preliminary  design. 


NOTATION 

A Aerodynamic  influence  matrix  relating  grid  element  force  and  boundary  condition  velocity  (normalwash) 

B Matrix  approximating  function  for  the  product  Hi^  A"^ 

B^  One  of  the  coefficient  matrices  used  in  defining  B 

b An  element  of  the  matrix  B 

C|<  One  of  the  coefficient  matrices  used  in  approximating  F 

e Base  of  natural  logarithms,  2.718... 

F Generalized  aerodynamic  influence  matrix  relating  generalized  forces  and  motions 

f Frequency,  Hertz 

G Power  spectral  density;  also  a generalized  aerodynamic  influence  matrix  relating  generalized  forces  and 

gust  velocities. 

H Frequency  response  function  relating  an  airplane  response  and  a gust  velocity 

H.|  2 3 Modal  matrices,  relating  generalized  motion  and  grid  element  normal  displacements:  (1 ) linear  at  center  of 
pressure,  (2)  angular  and  (3)  linear  at  normalwash  reference  point. 

i Imaginary  number  Also  an  index 

j An  index 

K A factor  used  in  turbulence  design  criteria.  Also  a modified  spherical  Bessel  function 

k An  index 

L Scale  of  turbulence.  Also,  a randomly  varying  structural  load 

Lsteady  average  value  of  a random  structure  load 

M Mach  number 


The  number  of  dimensions  of  a multiple-dimension  variable 
Scaled  Laplace  variable,  s/U„ 


Q 


Vector  (column  matrix)  of  generalized  aerodynamic  forces 


r 
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q 

Generalized  displacement  vector  (column  matrix) 

R 

Correlation  function,  transform  related  to  power  spectral  density 

Rk 

One  of  the  coefficient  matrices  used  in  approximating  G 

s 

Wing  semi-span 

s 

Complex  Laplace  variable 

T 

Rectangular  or  diagonal  matrix  of  gust  penetration  time  delays  (Laplace  transformed)  for  grid  elements  or 
groups  of  grid  elements 

\ 

1 

Uo 

Airspeed 

1 

1 

''9 

Lateral  gust  velocity  vector  (column  matrix) 

! 

'^9 

Vertical  gust  velocity  vector  (column  matrix) 

I 

X 

A superposition  of  responses  (for  example  structural  element  stress).  Also,  the  argument  of  the  function 
K^(x) 

! 

y 

Spanwise  coordinate.  Also,  a typical  response  (for  example  a bending  moment) 

i 

a 

Dummy  argument  of  R.  Also,  a vector  (column  matrix)  of  points  on  a multiple  dimension  sphere 

p 

A denominator  constant  used  in  defining  B and  b 

y 

A denominator  constant  used  in  approximating  G 

r 

Mathematical  gamma  function,  the  generalized  factorial 

/i 

Difference  operator,  requires  subtraction  as  dy  = yg-y.! 

V 

The  order  of  the  function  Ky(x) 

Variance  of  a random  variable 

T 

Time  variable  for  the  correlation  function  R 

U) 

Fourier  variable  or  frequency,  radians/second 

Scaled  turbulence  frequency,  2 rrf/Uo 

Oo 

Turbulence  parameter  = /"(1/2)r(5/6)/L  r(V3)  = 1/1.339L 

{ } 

A column  matrix  (vector) 

• 

1 1 

A rectangular  or  square  matrix 

! 

1.  INTRODUCTION 

j 

A description  of  analytical  methods  found  useful  and  practical  in  math  modeling  for  airplane  active  control  system  design 
and  performance  analysis  is  presented.  The  selection  of  methods,  emphasizing  lifting  surface  aerodynamic  theory  and 
statistical  performance  evaluation,  was  based  not  only  on  their  current  acceptance  by  engineering  specialities  but  also  on  the 
existence  of  format  or  cost  barriers  which  have  prevented  widespread  use  by  control  system  designers.  The  adaptation, 
extensions  and  simplifications  discussed  in  the  following  sections  have  been  used  in  practical  engineering  problems 
requiring  design  of  active  flight  control  systems. 

The  need  for  compatible  airplane  math  models  for  analysis  of  stability  and  control,  flutter,  structural  integrity  and  comfort 
has  been  heightened  by  two  factors:  first,  the  ability  to  improve  performance  in  each  of  these  areas  with  active  control 
systems,  and  second,  the  lack  of  frequency  separation  in  large  airplanes  between  these  areas  of  concern  which  means  all 
performance  measures  must  be  monitored  to  prevent  inadvertent  degradation.  The  experience  basis  for  the  chosen 
methods  is  outlined  by  the  following:  the  development  of  measures  of  structural  performance  (Reference  1 ),  the  design  and 
flight  testing  of  an  active  control  fatigue  damage  reduction  system  (Reference  2),  development  of  the  interface  between 
stress  analysis  and  statistical  loads  analysis  (Reference  3),  and  the  design  and  flight  testing  of  active  control  load  reduction 
and  ride  improvement  systems  (Reference  4)  and  a flutter  suppression  system  (Reference  5). 

i 
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The  following  paragraphs  and  appendices  contain  sufficient  information  for  the  techniques  to  be  used  by  aeroelasticians 
and  active  control  analysts.  The  primary  intent  is,  however,  to  encourage  the  active  controls  designer  to  ask  for  and  use  math 
models  which  are  adequate  for  simultaneous  analysis  of  all  airplane  dynamic  performance  requirements. 

2.  LIFTING  SURFACE  AERODYNAMIC  SOLUTIONS  IN  EQUATIONS  OF  MOTION 

Generally  useful  equations  of  motion  must  be  expressible  in  all  user  formats  including  differential  equations,  difference 
equations,  Laplace  transformed  equations,  Z-transformed  equations  and  Fourier  transformed  equations.  Linearity  of  the 
structural  and  aerodynamic  theory  is  assumed,  although  linearity  may  be  with  respect  to  a nonlinear  “steady-state " solution. 
This  section  is  concerned  with  the  aerodynami"  part  of  the  equations  and  the  unsteady  linear  aerodynamics  typified  by  the 
lifting  surface  theories.  These  are  calculable  by  pressure  series,  doublet  lattice,  Mach  box  and  other  methods.  This  theory  is 
not  directly  suitable  for  equations  of  motion  since  it  cannot  be  represented  by  a finite  number  of  state  variables  in  the  user  s 
transform  formats.  Also,  most  current  solution  methods  are  limited  to  numerical  evaluation  of  the  Fourier  transformed 
airloads  for  selected  frequencies. 

Given  a set  of  Fourier  transformed  airloads  covering  the  frequency  range  of  inputs  and  structural  modes,  the  use  of  an 
approximating  function  for  the  airloads  is  permissible  and  desirable.  The  function  can  fit  the  known  solutions  to  any  desired 
accuracy  and  can  serve  as  an  interpolator  for  other  frequencies.  This  function  must  be  physically  realizable  and  stable  (its 
time  transform,  the  impulsive  admittance  of  airloads  to  boundary  normal  wash,  must  be  nonanticipative,  real,  and  zero  at  time 
= infinity). 


Such  an  approximating  function,  since  it  will  not  accurately  fit  all  frequencies  from  zero  to  infinity,  cannot  be  presumed  to  be 
universally  usable.  For  example,  a Laplace  transform  of  the  function  would  not  be  accurate  in  the  vicinity  of  a zero  or  pole 
introduced  by  the  approximating  function. 

The  usefully  accurate  range  for  Laplace  arguments  is  not  obvious,  although  analytic  continuity  suggests  that  leaving  the 
imaginary  axis  (small  positive  or  negative  damping)  is  comparable  to  interpolation  along  the  imaginary  axis.  The  region  near 
the  imaginary  axis  is  of  greatest  physical  interest. 


Fourier  transform  lifting  surface  airloads,  whether  panel  airloads,  pressure  coefficients  or  othenwise  generalized  airloads, 
are  presumed  equal  to  the  product  of  an  influence  matrix  and  a matrix  of  local  normalwashes.  The  latter  matrix  is  due  to 
gusts,  control  surface  and  airplane  motions,  and  structural  deformations.  The  approximating  function  must  fit  the  elements  of 
the  influence  matrix.  Equations  of  motion  used  in  References  4 and  5 were  based  on  fitting  each  element  of  the  influence 
matrix  with  a function  of  the  form  shown  in  Equation  (1 ).  The  denominator  constants  are  noncritical  and  are  typically  chosen 
as  = <*>rnax^’'*“’o'  numerator  constants  are  selected  to  give  the  least  square  error  for  the  values  ofruat  which  b is 
known. 
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Appendix  I provides  the  computations  required  to  incorporate  these  approximating  functions  into  equations  of  motion  for 
generalized  freedoms.  For  each  generalized  displacement,  the  aerodynamic  plus  structural  terms  require  six  state  variables. 
The  usefulness  of  this  formulation  is  in  regarding  the  imaginary  variable  p as  a complex  scaled  Laplace  variable,  s/Uq.  When 
this  is  done,  eigenvalues  of  the  Laplace  characteristic  equation  appear  as  in  Figure  f . Airplane  rigid  body  and  vibration 
eigenvalues  appear  as  expected.  The  four  clusters  of  eigenvalues  on  and  near  the  negative  real  axis  are  due  to  the 
approximating  functions  and  therefore  not  physically  meaningful.  Q-R  iteration  with  complex  conjugate  shifts  (Reference  6) 
is  adequate  for  this  eigenvalue  problem. 

Equations  of  motion  in  the  form  shown  in  Appendix  I are  satisfactory  for  all  types  of  aeroelastic  analysis  - open  and  closed 
loop  flutter,  divergence,  gust  response,  maneuver  response,  stability  and  control  and  active  control  synthesis  It  is 
sometimes  useful  to  let  the  approximating  function  give  linear  piston  theory  results  when  the  frequency  approaches  infinity; 
however,  this  is  not  required  for  low  frequency  analysis  and  it  does  not  imply  knowledge  of  high  noninhnite  frequency 
response. 

The  differential  equations  deduced  from  the  Laplace  equations  are  used  for  designing  to  criteria  specified  as  time  histones 
and  analysis  of  control  limiting  and  saturation  and  closed-loop  limit  cycles  Other  forms  of  the  equations  are  readily 
determined  if  the  eigenvectors  and  eigenvalues  of  the  Laplace  equations  are  put  in  a real  uncoupled,  first-order  form  as 
illustrated  in  Section  9.67  of  Reference  7.  Eigenvectors  corresponding  to  the  approximating  function  eigenvalues  may  be 
deleted 


Figure  1 . Typical  Eigenvalues  of  the  Laplace  Characteristic  Equation  Using  Approximating  Functions 


3.  USE  OF  UNSTEADY  AERODYNAMIC  INFLUENCE  MATRICES  GENERATED  FOR  PLATE  DOUBLET 
AERODYNAMIC  MODELS 


The  use  of  active  control  surfaces  on  airplanes,  especially  several  small  surfaces,  leads  to  complicated  and  expensive 
airplane  math  models.  The  following  discussion  concerns  the  use  of  lifting  surface  aerodynamic  theory  through  finite  element 
aerodynamic  doublet  modeling,  the  conversion  of  steady-state  solutions  to  unsteady  solutions,  and  the  successful  use  of 
unsteady  elements  analogous  to  the  Woodward  steady-state  trapezoidal  element  (Reference  8). 

Lifting  surface  theory  is  derived  from  the  acoustic  wave  equations.  This  theory  is  valid  at  any  Mach  number  if  the  boundary 
conditions  imposed  do  not  require  significant  pressure  perturbations  compared  with  static  pressure,  or  equivalently,  do  not 
require  significant  velocity  perturbations  compared  with  the  speed  of  sound.  When  the  solutions  are  presented  in  the  Fourier 
transform  domain,  steady-state  lifting  surface  theory  appears  as  a special  case  when  the  frequency  is  zero.  Another  special 
case  occurs  when  the  boundary  conditions  on  a region  of  the  wing  are  constant  for  a radius  large  compared  with  the 
wavelength  of  the  motion;  then  the  solution  of  lifting  surface  theory  for  that  region  will  be  the  same  as  that  of  linear  piston 
theory. 


Finite  element  aerodynamic  solutions  require  dividing  the  airplane  into  a fine  enough  grid  so  that  the  user  will  need  to 
specify  only  one  boundary  condition  and  compute  only  one  aerodynamic  force  for  each  grid  element.  When  each  element 
influences  every  other  element,  a square  influence  matrix  must  be  computed.  Two  well-known  methods  of  computing  the 
influence  matrix  are  doublet  lattice  and  Mach  box,  see  (for  example)  Reference  9.  The  former  method  is  suitable  for  subsonic 
flow  only,  the  latter  only  for  supersonic  flow.  The  latter  requires  that  grids  must  be  Mach  line  oriented  and  changed  for  each 
new  analysis  Mach  number. 

A Kernel  function  (see  Reference  10)  giving  velocity  perturbations  due  to  an  aerodynamic  force  applied  at  a point  can  be 
readily  computed  at  any  Mach  number.  Hedman,  in  Reference  11,  analytically  integrated  the  steady-state  Kernel  along  a 
spanwise  line  and  constructed  a subsonic  influence  matrix  using  these  line  elements  (vortex  lattice).  Woodward,  in 
Reference  8,  integrated  the  Kernel  over  a trapezoidal  "plate " element  and  showed  the  possibility  of  using  a single  element  at 
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all  Mach  numbers.  The  authors  of  Reference  10  developed  a practical  unsteady  form  of  the  line  element  method  (doublet 
lattice)  and  recommended  an  extension  of  the  Woodward  method  to  unsteady  flow. 

The  author  uses  a simple  scheme  for  converting  vortex  lattice  influence  matrices  to  doublet  lattice.  The  ratio  of  the  normal 
component  of  the  unsteady  Kernel  to  the  steady  Kernel  is  a well-behaved  function  that  shows  gradual  changes  in  phase  and 
amplitude  over  regions  the  size  of  doublet  lattice  elements.  To  calculate  the  unsteady  (complex)  normalwash,  the  author 
merely  multiplies  the  integrated  steady  (real)  normalwash  by  this  ratio  obtained  at  a point  in  the  element  near  the  normal- 
wash  reference  point.  For  practical  size  grid  elements,  this  method  is  as  accurate  as  numerically  integrating  the  unsteady 
Kernel  along  a line  and  it  is  not  erratic  for  high  aspect  ratio  elements. 

A similar  scheme  was  successful  for  converting  steady  plate  element  influence  matrices  to  unsteady,  which  permits  the 
use  of  “plate  doublet  aerodynamic  models.”  The  point  in  the  receiving  element  where  the  unsteady  to  steady  Kernel  ratio  is 
computed  is  chosen  as  deep  as  possible  within  the  Mach  cone  of  the  sending  (pressure)  element.  It  is  likely  that  other  popular 
steady-state  solutions,  including  those  for  finite  thickness  wings,  could  be  converted  to  unsteady  solutions  by  the  same 
method. 

Three  interpretations  or  restrictions  of  the  Woodward  solution  were  required.  First,  an  influence  matrix  for  a particular  grid 
pattern  is  an  entity  in  itself,  independent  of  the  aerodynamic  element  used  to  derive  it,  and  consequently  does  not  imply  any 
particular  location  tor  the  center  of  pressure  within  the  element.  When  generalizing  forces,  the  vortex  lattice  convention  was 
used  subsonically  which  places  the  center  of  pressure  at  the  element  quarter-chord,  mid-span.  When  the  element  leading 
edge  was  supersonic,  the  center  of  pressure  was  placed  at  mid-chord,  mid-span  and  the  location  was  interpolated  tor 
in-between  Mach  numbers. 

Second,  the  location  of  the  normalwash  reference  point  within  the  element  was  always  0.8808  chord,  mid  span,  where 
0.8808  = 1/2  + 1/2  (e2  -1)/(e2  f 1).  This  value  was  chosen  so  a single  isolated  rectangular  element  of  aspect  ratio 
approaching  infinity  would  have  a subsonic  lift  coefficient  of  2 The  supersonic  lift  coefficient  of  such  an  element  will 

be  a//m2  . 1 regardless  of  the  normalwash  location.  A similar  analysis  of  the  line-doublet  element  will  show  the  subsonic 
normalwash  location,  if  at  mid-span,  must  be  one-half  chord  aft  of  the  pressure  line,  but  that  no  (real)  supersonic  location  is 
possible. 

Third,  although  the  Woodward  scheme  may  be  readily  programmed  for  trapezoidal  elements,  only  parallelogram  elements 
are  consistently  satisfactory.  In  calculating  the  steady  influence  matrix,  the  trapezoidal  grid  element  is  replaced  by  a 
parallelogram  pressure  element  which  has  the  same  chord  at  mid-span  and  the  same  sweep  at  mid-chord. 

Results  of  sample  computations  are  presented  in  Figure  2 for  the  AGARD  supersonic  planform  These  are  compared  with 
the  Mach  box  results  published  in  Reference  12.  This  method  gives  reasonably  close  answers  for  coarse  grids  and 
converges  smoothly  and  quickly  for  finer  grids.  Subsonic  solutions  are  equivalent  in  accuracy  to  doublet  lattice.  Solutions  at 
any  Mach  number  cost  about  the  same  as  doublet  lattice.  The  grid  chosen  can  be  used  at  all  Mach  numbers  without  change 
and  solutions  are  available  at  all  transonic  Mach  numbers.  Planar  doublet  aerodynamic  models  are  recommended  for  act.ve 
control  design  and  performance  evaluation. 

4.  THREE-DIMENSIONAL  TURBULENCE 

The  basic  concepts  of  dynamic  analysis  for  multiple  random  inputs  are  well  understood  and,  given  a suitable  mathematical 
model  of  the  input,  the  designer  of  active  control  systems  does  not  have  difficulty  with  design  synthesis.  In  the  important  case 
of  atmospheric  turbulence,  adequate  math  models  exist  but  are  cumbersome  and  not  widely  applied.  A convenient  and 
inexpensive  method  is  presented  for  modeling  three-dimensional  turbulence,  developed  from  the  work  of  Sawdy  in 
Reference  13. 

The  most  noticeable  improvement  in  analyses  that  use  three-dimensional  turbulence  is  in  calculation  of  the  variance  of 
generalized  mode  velocities  and  accelerations,  and  therefore  in  the  calculation  of  zero-crossing  and  maxima  rates.  The 
reason  for  change  is  that  parallel  components  of  turbulence  in  a plane  perpendicular  to  the  flight  path  are  only  partially 
coherent  in  three-dimensional  turbulence.  For  example,  the  vertical  component  of  turbulence  across  the  span  of  a straight 
wing  airplane  will  at  any  particular  instant  have  some  random  variation  instead  of  being  identical  as  in  one-dimensional 
turbulence  analysis.  The  coherence  between  a vertical  gust  at  the  left  wing  tip  and  a vertical  gust  at  the  right  wing  tip 
decreases  rapidly  as  frequency  increases.  If  this  decrease  is  neglected,  most  airplane  math  models  will  predict  infinite 
zero-crossing  rates  for  generalized  freedoms. 

The  spectral  density  of  a response  to  the  vertical  component  of  turbulence  is  given  (see  Reference  14,  eq.  3.178)  in 
Equation  (2)  where  Gyy  yj,  = cross  spectral  density  of  the  ith  and  jth  vertical  gust  as  shown  in  Figure  3,  and  where  H includes 

the  time  delay  caused  by  gradual  penetration  of  the  gust. 


G (f):  Hyyyjff)  G^W  ( ^ ^ ^ W 

' i I I I 


(2) 
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Figure  2.  Comparison  of  Mach  Box  and  Plate  Doublet  Force  Coefficients  for  the  AGARD  Supersonic  Wing 

The  summation  is  used  as  an  approximation  of  the  spanwise  integration  shown  in  Reference  13.  The  summation  extends 
across  the  full  span  of  the  airplane  and  vertical  gusts  excite  both  symmetric  and  antisymmetric  modes.  Since  the  gust 
penetration  time  delays  are  included  in  the  frequency  response  functions,  the  gust  cross  spectral  densities  are  real.  They  are 
given  for  a unit  standard  derivation  by  the  correlated  spectra  equation  (A31 ) of  Reference  1 3.  Sawdy's  equation  is  repeated 
here  as  Equation  (3). 


[l  + 8/3  (fi,  (^o 
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BOUNDARY  CONDITIONS  DETERMINED  BY 

Wj  AND  PENETRATION  TIME  DELAY 
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Figure  3.  Spanwise  Variation  of  Vertical  Gust 

In  Equation  (3),  is  the  spanwise  separation  distance.  Using  the  limit  shown  in  Equation  (4),  the  special  case  of  zero 
spanwise  separation  reduces  to  the  familiar  expression  in  Equation  (5). 

lim  K„(x):1/2r(i»)(2/x)''  (4) 


2L  [1^  8/3(fli/flo)^] 


The  gust  cross  spectral  density  may  be  identified  alternafely  as  in  Equation  (6), 


+ (-^)  )cos(2,r 


fr)dT 


where  Rw,vvj  — ^ 'S  the  von  Karman  autocorrelation  function  for  vertical  gusts  and  is  given  in  Equation  (7). 

The  T symmetry  is  apparent  in  this  form  and  is  the  reason  the  cross  spectral  density  is  real. 


Equation  (6)  is  not  recommended  for  evaluafing  Gy^yv  because  of  computing  cost.  However,  the  equation  is  valuable  for 
illusfrating  an  alternate  derivation  and  the  significance  of  the  forms  from  Reference  13.  A rapid  method  for  calculating  the 
modified  spherical  Bessel  functions  in  Equation  (3)  is  presented  in  Appendix  II. 

In  Figure  4,  which  is  based  on  Equation  (3),  the  cross  spectral  density  of  two  gusts  at  points  with  spanwise  separation  dy 
has  been  normalized  by  the  gust  autospectral  density  ( Ay  = 0).  Flag  1 is  plotted  at  dy/L  = .2  and  wL/Uo=4  and  is 
representative  of  the  turbulence  which  drives  large  airplane  low  frequency  modes  at  low  altitudes  (L  = approximately  150 
m.).  If  fhe  scale  of  turbulence  "L"  were  doubled,  the  point  would  be  replotted  at  Flag  2,  but  it  does  not  move  appreciably 
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upward  toward  the  1 -D  gust  value.  If  the  scale  of  turbulence  were  four  times  the  original,  the  point  would  be  at  Flag  3.  If  the 
scale  were  ten  times  the  original,  the  point  would  be  at  Flag  4,  where  the  spectral  ratio  is  still  not  significantly  higher  than  for 
the  original  scale.  A similar  trend  can  be  shown  for  higher  frequency  responses.  In  an  airplane  analysis,  only  responses  near 
zero  frequency  converge  to  one  dimensional  gust  values  for  reasonable  scale  lengths. 


See  Section  4 for  a discussion  of  the  flags. 


/4y/L  = 0 


4y/L  = .01 
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Figure  4.  Effect  of  Spanwise  Separation  on  Vertical  Gust  Cross-Spectral  Density 
5.  MULTIDIMENSIONAL  STATISTICAL  PERFORMANCE  CRITERIA 


The  use  of  multidimensional  structural  performance  criteria  is  further  complicated  when  those  criteria  are  statistical.  The 
structures  analyst  must  be  concerned  with  "Where  does  the  active  control  system  designer  stop?"  For  example,  does  he 
need  to  know  the  damage-to-stress  relationship,  or  the  stress-to-load  relationship,  or  even  the  load-to-state  variable 
relationship?  The  answer  depends  on  the  criteria.  The  following  is  offered  as  a tool  since  such  criteria  sometimes  are 
selected  at  the  convenience  of  the  analyst. 

If  X,  consisting  of  Xsteady  + AX,  is  any  superposition  of  "n"  responses  (each  consisting  of  Lsteady  + AL),  the  design 
criteria  for  X is  that  (AX)2  be  less  than  a factor  times  its  variance  in  unit  variance  turbulence,  as  in  Equation  (8), 

(AX)2<K<.*2 

Then  all  points  on  the  design  surface  in  space  L are  mappable  using  Equation  (9a).  In  Equation  (9a),  o ^ll  'S  a square 
symmetric  matrix  of  covariances  of  the  responses,  L,  and  any  of  its  symmetric  square  roots  is  acceptable.  The  matrix  (o^ 
contains  the  coordinates  of  a point  on  the  surface  of  an  n-dimensional  sphere  obtained  by  solving  Equation  (9b).  Thus,  the 
control  system  designer  can  work  in  the  space  L to  define  its  covariance  matrix,  knowing  the  X - to  - L relationship  can  be 
designed  later  in  such  a manner  that  X can  meet  its  criteria. 


I'-Irll.sTEADyl  + K^o’LL] 


1 


Equation  (9)  was  developed  in  Reference  3 to  provide  a design  loads/design  stress  interface.  Other  procedures  for  more 
complicated  criteria  are  outlined  in  the  same  reference  and  could  be  adapted  to  more  general  active  control  design  use. 
Figure  5 illustrates  a two-dimensional  case,  shear  and  bending  moment,  mapped  for  a mission-weighted  exceedance 
criteria.  This  figure  was  developed  without  knowledge  of  the  stress-to-load  relationship.  Using  this  figure,  engineers  can 
design  structure  which  has  acceptable  stresses  for  load  combinations  on  and  within  the  envelope  shown. 


Figure  5.  Example  Presentation  of  Two-Dimensional  Design  Envelope  for  Statistical  Criteria 


6.  RECOMMENDATIONS 

The  implementation  of  active  controls  in  flight  hardware  and  the  flight  attainment  of  predicted  airplane  structural  benefits 
are  achievable  with  today's  analytical  methods  and  technology.  Successful  design  is  based  on  first  predicting  the 
characteristics  of  feasible  hardware  and  structure  and  the  flight  environment  and  then  synthesizing  control  systems. 
Minimum  requirements  for  the  math  models  are  summarized  below. 

Control  hardware  math  models  should  include  a nonlinear  simulation  of  the  preliminary  design  actuator.  If  realistic 
tolerances  for  backlash  and  hysteresis  prevent  attainment  of  performance  goals,  smaller  control  surfaces  with  larger 
amplitude  motions  should  be  considered.  Control  system  math  models  should  permit  the  designer  to  trade  electrical  versus 
mechanical  functions.  For  example,  the  low-pass  characteristics  of  an  actuator  may  be  less  desirable  than  those  of  a 
sprecially  designed  electronic  circuit.  Low-pass  filtering  must  be  completed  within  the  range  of  validity  of  truncated  mode  math 
models. 

Airplane  math  models  should  be  available  which  are  adequate  for  analysis  of  flutter,  dynamic  loads,  and  stability  and 
control.  The  active  control  designer  must  synthesize  his  system  using  these  math  models  at  the  preliminary  design  stage. 
The  use  of  '-.implified  models  and  optimization  procedures  is  not  precluded,  but  a system  should  not  be  considered  a 
preliminary  design  until  its  performance  is  adequate  on  complete  math  models.  Providing  such  math  models  to  the  active 
control  designer  in  a form  he  can  use  has  been  the  subject  of  previous  sections. 

The  author  recommends  the  adoption  of  state-variable  equations  by  all  analysts  involved  with  airplane  dynamics  as  a 
standard  interface  and  means  of  information  exchange.  Improvements  in  math  modeling  theory  then  become  as  transferable 
as  data  updates  When  a common  format  is  established,  the  economy  and  convenience  of  more  complete  math  models  will 
be  self-evident. 
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APPENDIX  I - GENERALIZED  EQUATIONS  OF  MOTION  WITH  APPROXIMATING  FUNCTIONS 

The  generalized  forces  are  defined  in  Equation  (l-la,  b,  c)  for  the  interpolating  variable  "p".  The  influence  matrix  |A(p)|'^ 


{Q(p)|=  [F(p>]  |q(p){  f[G(p)]  1 ® 1 ' ** ' *^  ^ * U_ 

(l-1a) 

[F(p)]  : [A(P>  J’([H2]  + p[H3]) 

(Mb) 

[G(p)]=  [hi]^[a(p)]'^[t(p»] 

(Me) 

.Jit 


and  the  gust  penetration  matrix  |T  (p)|  are  known  only  for  discrete  values  of  p.  Define  an  approximating  function  as  in 
Equation  (1-2), 


[H,f  [A(pl]'»[B„] 


which  is  the  matrix  equivalent  of  Equation  ( 1 ),  where  the  | B|^  | are  chosen  for  a least  square  error  fit  of  the  known  values,  and 
where  the  are  preselected  as  in  Section  2.  Substituting  Equation  (1-2)  in  Equation  (Mb)  gives  Equation  (l-3a,  b,  c,  d). 


[F(p>]~  p[c.i]+[Co]  J 

[C-i]  -([Bo]  +[bi]  + • • ■ + [B4])  [H3] 


(l-3b) 


[Co]  ••  [Bq]  [”2] 


(|-3d) 


The  gust  coefficient.  Equation  (Me)  may  be  fitted  with  a separate  approximating  function  as  in  Equation  (l-4a,  b),  in  which 
|R,|  2 I provide  a least  square  error  fit.  The  author  lets  the  penetration  matrix 

[G(pI]~[Rq]+  [R2]  (l-4a) 


(l-4b) 


lT(p)|  represent  local  groupings  of  aerodynamic  elements  with  the  time  delay  to  the  leading  edge  of  the  grouping  being 
supplied  at  solution  time. 

APPENDIX  II  - NUMERICAL  EVALUATION  OF  BESSEL  FUNCTIONS  IN  CROSS-SPECTRAL  MATRIX  FOR 
VON-KARMAN  TURBULENCE 

The  use  of  Equation  3 in  Section  4 requires  evaluation  of  a large  number  of  modified  spherical  Bessel  functions  of  fractional 
order.  The  following  procedure  is  an  acceptable  compromise  of  speed  and  accuracy  for  the  range  of  orders  and  arguments 
significant  in  turbulence  analysis. 

One  of  the  integral  representations  of  the  function  shown  in  Equation  (11-1 ) is  chosen  according  to  the  values  of  the  order 
and  argument.  These  are  better  suited  for  numerical  integration  than  the  integral  in  Equation  (6),  Section  4.  These  formulas 
were  selected  from  those  published  in  Reference  15,  Paragraph  9.6.  The  integration  is  accomplished  numerically  using 
15-point  Laguerre  weighting  factors  and  abcissas  (see  Reference  15,  Table  25,9). 


Ci,(x):  # e cos  h(r/t)dt,»/<  1/2  or  x>1/2 

o 

r(t'+1/2)/  2 Y /cos(xt)dt 

\ * )/{1  + t2)»'+  1 


,r/>1/2  AND  X < 1/2 
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ABSTRACT 

Military  Specifications  (MILSPECS)  are  often  employed  by  the  USAF  procuring  authority  as  guidelines 
for  design,  development,  acceptance  testing  and  mission  application  of  military  aircraft.  The  MILSPECS 
must  usually  be  satisfied  by  fornuilat ions  of  the  aerodynamic  and  dynamic  analyses  that  are  consistent,  or 
eauivalent  if  not  identical.  When  control  configured  vehicle  considerations  are  involved,  however,  incon- 
sistencies resulting  from  analysis  expediency  or  previous  engineering  convention  may  occur.  In  this 
paper  YF-16,  C-5A,  B-52E  and  large  transport  aircraft  design  studies  and  flight  tests  provide  data  for  a 
discussion  and  numerical  illustration  of  these  inconsistencies.  It  is  concluded  that  they  may  be  mini- 
mized or  avoided  altogether  if  flight  control  specialists  become  more  familiar  with  restrictions  of 
present~day  unsteady  areodynamic  theory,  and  structural  specialists  increase  their  knowledge  of  modern 
dynamics  and  control  theory. 

1.0  INTRODUCTION 

Several  recent  flight  test  programs  sponsored  by  the  Flight  Control  Division  of  the  Air  Force  Flight 
Dynamics  Laboratory  (AFFDL)  have  demonstrated  the  possibility  of  improving  aircraft  performance  thru  less 
Inhibited  use  of  the  flight  control  system.  Within  the  USAF  and  much  of  the  aerospace  industry  this  tech- 
nological development  has  been  given  the  acronym  CCV,  the  Control  Configured  Vehicle.  One  advantage  of 
the  CCV  concept  is  that  several  of  the  structural  design  problems  may  be  solved  with  high-gain,  redundant, 
automatic  control  systems,  rather  than  redefining  the  structure  or  redistributing  mass.  In  fact  CCV  im- 
plies that  the  shape  of  the  aircraft  may  in  some  cases  be  determined  by  the  control  systems  required  to 
meet  mission  goals. 

A number  of  CCV  concepts  have  been  demonstrated  in  flight,  both  In  part  and  as  a whole  for  selected 
aircraft  configurations;  a substantial  research  effort  to  develop  tools  for  CCV  analysis  is  in  progress, 
it  will  therefore  be  assumed  in  this  paper  that  CCV  is,  or  soon  will  be,  a practicable  aircraft  design 
option  and  will  indicate  flight  control  systems  that  purposefully  alter: 

•Static  stability  and  other  handling  qualities  (8785B) 

• Maneuver  loads  (8860A,  8861A) 

• Fatigue  life  (8866a,  9^900) 

• Ride  quality  (9^900) 

• Flutter  mode  damping  (8870A,  9^90D) 

The  numbers  in  parentheses  refer  to  selected  Military  Specifications  (Refs,  1-6)  that  include  considera- 
tions of  interaction  of  the  control  system  and  structure. 

An  examination  of  these  MILSPECS  as  design  guidelines,  with  CCV  in  mind,  shows  that  each  recognizes 
the  cross-dlscipl inary  nature  of  these  control  systems.  A partial  listing  of  relevant  paragraphs  (Table 
1),  shows  the  degree  of  consideration  to  vary  from  specific  phase-  and  gain-margin  requirements,  in 
MILSPECS  8870A  and  9^900,  to  more  general  statements  the.**-  leave  interpretation  of  a particular  flight 
vehicle  to  the  discretion  of  the  contractor  and  the  procuring  authority.  Clearly,  it  is  the  intent  of 
the  MILSPECS  to  guarantee  to  the  procuring  authority  that  interaction  of  the  flight  control  system  with 
the  structure  of  conventional  and  CCV  aircraft  is  properly  considered  in  design,  development,  acceptance 
testing,  and  mission  application. 

It  is  not  the  objective  of  this  paper  to  delve  into  the  origin  and  background  of  the  MILSPECS,  nor  to 
address  a specific  redundancy  or  omission  in  current  specifications.  Its  purpose  is  to  indicate  several 
problems  of  analysis  consistency  that  arise  as  the  individual  flight  control  and  structural  engineering 
disciplines  complete  their  tasks  in  fulfillment  of  design  guidelines  such  as  the  MILSPEC  requirements. 
These  inconsistencies  In  aerodynamic  and  dynamic  formulations  occur  in  practice  because  each  engineering 
discipline  introduces  different  mathematical  approximations  and  because  analyses  require  certain 
expendiencles. 

Inconsistency  in  aerodynamic  formulation  appears  to  occur  primarily  because  each  of  the  disciplines 
has  a different  responsibility  in  aircraft  design.  The  problem  is  particularly  acute  as  it  pertains  to 
interaction  between  the  control  system  and  the  structure.  The  flight  control  specialist,  who  is  charged 
with  responsibility  for  control  system  synthesis,  employs  the  unique  insight  afforded  by  aerodynamic 
approxlmat ions  in  the  frequency,  Laplace,  and  time  domains.  Recent  design  studies  of  CCV  configurations 
indicate  that  in  the  future,  time  domain  aerodvnamic  approximation  will  assume  added  importance  as  very 
complex  multi-loop  control  systems  are  synthesized  via  modern  control  theory.  The  structural  specialist, 
who  is  accountable  primarily  for  control  system  analysis,  and  not  synthesis,  employs  most  commonly  aero- 
dynamic approximations  in  the  frequency  domain.  There,  an  application  of  classical  control  theory  Is 
often  adequate  to  demonstrate  compliance  with  design  guidelines.  It  is  clear  that  analysis  inconsis- 
tencies can  occur  if  each  discipline  individually  utilizes  its  own  approximation  of  the  aircraft  aero- 
dynamics in  these  various  mathematical  domains. 


I 
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Inconsistency  in  the  dynnmicnl  description  of  the  aircraft,  ylven  a ccMiwion  aerodynamic  formula!  ion, 
usually  arises  frcjm  (I)  procedures  used  to  reduce  the  number  of  elastic  decrees  of  freedom  (OOF);  (2) 
deletion  of  the  Initial  condition  of  state  In  the  perturbat Ion-dynamics  ions ; and  (3)  exclusion  of 

the  forward-speed  dcqree  of  freedom.  Studies  of  s i ng le- i nput/mul I i pi e-outpul  CCV  control  systems  (c.g., 
flutter  mode  control  alone)  have  not  adequately  pictured  the  overall  system  dynamics  problems,  particu- 
larlly  those  involving  handling  qualities  and  several  CCV  concepts.  Also,  previous  studies  of  multi- 
Input/mul t ipic-output  CCV  control  systems,  Intended  for  conventional  aircraft,  have  not  highlighted  the 
aforementioned  dynamics  problcrus  because  of  the  inherent  natural  aerodynamic  stability  of  the  aircraft 
and  the  conventional  dynamic  approximat ions  in  the  case  of  a relatively  large  frequency  separation 
between  the  mean-axis-system  modes  (rigid  body  modes)  and  the  structural  dynamic  modes. 

The  approach  followed  herein  is  to  review  briefly  the  CCV  analysis  and  synthesis  task,  in  section 
2.0,  prior  to  discussing  In  sections  3 0 and  ij.O  the  i neons  i s tency  that  arises  in  aerodynamic  and  the 
dynamic  formulations.  The  reviev;  is  considered  in  order  to  define  a general,  modern  control,  state 
space  terminology  which  v/ill  condense  the  description  of  CCV  aircraft  into  a compact  set  of  matrix 
differential  equations  that  may  then  be  subsequently  referred  to  in  the  discussion.  Classical  and 
modern  control  theories  are  assumed  to  be  well  understood  so  that  emphasis  may  be  placed  on  detailed, 
numerical  aircraft  design  or  flight  test  results  for  the  B-52E,  C-5A,  YF-16,  and  a large  transport- 
category  aircraft.  These  numerical  results  are  presented  graphically  in  terms  of  Bode,  Nyquist,  root 
locus  and  time  history  data  plots.  Section  5-0  contains  concluding  remarks.  Throughout  this  discussion, 
one  of  the  main  goals  of  the  Conference  Committee  on  the  Structural  Aspects  of  Active  Controls  should  be 
recal led: 

"...  to  consider  what  the  make  up  or  ingredients  should  be  for  good  and  appropriate  criteria  and 
sped  f icat  ions." 

2.0  BACKGROUNO:  CONTROL-SYSTEM  ANALYSIS  AND  SYNTHESIS 

The  synthesis  of  a CCV  or  of  a conventional  control  system  may  be  accomplished  by  classical  control 
theory,  modern  control  theory,  or  more  likely  a blending  of  the  two,  since  performance  criteria  are 
specified  In  both  frequency  and  time  domains.  In  synthesis,  the  problem  is  the  development  of  a control 
system  guaranteeing  that  relative  stability,  steady-state  accuracy,  transient  response  and  frequency 
response  of  the  control  system  are  adequate  to  satisfy  the  performance  criteria.  In  analysis  the  control 
system  already  exists  and  analytical  or  experimental  methods  are  employed  to  measure  the  performance  of 
the  controlled  dynamical  system.  A cursory  review  of  design  guidelines  such  as  the  MILSPEC  paragraphs  in 
Table  1 indicates  that,  for  the  most  part,  the  structural  specialist  presently  acts  as  an  analyst  who  Is 
supplied  a control  system  as  a portion  of  his  dynamical  system.  If  the  structural  specialist  finds  a 
structural  problem  induced  by  the  control  system,  the  only  practical  options  are  (1)  to  send  the  control 
system  back  to  the  flight  control  staff  for  redesign  (synthesis);  (2)  to  alter  structural  layout  or  mass 
distribution  to  compensate  for  the  controller;  or  (3)  to  alter  the  control  system  himself.  If  the  two 
disciplines  develop  identical  or  equivalent  aerodynamic  and  dynamic  formulations  for  the  unaugmented 
aircraft,  the  first  choice  is  the  logical  course  of  action.  If  it  is  determined  that  the  formulations 
differ,  the  second  two  options  might  seem  viable;  they  are  certainly  not  advisable  because  the  control 
system  design  Is  likely  to  diverge. 

2.1  Definition  of  Problem 

Specifically,  in  synthesis  the  goal  is  to  establish  a relationship  between  measurements,  and 
controls  6_,  of  the  form  £(t)  = K*  ^(t)  where  K*  is  a matrix  of  control  system  gains.  The  measurements  and 
responses  are  related  algebraically  to  a first  order,  ordinary  differential  equation  defining  the  states, 
X,  that  Is,  in  the  modern  control,  state  space  form: 


State  equation; 


x(t)  • A x(t)  + B ^(t)  + n| (t) 


Measurement  equation:  ^(t)  • M )<(t)  + ^ ^(t)  (lb) 

Response  equation:  “ C )<(t)  + D ^(t)  (Ic) 

Control  equation:  6(t)  » K*^(t)  (Id) 

Here  the  matrix  equations  may  be  nonlinear,  with  time-varying  coefficients  in  the  matrices  A,  B,  C,  D,  M, 
and  K*.  Typically,  the  states,  )^,  consist  of  reference-axis-system  motions,  kinematics,  elastic 
deformations,  and  other  variable  associated  with  specifics  of  the  control  system  and  its  modelling, 
The  controls  consist  of  aerodynamic  control  surface  motions,  engine  thrust  and  other  means  of  pur- 
posefully inducing  forces  and  moments  on  the  aircraft.  The  measurements,  consist  of  accelerations, 
rotation  rates,  air  data,  inertial  planform  data  and  other  sensed  quantities.  The  responses,  usually 
consist  of  structural  moments,  shears  and  stresses,  model  following  errors,  control  rates,  tracking 
errors  and  other  outputs  of  the  system.  The  term  n_j(t)  is  process  noise  associated  v/ith  atmospheric 
gusts  and  turbulence  or  pilot  commands,  and  measurement  error  induced  by  the  electronics  or 

mechanics  of  the  sensors. 


It  is  assumed  that  equation  (I)  defines  a dynamical,  controlled  system  that  is  analog  (or  continuous). 
A state-space  form  of  equation  (I)  also  exists  for  a discrete-data  digital  system;  it  will  not  be  discussed 
as  most  modern  control  theory  applies  equally  well  to  digital  control  systems,  and  because  many  digital 
control  systems  result  from  "digitization"  of  previously  designed  analog  systems.  For  all  future 
discussions  it  will  be  assumed  that  equation  (I)  has  constant  coefficients  and  has  been  linearized  in  the 
Independent  variables  of  state,  )«,  and  controls,  It  should  be  recognized  that  the  linearized,  constant" 
coefficient  problem,  although  the  most  common  in  the  literature,  is  only  a portion  of  the  total  CCV 
analysis  and  synthesis  task.  In  actual  design  a substantial  amount  of  engineering  time  must  be  spent  in 
accounting  for  non-linearities  in  aerodynamics,  structure,  and  dynamics  and  control;  these  non-lincar- 
Itles  are  not  discussed  herein.  It  Is  also  assumed  tliat  the  structural  and  the  flight  control  specialists 
have  agreed  upon  the  significant  terms  constituting  the  control  equations,  but  that  each  may  have  a 
different  approximation  of  the.  state,  measurement  and  response  equations  owing  to  differences  in  aero- 
dynamic or  dynamical  formulation. 
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2.2  Classical  Control  Theory 

In  classical  control  theory  determination  of  the  control  equation  and  K*  consists  of  a trial-and- 
error  application  of  analysis  procedures  until  the  performance  criteria  are  satisfied.  The  procedures 
are  largely  a graphical  portrayal  of  Bode,  Nyquist,  and  root  locus  stability  criteria  for  open-loop  and 
closed-loop  systems.  The  measurement  noise  is  usually  ignored  in  all  but  specialized  studies. 

• Bode  and  Nyquist.  The  Bode  graphical  method  determines  relative  stability  information.  The 
Nyquist  method  is  part iculari ly  useful  in  that  it  treats  stability  of  both  minimum  and  non-minimum 
phase  system.  In  both  graphical  procedures  the  open-loop  transfer  function  is  considered,  and  the 
phase  and  gain  margins  (the  design  criteria  found  in  MiLSPEC  8870)  are  determined  from  a plot  of 
the  amplitude  and  the  phase  angle  of  the  transfer  function,  with  frequency,  u>,  varying  from  zero 
to  infinity.  One  advantage  of  these  methods  is  that  they  may  be  used  with  frequency-dependent, 
unsteady  aerodynamic  forces.  A disadvantage  is  that  a mul t i p le- i nput/mu 1 1 i pi  e-output  system 

(typ  ical  of  Cev)  must  be  analyzed  one  loop  at  a time  for  its  stability  characteristics;  thus,  a 
change  in  a compensating  control  element  requires  reanalysis  of  ail  the  loops.  Once  the  system  is 
designed  to  be  stable,  Nichol's  chart  may  be  used  to  convert  the  open-loop  magnitude  and  phase 
relationship  to  closed-loop  frequency-response  characteristics. 

• Root  Locus.  This  method  is  also  graphical  and  requires  that  the  analyst  define  his  system 
[equation  (l)J  in  the  Laplace  domain.  Both  the  open-  and  closed-loop  characteristics  of  a 

mul t i pie- input/mul t i pie-output  system  may  be  determined.  Much  of  the  information  contained  in  the 
Bode,  Nyquist,  and  Nichol's  chart  may  be  calculated  directly  from  the  root  locus  definition  of  the 
dynamical  system.  The  advantages  of  the  method  are  substantial  in  that  closed- loop  stability  can 
be  directly  determined  by  working  with  the  open-loop  transfer  functions  of  the  overail  system. 
Transient  responses  and  steady  state  accuracy  may  also  be  easily  determined  by  an  inverse  Laplace 
transform  to  the  time  domain.  The  effect  on  the  stability  of  the  dynamical  system  of  a gain  change 
in  one  loop  of  the  control  system  may  be  determined  directly.  In  addition,  a root  locus  on  the 
dynamic  pressure,  on  sensor  location  or  type  and  on  the  physical  parameters  of  the  system  (such  as 
those  shown  and  discussed  subsequently  in  Figure  II)  make  the  method  ideally  suited  to  synthesis 
by  reanalysis.  The  method  requires  that  unsteady  aerodynamic  information  in  equation  (I)  be 
expressed  as  constant  coefficient,  ordinary  differential  equations  in  the  Laplace  domain.  This 
disadvantage  is  discussed  subsequently  in  section  3.0. 

2.3  Modern  Control  Theory 

In  modern  control  theory  synthesis  of  a CCV  control  system  proceeds  more  directly  than  in  classical 
control  theory.  That  is,  the  t r ia I -and-error  analysis  and  reanalysis  calculations  are  at  a minimum.  It 
is  first  necessary  that  equation  (I)  and  the  desired  performance  criteria  be  expressed  as  responses,  r^, 
in  the  time  domain.  Before  the  gains,  K*,  may  then  be  determined,  an  index  of  performance,  J(r^,t),  is 
mininized  to  determine  the  optimal  gains,  K,  subject  to  constraints  on  the  states,  and  controls,  In 

CCV  synthesis  studies  sponsored  or  conducted  by  the  AFFDL,  the  specific  form  of  J most  commonly  used  has 
been  the  quadratic  performance  index: 

J(r,r)  * F r,r  _tT  n HrT 
~ 0 

where  Q is  a weighting  matrix  selected  a priori  and  E is  the  expectancy  operator  required  for  systems  with 
nondeterministic  inputs,  (t).  More  specifically,  the  process  employed  to  determine  optimal  gains,  K, 
that  relate  the  controls  £ to  the  states,  ><,  involves  an  iterative  solution  for  the  Lagrange  multipliers, 
P,  using  the  matrix  Ricatti  equation  expressed  by: 

0 - A^P  + PA  + Q - PEP 

where  A - A - B (D^QD)"'  O^QC 
i - B (D^QO)''  B^ 

Q - c’^QC  - C^QD  (D^QO)"'  O^QC 


then,  knowioQ  P^, 


and  Final ly , 


K * -(D^QO)’^  [D^QC  + B^P] 


K X 


Only  in  the  simplest  systems,  e.g.,  rigid  aircraft,  is  it  possible  to  measure  the  states  x as  feed- 
back. In  general  the  measurements,  include  a combination  of  mean-axis-system  motions  as  weTl  as 
structural  deform.it  ions  and  for  this  reason  a second  set  of  gains,  K*,  using  the  realistic  measurements, 
are  then  calculated.  The  procedure  for  doing  this  varies  and  is  dependent  upon  the  experience  of  the 
flight  control  engineer.  One  method  used  at  the  AFFDL  consists  of  augmenting  the  measurement  equation 
with  states  not  directly  measureable,  and  then  Inverting  the  measurement  equation  to  arrive  at  a relation- 
ship between  states,  x,  and  augmented  measurements,  . This  new  relationship  Is  substituted  in  the 
response  equation  and  the  index  of  performance  in  order  to  arrive  at  the  second  optimization  problem. 
Finally,  an  Incremental  gradient  aiggrithm  based  upon  the  previously-determined  weighting  matrix,  Q,  and 
optimal  gains,  K,  is  applied.  The  process  of  arriving  at  the  best  practical  gains,  K* , is  partially 
trial-and-error  due  to  numerical  stability  problems  that  develop  if  the  incremental  gradient  method  begins 
too  far  from  the  global  minimum  of  the  problem. 


The  advantages  of  optimal  control  are  substantial  since  a stable,  mul tiple-input/mul tiple-output 


CCV  control  system  can  be  designed  directly  In  the  optimal  sense  with  a relatively  small  manpower  Invest* 
ment.  This  means  a particular  component  of  the  responses,  such  as  wing  root  bending  nKwient,  can  be 
treated  without  aggravating  other  responses  such  as  those  associated  with  handling  qualities  or  flutter. 

One  of  the  disadvantages  of  optimal  control  methods  is  that  it  is  difficult,  in  a stable  numerical  sense, 
to  enforce  state,  control  and  moasuretiiont  constraints  except  by  a trlal-and-crror  method  involving  Q. 

Also  it  is  necessary  that  the  unaugmented  system  be  stable  at  the  beginninq  (usually  easy  to  accomplish) 
and  it  is  not  generally  possible  to  enforce  relative  stability,  frequency  domain  stability  margins  even 
though  the  system  is  guaranteed  stable  in  an  absolute  sense.  Additionally,  all  aerodynamic  information 
must  be  expressed  In  the  time  domain. 

The  importance  of  maintaining  consistency  In  the  aerodynamics  and  dynamics  formulations  is  clearly 
Illustrated.  The  optimal  or  suboptimal  control  system  that  evolves  depends  upon  all  the  information  pre- 
sented In  the  state,  measurement  and  response  equations.  If  the  flight  control  specialist  accomplishes 
Its  synthesis,  6_  * K*^,  using  one  set  of  numerical  values  for  A,  B,  C,  0,  and  M,  and  the  structural 
specialist  then  uses  other  numerical  values  to  analyze  this  CCV  aircraft  according  to  MILSPEC  require- 
ments, the  resulting  design  is  sure  to  be  erroneous. 

3.0  aerodynamic  FORMULATIONS 

The  riqht-hand  sides  of  equations  (la),  (Ib),  and  (Ic)  contain  terms  associated  with  aerodynamic 
forces  that  act  upon  the  vehicle.  These  forces  are  known  to  be  proportional  to  the  mean-axis-system 
motion,  x^,  elastic  deformations,  controls,  and  the  atmosphere,  There  are  numerous  theoreti- 

cal, semi empi r i ca 1 and  exper (mental  approximations  of  those  aerodynamic  forces  reported  in  the  literature. 
The  Woodward  method,  vortex  lattice,  strip  theory,  doublet  lattice,  Rho's  method,  kernel  function,  piston 
theory,  spline  aerodynamics,  etc.,  are  just  some  of  the  highly  specialized  aerodynamic  methods  that  are 
employed  in  CCV  analysis  and  design.  Of  these,  the  steady  aerodynamic  methods  employed  to  calculate 
Initial  values  of  the  states,  measurements  and  responses  do  not  at  present  appear  to  introduce  consistency 
problems  because  a great  deal  of  experimental  data  and  user  experience  have  been  accrued. 

On  the  other  hand,  theoretical  computation  of  the  unsteady  aeoodynamic  forces  that  are  very 
important  to  the  CCV  dynamic  analysis  and  synthesis  problem  appears  to  present  the  designer  with  two 
choices  that  may  not  enforce  consistency  between  the  methods  of  the  structural  and  flight  control  staffs. 
With  one  choice,  the  unsteady  aerodynamics  are  estimated  in  the  frequency  domain,  resulting  in  frequency- 
dependent-coefficient,  ordinary  differential  equations  that  are  most  suitable  for  analysis  via  classical 
control  methods.  This  means  absolute  and  relative  stability  and  frequency  response  of  the  open-  and 
closed-loop  system  may  be  estimated,  thus  satisfying  many  of  the  structural  design  MILSPEC  requirements 
cited  in  Table  I.  With  the  second  choice,  unsteady  aerodynamic  forces  are  estimated  in  the  Laplace 
domain;  the  resulting  linearized  equations  are  then  transformed  to  the  time  domain  where  control 
system  analysis  and  synthesis  via  modern  control  theory  may  be  accomplished.  The  transforrat ion  of  these 
constant-coeff icient,  linearized  equations  from  the  Laplace  to  the  frequency  domain  Is  well  known.  How- 
ever, frequency-dependent-coefficient  equations  in  the  frequency  domain  may  not  be  transformed  to  the 
Laplace  and  time  domains  without  some  approximation.  The  approximation  involved  in  passing  from  frequency 
to  Laplace  and  time  domains  introduce  consistency  grob}€(r\s  between  the  methods  of  structural  and  flight 
control  specialists. 

One  means  of  identifying  potential  inconsistency  is  to  consider  recent  CCV-related  in  which 

the  aerodynamic  methodology  Is  isolated  and  its  impact  on  CCV  studied:  Four  recent  studies  cone  to  mind: 

• Cunningham  (Ref.  7)  investigated  the  problem  of  extending  frequency  dependent  aerodynamics  to  the 
Laplace  and  time  domain  of  analysis;  the  objective  was  development  of  an  unsteady  aerodynamic 
methodology  compatible  with  1 inear-systems  root  locus  analysis  and  nonlinear-system  time  history 
analysis.*  The  procedure  is  to  fit  f requency7dependent  generalized  aerodynamic  forces  with 
Tschebychev  polynomials,  then  Fourier  transform  to  the  time  domain  to  get  Indicia!  functions, 
approximate  the  Indiclal  function  by  a second  series  of  Tschebychev  polynomials,  and  finally 
Laplace  transform  that  series.  The  result  Is  a polynomial  with  terms  having  first  order  poly- 
nomials In  the  Laplacian  variable  as  their  denominators.  One  problem  arising  from  this  procedure 
Is  that  a large  number  of  poles  are  Introduced  for  each  retained  dynamic  degree  of  freedom;  for 
Instance,  in  one  application  a total  of  68  poles  result  when  three  mean-axis  degrees  of  freedom  and 
one  structural  mode  degree  of  freedom  are  retained  and  unsteady  aerodynamic  force  indicial  functions 
are  approximated  by  a five  term  polynomial.  The  large  number  of  poles  Introduces  numerical  error 
when  Che  roots  of  the  equation  are  determined  for  root  locus  studies;  the  practical  upper  limit  for 
acceptable  numerical  accuracy  is  100-150  complex  number  roots.  The  test  of  this  method  on  the 
YF-16  missiles-on  roll-loop  stability,  discussed  In  section  k.l,  indicate  that  the  root-locus 
classical  method  predicts  Instability  at  6.11  Hz,  and  at  a lower  roll-loop  gain  than  the  6.50  Hz 
Instability  observed  in  flight  test.  The  root  locus  analysis,  however,  matches  the  results  of 
previous  analytical  analyses  In  the  frequency  domain. 

• Roger  (Ref.  8)  approximated  the  frequency-dependent  aerodynamic  forces  of  the  doublet  lattice 
method  as  a rational  polynomial  function  of  the  Laplace  variable  subject  to  the  restrictions  that 
(I)  it  must  have  complex-conjugate  symmetry;  (2)  the  denominator  roots  must  lie  in  the  left  half 
of  the  complex  plane;  and  (3)  it  must  approximate  the  values  of  the  complex  coefficient  when 

s ■ 0 ♦ jw  for  those  values  of  ui  analyzed  using  the  doublet  lattice  theory.  The  6-52E  CCV  program 
analysis  and  synthesis  tasks  employ  this  method.  Another  paper  In  the  conference  discusses  the 
method  In  more  detail. 

• Brune  (Ref.  9)  modified  the  steady  aerodynamic  method,  developed  at  the  Boeing  Company  by  Woodviard 
et  a1,  to  Include  a first-order,  reduced-frequency,  unsteady  aerodynamic  approximation  to  be  used 
in  the  frequency,  Laplace  and  time  domain.  This  approximation  is  applicable  to  general  three- 
dimensional  configurations  at  subsonic  and  supersonic  speeds  (such  as  swept-wing  aircraft  Idealized 
In  Figure  1).  The  method  Is  theoretical  ly  valid  for  reduced  aerodynamic  frequencies,  k ■ <a>c/21>), 
of  order  much  less  than  1.0,  much  less  than  C(I-m2)/m^3  at  subsonic  speeds,  and  much  less  than 

at  supersonic  speeds.  In  contrast,  doublet  lattice  and  other  unsteady  aerodynamic 
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methods  are  .theoret  Icnlly  valid  Tor  a reduced  frequency  of  the  order,  !T(I'M'^)/M^  at  subsonic  speeds 
and  of  the  order  at  supersonic  speeds.  A numerical  quant i t I zat Ion  of  these  boundaries, 

as  well  as  a brief  summary  of  the  Woodward  approach  to  the  aerodynamics,  is  presented  in  Figure  2. 
When  pressure  distributions  estimated  by  doublet  lattice  theory  for  various  aspect  ratio,  wing- 
alone  configurations  are  compared  to  similar  data  estimated  by  this  low  frequency  approx imat ion , 
boundaries  of  appi icabt I i ty  may  be  established  (Ref.  9) • These  boundaries  for  aspect  ratio  (AR) 

2 and  8 wings  at  subsonic  Mach  numbers  are  shown  in  Figure  2.  Note  that  the  method  is  most  theo' 
rctically  valid  for  slender,  low  aspect  ratio  configurations.  If  the  reduced  frequency  required 
for  a 0 Hz  to  10  Hz  flight  control  system  for  the  C“9A  or  B-l  are  plotted  on  Figure  2,  several 
important  points  would  be  noted.  First,  from  a theoretical  standpoint  alone  it  would  appear  that 
most  methods  are  barely  adequate  at  subsonic  speeds  and  are  inadequate  at  transonic  speeds.  Also, 
the  low  frequency  approximation,  although  of  questionable  accuracy  at  subsonic  speeds  when 
compared  with  other  theories,  may  be  adequate  for  CCV  analysis  and  synthesis  at  supersonic  speeds. 
Here,  low  aspect  ratio  and  low  reduced  frequencies  are  the  rule. 

The  most  exhaustive  study  of  the  low  frequency  method  for  CCV  applications  has  been  on  an  aspect 
ratio  6.96  wing-alone  configuration  with  the  structure  and  (Mach  number  0.80)  aerodynamics  repre- 
sentative of  a largo  transport  aircraft.  This  study,  conducted  by  Kroll  and  Hiller  (Ref.  10) » Is 
limited  to  an  Investigation  of  frequency  characteristics  of  unsteady  aerodynamic  terms  In  transfer 
functions  that  relate  responses  and  measurements  to  vertical  turbulence  and  Inboard  aileron 
harmonic  motion.  The  data  are  presented  as  power  spectral  density  plots  and  I lnear*'ampl Itude  Code 
plots  of  the  unaugmented  aircraft  transfer  functions.  The  reference-axis-system  motion  states, 
are  reduced  and  then  transformed  to  be  the  inertial  coordinate  states:  vertical  displacement, 

z,  and  rotation  about  the  center  of  mass,  0.  Six  invacuum  modes  are  retained  In  the  states, 

In  this  truncated  modes  formulation  the  Initial  conditions  are  wlngs-level  flight,  parallel  to  the 
flat  earth;  also,  gravity  terms  are  neglected.  The  aerodynamic  methods  employed  are: 

1 . Woodward  method 

2.  Woodward  method,  modified  by  Kussner/Wagner  C'i'(w)/^(u})]  functions. 

3.  Woodward  method,  modified  by  low  frequency  corrections  on  the  states 

4.  Woodward  method,  modified  by  low  frequency  correction  on  states  Qusts 

and  controls 

5.  Strip  theory  method  corrected  with  Kussner/Wagner  f'f'(w)/^(w)!]  functions. 

6.  Rho  kernel -funct ion  unsteady  aerodynamic  method. 

7.  Vortex  spline  kernel-function  unsteady  aerodynamic  method. 

8.  Doublet  lattice  unsteady  aerodynamic  method. 

Figure  3 contains  selected  results  of  the  study  for  flight  in  vertical  turbulence.  In  this 
analysis  of  wing-tip  vertical  acceleration,  a2,  and  wing-root  bending  moment,  the  vortex  spline, 
doublet  lattice,  and  Rho  methods  provide  nearly  identical  results;  thus,  the  comparisons  are  not 
presented  here.  However,  when  the  various  formulations  of  the  Woodward  methods  and  the  strip 
theory  method  are  compared  with  the  Rho  method,  the  differences  are  moderate  to  large,  with  the 
largest  discrepancies  occurring  at  frequencies  greater  than  2 Hz.  Possible  reasons  for  this  are  a 
poor,  low  frequency  approximation  of  the  more  exact  unsteady  aerodynamic  theories,  or  the  intro- 
duction of  inaccuracies  when  the  Kussner/Wagner  functions  are  applied  to  the  states  of  a three- 
dimensional  wing.  The  results  of  the  report  concerning  control  surface  transfer  functions  should 
be  reviewed  with  some  care,  as  the  Woodward  methods  were  not  carefully  "tuned"  to  control  power 
calculations  as  the  doublet  lattice  method,  and  because  the  strip  theory  method  includes  empirical 
corrections  while  the  other  methods  do  not. 

• Disney,  Hargrave,  and  Hollenbeck  (Ref.  II)  analyzed  and  synthesized  in-part  the  Active  Lift  Dis- 
tribution Control  System  (ALDCS)  for  the  C-5A  using  a semi-empirical,  unsteady  aerodynamic  strip 
theory  method.  Good  correlation  with  flight  test  data  is  achieved  and  indicated  by  the  Bode 
amplitude  plot  (Fig.  4)  of  the  aileron  open-loop  transfer  function  of  the  ALDCS  aircraft.  Schv/anz 
and  Stockdale  (Ref.  12)  employed  the  low  frequency  correction  to  the  Woodward  method  to  analyze 
this  same  transfer  function  at  another  C-5A  mass  distribution;  the  results  are  compared  v/ith  the 
strip  theory  method  in  Figure  5.  Here  it  is  seen  that,  reasonable  correlation  in  phase  and  amplitude 
exists  up  to  4.5  Hz  (k^O.5);  above  4.5  Hz  the  two  methods  approximate  structural  dynamics  differ- 
ently. The  computer  method  used  by  Schwanz  and  Stockdale  is  the  FLEXSTAB  program  (Ref.  9);  the 
aileron  unsteady  aerodynamic  forces,  6 (w)  , are  neglected,  but  care  has  been  taken  to  correct  the 
aerodynamic  forces  of  the  control  surTaces  to  match  experimental  data. 

From  this  discussion  it  may  be  inferred  that  an  unsteady  aerodynamics  method  that  is  theoretically 
exact  in  the  Laplace  domain  of  analysis  would  help  to  enforce  consistency  between  the  methods  of  the 
structural  and  flight  control  analysts.  To  be  of  greatest  utility,  the  new  theory  should  possess  charac- 
teristics permitting  exact  trans format  ions  of  the  linearized  form  of  equation  (I)  between  the  frequency, 
Laplace  and  time  domains  of  analysis.  Until  a solution  is  found,  flight  control  engineers  must  continue 
to  question  and  test  the  unsteady  aerodynamic  approximations  they  employ  In  design,  and  structural 
spe'  ^alists  must  continue  to  appreciate  the  role  of  unsteady  aerodynamics  in  Che  special  synthesis  problem 
that  the  flight  control  specialist  solves  for  the  CCV  aircraft. 

49  DYNAMICAL  SYSTEM  FORMULATION 

ixperi^nce  wlthlr  r,*e  Flight  Control  Division  at  the  AFFOL  has  Indicated  that  the  crux  of  the  incon- 
sistency problem  In  dynamic  system  description,  when  the  aerodynamic  formulations  are  consistent.  Is  the 
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formulation  and  approx imat Ion  of  the  equations  of  ntotion  contained  in  equation  (la).  These  equations  of 
motion  relate  the  states,  and  to  measurements,  responses,  controls  and  Inputs.  In  order  to 

subsequently  explore  these  possible  inconsistencies  in  the  formulation  and  application  of  the  equations 
of  motion,  i^  is  best  to  express  them  in  linearised,  second'order , ordinary  differential  equation  form. 
This  form  separates  aerodynamic  terms  from  the  Inertial,  dampinq  and  stiffness  terms.  Using  the  non* 
inertial,  n>ean  axis  as  the  body-fixed,  reference  axis  system  and  the  notation  of  References  3,  13  *Jnd  |l«, 
these  equations  may  be  written  as: 

M*ir  + -I’  ** 

mu-fdu  + ku  - x^,  i.  M|  • (2b) 
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Here  the  subscript  I denotes  initial  value,  while 
c0| , is  cosine  0|,  etc. 
s0| , Is  sine  0| , etc ■ 

H|  Is  the  total  mass. 

Ijjl  is  the  total  moments  and  products  of  inert! 
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f is  the  matrix  of  mean-axis  system  modes. 

^ Is  the  matrix  of  Invacuum  modes. 

f^p  are  the  distributed  aerodynamic  forces  due 
and  turbulence. 


I about  the  non-inertial  axis  system. 

Euler  angles  in  the  rotation  sequence  of  Reference 

components  of  axis-system  translation, 
components  of  axis-system  rotation  rate. 

to  changes  in  state,  controls  or  atmospheric  gusts 


I).  I Reduction  in  Elastic  Deformation  States 

An  examination  of  equation  (2)  indicates  than  it  can  consist  of  a very  large  number  of  second-order 
equations  for  y.  Common  sense,  augmented  by  theoretical  and  numerical  accuracy  considerations,  dicta'es 
that  equation  (2)  should  be  reduced  in  size  to  its  absolute  minimum.  Typically,  in  this  reduction  process 
the  analyst  first  determines  the  frequency  range  over  which  the  CCV  control  system  will  be  required  to 
solve  design  problems,  e.g.,  from  0 Hz  to  20  Hz.  In  this  example,  the  problem  of  state  reduction  then 
simplifies  to  eliminating  higher  frequency  states  while  maintaining  reasonably  correct  dynamics  from  0 Hz 
to  20  Hz  excitatjr.1*  frequencies,  and  appropriate  asymptotic  behavior  as  excitation  frequency  exceeds 
20  Hz  and  approaches  infinity.  There  are  three  state  reduction  procedures  reported  in  the  literature: 

Residual  flexibility  - the  invacuum  mode  shapes  representing  tne  elastic  motion  are  separated  into 
"retained"  and  "deleted"  modes.  The  deleted  modes  are  purged  dynamically  but  retained  stati- 
cally as  static  elastic  corrections  to  the  remaining  states.  The  correction  factors  are  related 
to  the  retained  invacuum  modes  and  the  flexibility  matrix  or  stiffness  matrix  of  the  structure. 

Truncated  modes  - the  deleted  modes  of  the  residual  methods  are  not  represented  by  the  static 
elastic  correction  factor.  This  is  the  most  common  method  reported  in  the  literature. 

Static  elastic  or  quasi-static  - Che  motions  of  the  structure  relative  to  the  mean  axis  are  assumed 
to  be  in  phase  with  the  axis  system  motions. 

Only  recently  have  these  state  reduction  procedures  been  tested  numerically  on  realistic  CCV  hardv/are 
design  cases: 

* Schwanz  and  Stockdale  (Ref.  12)  studied  the  C-5A  ALDCS  design  problem  using  equation  (2)  as  mech- 
anized in  the  FLEXSTAB  computer  software  of  Reference  9-  Figure  6 presents  the  effects  of  resi- 
dual ization  and  truncation  on  the  short  period  and  phugoid  dynamics  of  the  unaugmented  C-5A  air- 
craft. Actual  data  points  are  noted  at  3,  7,  and  l'3  retained  structural  modes;  the  straight  lines 
are  drawn  only  to  facilitate  visualization  of  the  trend.  As  shown,  the  residual  formulation 
predicts  the  mean-axis-system  dynamic  characteristics  nearly  independently  of  the  number  of  modes 
retained  while  the  truncated  formulation  does  not.  Also  note  Chat  the  residual  formulation 
correctly  coverages  to  the  static  elastic  formulation  for  zero  retained  modes,  while  the  truncated 
formulation  does  not.  The  reason  for  the  improved  accuracy  is  that  in  the  residual  formulation, 
the  generalized  aerodynamic  forces  as  well  as  the  mean-axis-system  stability  derivatives  are 
numerically  dependent  upon  the  number  of  modes  retained,  whereas  in  the  truncated  formulation 
these  assume  a constant  value  independent  of  the  number  of  retained  modes.  As  as  example.  Figure  7 
presents  the  stability  derivatives  that  are  used  in  the  3,  7,  and  13  retained  structural  mode 
study  Summarized  in  Figure  6. 

* Konar,  Mahesh,  Stone,  and  Hank  (Ref.  15)  investigated  the  effect  of  res idua I i za t ion  and  truncation 
upon  the  design  of  an  optimol  ALDCS  control  system  fpr  the  C-5A.  Of  interest  here  is  the  effect 
of  residual  and  truncation  mathematical  formulations  upon  the  time-domain  responses  of  the  closed- 
loop  dynamical  system.  These  data  are  calculated  by  first  updating  the  FLEXSTAB  math  model  of  the 
C-5A  (Ref.  12)  with  experimental  data  to  match  flight  test  results  as  closely  as  possible.  Then, 
an  optimol  ALDCS  control  system  is  designed  (as  discussed  in  section  2.0)  for  the  residual  flexi- 
bility case  in  which  2k  unougmented  aircraft  states  are  retained.  Finally,  given  that  control 
system,  time  histories  of  selected  responses  are  calculated  for  the  case  in  which  only  the  sensor 
equations  and  then  the  sensor  plus  state  equations  are  truncated.  Typical  responses  are  presented 
In  Figures  8 and  9.  As  shown,  the  pitch  rate  of  the  center  of  mass  and  the  wing-root  bending 
moment  are  more  sensitive  to  sensor  truncation  than  to  state  truncation. 

* Peloubet,  Haller,  Cunningh.am,  Cwach,  and  Watts  (Ref.  7)  analytically  studied  the  effect  of  resi- 
duallzation  and  truncation  on  a flight-test  observed  instability  in  the  YF-16  handling  qualities 
augmentation  system.  The  problem,  mentioned  previously  in  section  3-0  was  corrected  In  the  YF-lb 
In  197^,  but  is  of  continuing  research  interest  because  it  is  difficult  to  thecret ica 1 1 y predict 
the  phenomena.  Tables  2 and  3 present  selected  results  from  Reference  7-  Here  the  analysis  Is  In 
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the  frequency  dom.nn  with  experimental ly-corrccted,  doublet  lattice  theory  providinq  the  unsteady 
aerodynamics  Information.  As  shown  In  Table  2,  the  low  frequency  (u»  * 0.153  Hz),  residual  flexi- 
bility measurements  for  sensors  near  the  center  of  mass  of  the  unauqntented  aircraft  have  a magni- 
tude and  phase  angle  nearly  Independent  of  the  number  of  retained  structural  degrees  of  freedom. 

In  contrast,  the  truncated  iiH)dcs  formulation  provides  transfer  function  amplitudes  and  phases 
that  are  dependent  upon  the  number  of  retained  states.  In  the  truncated  formulation  of  the  un- 
augmented  YF-16  dynamics,  It  takes  six  clastic  degrees  of  freedom  to  approach  the  residual 
flexibility  formulation  results  for  only  one  retained  structural  degree  of  freedom.  When  the 
lateral-directional  handling  qualities  augumentat ion  system  is  included  In  the  analysis  the 
results  are  somevd^at  mixed.  In  the  case  of  yaw-loop  closed  and  roll-loop  open  for  the  wing-tip 
missiles-on  (Table  3a),  the  residual  formulation  solution  for  the  phase-crossover  gain-margin 
converges  much  more  rapidly  than  does  the  truncated  modes  formulation.  Based  upon  convergence 
alone,  the  missiles-off  conf Igurat Ion  (Table  3b)  also  appears  to  be  bes-t  represented  by  the  resi- 
dual flexibility  formulation. 

From  these  results,  it  may  be  concluded  that  the  truncated  and  the  residuallzed  forms  of  equations 
(1)  and  (2)  describe  different  system  dynamics  for  the  unaugmented  and  augmented  aircraft.  Of  course, 
the  difference  between  the  two  formulations  may  be  large  or  small  depending  upon  the  aircraft  and  its 
mission.  A review  of  the  literature  indicates  that  the  structural  spcalalist  currently  perfers  the 
truncated  formulation,  while  the  flight  control  specialist  may  employ  both  formulat ions . Clearly,  incon- 
sistencies may  arise  if  each  discipline  approaches  the  design  requirements  with  Its  own  formulation  of 
the  dynamics. 


A. 2 Representation  of  Initial  Conditions 

It  is  standard  practice  in  the  flight  control  studies  to  represent  the  dynamics  of  the  mean-axis- 
system  by  states  and  as  indicated  in  equation  (2).  In  contrast,  structural  specialists  most  often 
approximate  these  states  as  the  so-called  "rigid  body  modes": 


[xz0yij*jj  , where  x,  y,  and  z are  Inertial  coordinates. 


An  inspection  of  the  terms  in  equation  (2)  Indicates  that  initial  conditions  of  motion  (quantities 
that  a“e  subscripted  with  "I")  play  a prominent  role  in  the  numerical  calculat  ions . Of  further  interest 
then  is  the  consideration  of  initial  conditions  in  state  equations  with  independent  coordinates  z^. 

This  consideration  of  intitial  conditions  is  best  approached  by  observing  the  kinematic  expression 
in  equation  (2c)  provides  a partial  transformation  between  ^ and  namely  (Ref.  \U): 
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n 0 -s*  cOc» 

If  equation  (2c)  is  augmented  by  an  additional  kinematic  expression  between  u,  v,  w and  x,  y and  z,  a 
complete  expression  may  be  developed  to  transform  (2)  from  and  ^ space  to  and  space.  This 

additional  kinematic  expression  is  the  well  known  trajectory  equations  used  in  stability  an^  control 
("flight-path"  equations  of  Reference  14): 
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These  kinematic  expressions  may  be  combined  to  yield  the  desired  transformation  of  state  (Ref.  16): 
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Note  that  the  C runs  forms t ion  contains  both  initial  conditions  of  motion  as  well  as  perturbation  values. 
Once  the  initial  conditions  arc  specified  in  space,  the  general  form  of  the  transform  becomes: 

X - Tlz.  + T2z  (4) 


where  Tl  and  T2  may  be  nonlinear  in  the  kinematic  states,  x . 

As  examples: 

• The  transformation  for  the  case  of  steady,  straight,  wings-ievel,  climbing  flight  is  defined  using 
stability  axis  (W|-0)  variables  by  setting  Vj , Pj , , Rj,  f,  equal  to  zero  In  equation  (3); 
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U|  and  0|  are  nonzero  constant  values.  In  this  case  Tl  and  T2  are  independent  of  ^ and  sparse: 
u - X 

w - xsO|  + zcO|  + U|0 
q - 6 

V • y-i(i  U|C0|  (5) 
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• The  transformation  for  the  case  of  steady,  wings-level,  Ig  straight  and  level  flight  has 
only  U|  nonzero.  Thus,  for  stability  axes: 

u “ X 

w = 2 + U|6 

<1-0  (6) 

V = y - U|i(/ 
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Conventional  and  CCV  flight  control  systems  are  designed  to  provide  satisfactory  performance  for 
many  initial  conditions  besides  steady,  level,  Ig  climbing  flight.  This  is  accomplished  using 
equations  of  motion  such  as  indicated  in  equation  (2).  However,  it  has  been  the  convention  of  many 
structural  dynamiclsts  to  denionstrate  compliance  with  certain  of  the  MILSPECS  using  equations  of  motion 
in  and  u space  that  apply  only  to  steady  Ig  straight  and  level  flight.  This  formulation  can  be  visu- 
alized by  first  equating  all  initial  conditions  in  equation  (2a)  to  zero  except  Uj  and  then  transforming 
using  equation  (6).  It  should  be  clear  that  the  solution  of  equations  (2)  for  the  states,  measurements 
and  responses  of  the  unaugmented  aircraft  will  be  different  for  flight  control  and  structures  staffs 
unless  great  care  is  taken  to  account  for  the  initial  conditions  in  a consistent  fashion. 

Consider  the  importance  of  initial  conditions  on  the  horizontal  tail  torsion  due  to  a frequency 
sweep  by  the  outboard  aileron  as  presented  in  figure  10.  These  data  presented  in  the  form  of  Bode 
logarithmic  amplitude  and  linear  phase  angle  plots,  are  estimated  theoretically  for  an  unaugmented 
B-52/L7  aircraft  configuration  using  the  Woodward  method  corrected  for  the  low  frequency  aerodynamics  on 
the  state  only.  The  FLEXSTAB  computer  program  defined  in  Reference  9 is  the  computer  mechanization  of 
equation  (2).  The  transfer  function  that  is  presented  is  typical  of  a response  that  might  be  considered 
in  designing  a maneuver  load  control  system  for  the  wing.  As  indicated  in  Figure  10,  the  horizontal  tall 
torsion  for  these  two  Initial  conditions  is  substantially  different  in  both  phase  and  amplitude  at  low 
frequencies  of  the  order  of  the  phugoid  where  a reduced  static  stability  or  handling  qualities  control 
system  might  be  required. 

4.3  Forward-Speed  Degree  of  Freedom 

It  is  often  apparent  in  structural  analyses,  particulatly  those  associated  with  flutter,  that  the 
structural  dynamicist  has  neglected  the  forward-speed  degree  of  freedom,  u (or  x) , while  the  flight  control 
specialist  has  not.  The  neglect  of  forward-speed  effect  means  that  the  first  equation  In  (2a)  is  elimi- 
nated. This,  affecting  consistency  of  analysis  between  structural  and  flight  control  specialists,  was 
justified  in  the  past  on  the  basis  that  there  Is  a large  frequency  separation  between  the  phugoid  and 
short  period  modes;  it  was  suff Iclent^for  the  structural  dynamicist  to  represent  only  the  short  period 
dominated  motions,  w and  q (or  z and  6).  However,  reducing  the  static  longitudinal  stability  can  cause 
a coalescence  of  short  period  and  phugoid  motions  into  two  aperiodic  modes  plus  a third  oscillatory  mode, 
a phenomenon  that  can  not  be  realistically  approximated  if  the  speed  terms  are  neglected. 

In  order  to  illustrate  numerically  the  importance  of  the  speed  degree  of  freedom,  a root  locus  on 
Cfn  (the  measure  of  longitudinal  static  stability)  of  an  unaugmented  B-52E/L7  configuration  Is  calculated 
for  the  Ig  wings-level  flight  case  discussed  in  Figure  10.  In  Figure  II  the  data  are  presented  for  a 
variation  of  from  -0.017/deg.  (statically  stable)  to  +0. 017/deg.  (statically  unstable).  Retention  of 
the  forward-speed  degree  of  freedom  results  In  the  presence  of  the  phugoid  mode.  As  is  varied  this 

mode  first  moves  to  the  real  axis  and  then  coalesces  with  the  aperiodic  short  period  mode  to  form  the 
third  oscillatory  mode  and  an  unstable  aperiodic  mode  with  a very  small  t ime-to-double  amplitude.  If  the 
forward-speed  degree  of  freedom  Is  neglected,  the  coalescence  of  modes  to  form  the  third  oscillatory  mode 
would  not  occur  and  the  unstable  short  period  aperiodic  mode  would  possess  a different  dynamic  character- 
istic. A reduced  static  stability  or  other  handling  qualities  control  system  is  of  course  sensitive  to 
ail  these<modes.  It  is  also  Influenced  by  and  alters  other  higher  Frequency  dynamic  modes  owinq  to 
(1)  structural  feedback  through  the  control  surface  actuator;  (2)  sensor  measurement  of  higher  frequency 
structural  modes;  and  (3)  eigenvector  coupling  between  states  of  the  compensating  control  system  and 
invacuum  structural  modes. 

5.0  CONCLUDING  REMARKS 

The  Intent  of  design  guidelines  to  enforce  consistency  In  structural  and  flight  control  analysis  and 
synthesis  has  been  discussed.  As  Indicated,  CCV-typc  control  systems  are  mu! t iple- 1 nput/mul 1 1 pie-output 
systems  in  most  applications,  requiring  that  the  aerodynamic  formulation  be  expressed  In  the  frequency, 
Laplace  and  time  domains  In  order  to  demonstrate  compliance  with  MILSPECS  and  similar  guidelines.  Also, 
as  shown,  the  mathematical  formulation  describing  the  dynamical  system,  particular! ly  the  equations  of 
motion  of  unaugmented  aircraft,  must  not  possess  unjustifiable  or  unproven  approximations  and  must  define 
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the  system  from  very  low  to  moderntely  large  frequencies.  In  order  to  enforce  consistency  in  structural 
and  flight  control  analyses  It  is  apparent  from  the  examples  given  that  the  flight  control  specialist 
should  become  increasingly  familiar  with  unsteady  aerodynamic  theories  and  computer i red  methods.  In 
addition,  it  is  evident  that  the  structural  specialist  should  become  equally  familiar  with  the  more 
precise  system  dynamical  formulations  that  the  flight  control  specialist  uses  mat ter-of-fac t I y . 


Ingredients  for  good  and  appropriate  criteria,  one  goal  of  the  Conference  Committee  on  Structural 
Aspects  of  Active  Controls,  may  well  be  found  in  the  state  space  form  of  the  equations  describing  the 
controlled  dynamical  system,  c.g.,  equations  (la  - Id).  This  formulation,  expressed  in  the  time  domain, 
enforces  consistency  between  flight  control  and  structural  specialists'  analyses  by  clearly  defining 
technical  responsibility  for  each  state,  measurement,  response  and  control  equation.  The  solution  of 
these  equations  provides,  simultaneously,  direct  and  indirect  numerical  evaluations  of  significant 
parameters  that  are  of  interest  to  the  specialists  who  must  currently  interpret  the  existing  MILSPECS  and 
other  guidelines.  Thus,  the  state-space  equations  are  a numerical  realization  of  the  design  guidelines 
that  interrelate  separate  disciplines  in  the  CCV  analysis  and  synthesis  task. 
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6 . 6 Effects  of  norop  I.Vi  1 1 c i 1 Y ■ control  cqiii  pniciU  , .iml  s tructiir.il  dyn.Tiiics.  Jl.  3-5.  }•? 

Since  acroe  las  t i c i ty , control  eguipment,  anil  structural  dyn.nuics  m.iy  exert 
an  Important  Influence  on  the  airplane  flyinq  qualities,  such  effects  should 
not  be  overlooked  in  calculations  or  analyses  directed  toward  Investigation 
of  compliance  with  the  requirements  of  this  specification. 


3.5  De forma t ions ■ The  cumulative  effects  of  elastic,  pcrmcnent,  or  thermal  3-8,  8.2 
deformations,  acting  singly  or  together,  which  result  from  application  of 
landing  loads,  fatigue  loads,  and  limit  loads  shall  not: 

a.  Inhibit  or  degr.ide  the  mechanical  operation  of  the  airplane 

b.  Affect  the  airplane's  aerodynamic  characteristics  to  the  extent 
that  performance  guarantees  or  flying  qualities  requirements  cannot 
be  met 


c.  Require  repair  or  replace  of  parts. 


3.1k  Stability  augmentation  devices.  The  effect  of  stability  augmentation 
and  automatic  control  devices  shall  be  included,  where  applicable.  The 
design  conditions  shall  apply  for  operative,  inoperative,  and  transient 
modes.  In  failure  cases,  any  variation  or  reduction  of  the  device  capability 
shall  be  included  in  the  structural  design  condition. 

3.22.2.1'  Airplane  turbulence  response.  Airplane  turbulence  response  shall 
be  characterized  by  the  response  parameters  A and  Nq  which  are  determined  by 
a dynamic  analysis  of  the  airframe.  The  dynamic  analysis  shall  incorporate 
rigid  body  motion,  significant  flexible  degrees  of  freedom,  the  flight 
control  system,  and  the  stability  augmentation  system.  Augmentation  system 
requirements  shall  be  consistent  with  MIL-F-8785.  Significant  effects  shall 
be  represented  by  equivalent  linear  representation  which  shall  be  demon- 
strated to  be  conservative.  In  conjunction  with  the  dynamic  analysis,  when 
specified  in  the  application  contract,  a flight  simulation  shall  determine 
the  effects  of  control  inputs  on  structural  loads  when  encountering  extreme 
value  turbulence.  The  turbulence  levels  shall  be  of  sufficient  severity  to 
cause  structural  limit  load  occurrences.  Adverse  coupling  shall  be 
corrected.  The  following  definitions  shall  apply. 


3.1.11.2  St  I f fness . The  stiffness  of  flight  control  systems  shall  be 
sufficient  to  provide  satisfactory  operation  and  to  enable  the  aircraft 
to  meet  the  stability,  control,  and  flutter  requirements  as  defined  in  the 
applicable  portions  of  MIL-F-8785,  MIL-A-8870,  MIL-F-83300  and  MIL-A-8865. 
Normal  structural  deflections  shall  not  cause  undesirable  control  system 
inputs  and  outputs. 


3.2.1  Aur.ientatton  systems.  For  airplanes  with  augmentation  systems,  the 
flutter  margins  and  damping  requirements  of  3.1  shall  be  met  both  with  the 
system  inoperative  (system  off)  and  with  the  system  operating.  In  parti- 
cular, at  speeds  up  to  V|_,  the  operating  system  shall  be  stable  with:  (I)  a 

gain  margin  of  at  least  i6  dB,  and  (2)  separately,  a phase  margin  of  at 
least  ±60*. 


3.2.6,  3-2.7,  k.2. 
4.2.3,  4.2.6 


4. 1.1. 7 Servo-control  analyses.  The  dynamic  characteristics  of  control 
surface  actuating  systems  such  as  servo  boost,  fully  powered  servo  control, 
and  other  types  shall  be  included  in  the  flutter  analyses.  The  effect  of 
high  temperatures  on  the  dynamic  characteristics  of  the  actuating  systems 
including  the  hydraulic  fluid  shall  be  included.  Augmentation  systems 
which  may  alter  the  dynamic  response  of  the  airplane  shall  also  be  in- 
cluded in  the  flutter  analyses;  the  method  and  approach  used  in  these 
analyses  shall  be  subject  to  review  and  approval  of  the  procuring  activity. 


3. 1.1.2  Service  loads  spectrum.  The  service- loads  spectrum  is  derived 
from  a collection  of  load  spectra.  Each  loads  spectrum  in  this  collection 
shall  define  the  expected  (average)  number  of  load  cycles  according  to  load 
magnitude  for  a given  source  of  repeated  loads.  The  loads  spectrum  for 
each  significant  source  of  repeated  loads  shall  bo  based  on  a realistic 
interpretation  of  the  design  usage.  The  contractor  shall  include  all 
significant  sources  or  repeated  loads.  The  source  of  repeated  loads  may 
include,  but  not  be  limited  to,  ground  handling  and  taxiing  operations, 
landing  operations,  flight  maneuvers,  atmospheric  turbulence,  in-flight 
refueling,  autopilot,  inputs,  cabin  pressurization,  buffeting,  terrain- 
following, and  the  ground-to-ground  cycle. 


AEROPyWAMIC  CONFiGURATiOW 


AERODYNAMIC  COMPONENTS 


FIGURE  2.  Small  Perturbation  Expansion  Parameters  and  Theoretical  Range  of  Applicability 
of  Low  Reduced  Frequency  Unsteady  Aerodynamics  Method  (Ref.  9) 
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FIGURE  3.  Effect  of  Aerodynamic  Formulation  on  Wing  Tip  Acceleration  and  Wing'Root 

Bending  Moment  Response  Amplitude  Due  to  Vertical  Turbulence,  M-0.80  (Ref.  10) 
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FIGURE  k.  Correlation  of  Unsteady  Strip  Theory  with  C-5A  Flight  Test  Data, 
Aileron  Loop  Open,  Elevator  Loop  Closed  (Ref.  12) 
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FIGURE  5.  Correlation  of  Unsteady  Strip  Theory  Hethod  with  Low  Frequency 
Unsteady  Aerodynamic  Correction  to  Woodward  Hethod,  Aileron 
Loop  Open,  Elevator  Loop  Closed  (Ref.  12) 
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FIGURE  6.  Effect  of  Retained  Structural  Invacuum  Modes  on  C-5A 
Mean-Axis  System  Dynamics  (Ref.  12) 
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Pitch-Rate  Response  Due  to  Inboard  Elevator  Step  (Ref.  15) 
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TABLE  2.  YF-16  Unaugnicnted  Airplane  Sensor  Response  at  Lowest  Analysis  Frequency  (Ref.  7) 

Table  2a.  MIsslles-on,  M=>0.30,  Altitude  - 20000  FT. 
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TABLE  3.  YF-16  Roll  Rate/Alleron  Open-Loop  Frequency  Response  at  -I80  Degrees  Phase,  Yaw-Loop  Closed 

(Ref.  7) 

Table  3a-  Mlssiles-on,  M»0.90,  Altitude  = 20000  FT. 
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YC-lU  CONTROL  SYSTEM  REDUNDANCY 


William  T.  Hamilton 
Vice  President  - Engineering 
Boeing  Aerospace  Company 
P.O.  Box  3999,  M.  S.  85-86 
Seattle .Washington  9812U 


Summary 

The  YC-lU  is  the  Boeing  entry  in  the  USAF  Advanced  Medium  STOL  Transport  (AMST) 
program.  The  task  of  operating  a large  Jet  aircraft  into  and  out  of  a semi-prepared, 
2,000  feet  long  airstrip  with  a 27,000  pound  payload  presents  an  unusual  flight  control 
challenge.  The  YC-lU  answers  this  challenge  using  an  advanced  flight  control  system  that 
includes  digital  computers.  Excellent  STOL  flying  qualities  have  been  achieved  through 
control  wheel  steering  and  speed  hold  modes.  Fail  operational,  fail  safe  performance  is 
provided  by  a triplex  flight  control  system.  Aircraft  dynamics  following  an  engine 
failure  are  docile  and  do  not  require  Immediate  pilot  attention  or  unusual  skill.  The 
superior  capability  of  digital  computers  to  perform  logic  functions  enables  a comprehen- 
sive, semi-automated,  preflight  test.  Failures  are  detected  and  identified  to  the  Line 
Replaceable  Unit  (LRU).  The  YC-lU's  use  of  redundant  digital  computers  in  the  flight 
control  role  is  a first  for  an  aircraft  designed  to  demonstrate  operational  use. 


YC-lU  Configuration 

The  YC-lU,  as  shown  in  Figure  1,  incorporates  many  advanced  technology  and  design- 
to-cost  features.  Powered  lift  is  produced  by  directing  the  exhaust  of  two  large, 
modern,  hlgh-bypass-rat io  turbofan  engines  over  circular  cross-section  flaps  by  the 
Coanda  effect.  Studies  by  NASA  and  Boeing  have  confirmed  the  superior  performance  of 
Upper  Surface  Blown  (USB)  flaps  and  the  Coanda  effect  in  generating  powered  lift.  See 
References  1 and  2.  The  engines  are  located  close  to  the  fuselage  to  insure  adequate 
control  following  an  engine  failure.  A twin  engine  configuration  was  selected  because  in 
Boeing's  opinion  engine  procurement  and  maintenance  costs  are  lower  than  those  for  a four 
engine  aircraft,  and  a modern  engine  with  the  desired  thrust  level  was  available.  USB 
flap  deflections  can  be  modulated  at  high  rates  in  response  to  Control  and  Stability 
Augmentation  (CAS)  commands  for  airspeed  control.  Cruise  speeds  in  excess  of  M = 0.7 
with  an  unswept  wing  are  achieved  through  the  use  of  supercritical  airfoil  sections. 
Superior  rough  field  operation  is  attained  from  a long  stroke,  high  energy  absorbing, 
high  flotation,  lever  action  landing  gear.  Deslgn-to-cost  is  evidenced  by  the  unswept 
wing,  constant  sections  in  the  fuselage  and  vertical  tail,  and  identical  structure  in  the 
right  and  left  horizontal  stabilizer  panels'. 


Figure  1 Boeing  advanced  medium  STOL  transport  prototype  - YC-lU 


Concern  ie  often  expressed  regarding  a twin-engine  aircraft's  performance  with  a 
failed  engine.  Actually,  an  aircraft  with  two  engines  is  designed  to  have  the  same 
residual  performance  with  an  engine  failed  as  a four-engine  aircraft.  As  a result,  this 
twin-engine  aircraft  has  50  percent  more  Installed  thrust.  This  larger  thrust  installa- 
tion is  economically  sound,  since  engine  acquisition  and  ownership  costa  are  lower  for 
two  large  engines  than  for  four  smaller  ones.  The  normal  YC-llt  flight  control  system 
inherently  controls  the  dynamics  following  an  engine  failure  so  Immediate  pilot  action  is 
not  required;  nor  is  unusual  piloting  skill  needed.  In  addition,  advantage  was  taken  of 
the  flight  control  system's  ability  to  improve  STOL  performance. 


Key  design  flight  conditions  are  listed  in  Figure  2.  Short  takeoff  and  landing  (STOL) 
is  characterized  by  an  86  Kt  approach  speed  on  a 6°  glide  slope.  The  operating  lift 
coefficient  W > 3.6  is  over  twice  that  of  conventional  Jet  transports.  Cruise  at 


qS 

high  Mach  number  and  altitude  means  that  all  of  the  transonic  aerodynamic  phenomena  must 
be  considered  in  the  design.  Finally,  low-altitude,  high-speed  dash  imposes  requirements 
to  operate  at  high  dynamic  pressures.  The  ratio  of  maximum  to  minimum  operating  dynamic 
pressures  is  about  twice  that  of  conventional  aircraft  and  50  percent  more  than  supersonic 
transports . 


MACH.  NO.  JET  TRANSPORT  CRUISE 


86  KTS  lEAS)  LANDING 
6°  GLIDE  SLOPE 


Figure  2. 


YC-lU  key  design  flight  conditions 
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The  electrical  flight  control  system  centers  on  triple 
computers  (Ref.  1* ) . Sensor  inputs  and  servo  drive  outputs 
face  units.  All  cross-channel  data  used  for  redundant  oper 
optics.  Fault  detection  and  channel  Isolation  is  automatic 
Systems  Ltd.  of  Rochester,  Kent,  England  designed  and  built 
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The  aircraft  also  has  triple  digital  air  data  computers,  an  inertial  navigation 
system,  and  an  electronic  attitude  director  indicator  (EADI).  The  EADI  displays  attitude, 
flight  path,  altitude  and  airspeed  information  supplemented  by  background  real-world 
television  display.  Electronic  systems  are  used  extensively  to  enhance  YC-lU  performance 
(Ref.  5). 
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Figure  4,  Flight  control  surfaces 

From  Its  Inception,  Boeing  has  required  superior  flying  qualities  for  the  YC-lU. 
Powerful  controls  are  provided  to  enable  quick  aircraft  response  to  control  commands 
during  low  airspeed,  STOL  operation.  Since  most  aircraft  accidents  occur  during  takeoffs 
and  landings,  and  are  generally  related  to  flying  qualities  characteristics,  not  engine 
failure,  the  YC-l4  was  required  to  be  easy  to  fly  in  these  critical  flight  regimes.  The 
pilot  is  able  to  direct  his  attention  principally  to  outside  situations  inherent  in 
tactical  airlift. 


STOL  aircraft  typically  experience  powerful 
basic  stability  and  control  characteristics,  and 
thrust  curve.  This  is,  increased  thrust  is  requi 
condition  opposite  to  that  of  most  Conventional  T 
a pilot  attempts  to  control  the  flight  path  angle 
control  column,  flight  path  instability  results, 
craft  has  to  control  attitude  with  control  column 
While  the  backside  control  technique  is  feasible, 
concluded  at  Boeing,  not  consistent  with  effectiv 
YC-14  is  designed  to  be  controlled  with  normal  pi 
flap's  high  rate  performance,  the  YC-14  is  easier 
because  of  its  superior  airspeed  hold  mode.  The 
approach  is  simply  to  point  the  aircraft  to  the  d 
control  column  and  control  wheel  commands.  Airsp 
accurately. 


propulsive  lift  Interaction  with  their 
fly  on  the  "backside"  of  the  speed- 
red  to  fly  at  a lower  airspeed,  a 
akeoff  and  Landing  (CTOL)  aircraft.  If 
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Ordinarily,  a pilot  flying  a STOL  air- 
and  sink  rate  with  engine  throttles, 
the  control  task  is  demanding,  and  as 
e military  operation.  As  a result,  the 
lot  techniques.  In  fact,  with  the  USB 
to  control  than  contemporary  aircraft 
pilot's  control  task  during  landing 
esired  touchdown  point  using  one-hand 
eed  is  maintained  automatically  and 


Fail-operational/fall-safe  performance  was  required  to  provide  desired  mission 
reliability.  In  the  simulators  the  YC-l4's  STOL  flying  qualities  with  CAS  inoperative 
are  rated  Level  2,  Figure  5,  "Flying  qualities  adequate  to  accomplish  the  mission  Flight 
Phase,  but  some  increases  in  pilot  workload  or  degradation  in  mission  effectiveness,  or 
both  exists."  However,  without  CAS,  the  pilot  must  use  the  "backside"  control  technique 
in  the  critical  STOL  landing  flight  regime.  It  was  not  considered  prudent  to  require  this 
increased  pilot  workload  for  tactical  airlift.  With  fail-operational  systems,  CAS-off 
operation  becomes  so  infrequent  that  diversion  to  an  airfield  with  a longer  runway,  where 
a normal  CTOL  landing  can  be  made,  is  reasonable.  Aircraft  and  crew  safety  benefit  with- 
out unacceptable  degradation  of  mission  effectiveness. 

Fail-operational/fall-safe  performance  was  achieved  through  redundant  control  surfaces 
and  systems.  At  least  three  aerodynamic  control  surfaces  actuated  by  three  independent 
Hydraulic  Systems  are  available  for  each  control  axis.  Three  CAS  channels,  each  powered 
by  a different  Electrical  System,  are  also  provided  for  each  control  axis.  Mechanical 
elements  in  the  Flight  Control  System  (FCS)  enable  safe  controllability  with  a total 
Electrical  System  failure. 
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Figure  5. 


Pilot  ratings  of  YC-ll|  flying  qualities. 


STOL  landings 


Longitudinal  Flight 


Control 


Pilot  control  commands  are  transmitted 
as  shown  in  Figure  6.  During  flaps-down  op 
commands  are  transmitted  through  the  FCE  to 
plus  a SECS  model.  Since  no  off-the-shelf 
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Figure  6.  Longitudinal  flight  control  system  integration 
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compute  and  monitor  commands.  The  C-l4  will  use  the  same  serie 
flaps-down  flight.  Elevator  actuators  are  modified  7^7  compone 

Pilot's  commands  are  transmitted  electrically  to  wing  spoi 
USB  flaps,  used  to  control  airspeed  during  landing  approach  in 
powered  lift,  are  commanded  electrically  by  the  FCE  with  inputs 
and  aircraft  dynamics  sensors.  Electrical  commands  are  also  se 
airspeed  control. 
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The  all-movable  stabilizer  provides  longitudinal  trim.  7^7  trim  control  modules  and 
a 707  Jackscrew  are  used.  Automatic  trim  is  commanded  by  the  FCE.  When  the  FCE  is 
disengaged,  the  pilot  can  command  trim  electrically  by  a thumbswitch  on  the  control  wheel, 
or  mechanically  by  a level  located  on  the  aisle  stand. 


Lateral  Flight  Control 


Pilot  control  commands  are  transmitted  electrically  and  mechanically  to  ailerons  and 
spoilers,  as  shown  in  Figure  7.  Sensor  feedbacks  are  transmitted  through  the  FCE  to 
aileron  and  spoiler  SECS  during  all  flight  conditions.  Aileron  actuators  are  components 
used  on  the  E-2A  airplane.  Inboard  spoilers  and  their  actuators  are  7^7  units;  outboard 
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spoilers  and  their  actuators  are  from  the  727.  An  electric  actuator,  controlled  by- 
switches  on  the  aisle  stand,  recenters  the  feel  mechanism  for  trim. 
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Figure  7.  Lateral  flight  control  system  integration 
Directional  Flight  Control 

Pilot  control  commands  are  transmitted  mechanically  from  rudder  pedals  to  the  power 
control  actuators,  as  shown  in  Figure  8.  CAS  inputs  are  series  summed  at  the  actuators. 
Rudder  actuators  are  modifications  of  those  used  on  the  DC-10.  Trim  is  accomplished  by 
mechanically  recentering  the  feel  system. 


Figure  8.  Directional  flight  control  sygtem  integration 


Aircraft  Control  Response 
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Figure  9.  Response  to  Pitch  command,  STOL  landing  configuration 
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Speed  Hold 

The  aircraft's  response  to  a speed  command  Increase  of  10  knots  vithout  other  pilot 
action  is  presented  in  Figure  10.  As  shown,  the  speed  increase  is  smooth,  attaining 
of  the  commanded  value  in  12  seconds.  The  USB  retracts  automatically  from  32  degrees  to 
Uo  degrees.  Thrust  increases  from  32,000  to  i*8,000  pounds  during  initial  acceleration 
and  then  reduces  to  a steady  state  value  of  U0,000  pounds.  Since  pitch  attitude  is  held 
constant,  steady^state  flight  path  angle  must  increase  as  speed  is  increased.  Consequently, 
the  pilot's  task  to  acquire  a new  airspeed  consists  of  selecting  the  desired  airspeed  on 
the  CDP  and  controlling  flight  path  angle. 
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Figure  10. 
Engine  Out 


Normal  response  to  speed  command,  STOL  landing  configuration 
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Figure  H.  Response  to  engine  failure  with  pilot  "hands  off”  of  controls 
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Systems  used  normally  during  every  landing  are  satisfactory  for  controlling  the  YC-lit 
to  a STOL  landing  with  an  engine  failed.  However,  the  aircraft  Is  in  a fairly  high  drag 
configuration  with  spoilers  used  for  lateral  trim  and,  as  a result,  the  glldeslope  could 
be  flattened  only  from  -6  to  -U  degrees.  Although  -U  degrees  performance  is  adequate  to 
complete  a STOL  landing  under  good  conditions,  flight  simulator  tests  showed  it  was  not 
forgiving  if  maneuvering  was  required  to  correct  for  an  off-nominal  approach.  To  improve 
engine-out  performance,  pressure  sensors  in  each  engine  transmit  signals  to  the  FCE  for 
failure  detection.  The  FCE  then  commands  flap  retrim  to  reduce  lateral  control  spoiler 
deflection.  With  the  resultant  reduced  drag,  the  YC-l4  can  be  pulled  up  to  level  flight 
at  STOL  landing  speed  with  a failed  engine.  A landing  can  be  readily  completed  even  when 
maneuvers  are  required,  or  with  go-around  configuration  selected,  a positive  climb  gradient 
can  be  attained  without  air-speed  increase.  The  improved  performance  with  automatic  flap 
retrim  thus  increases  the  options  available  to  the  pilot. 

Conclusions 

The  YC-lU  establishes  a new  level  for  tactical  transport  flight  control.  Its  flying 
qualities  are  excellent  for  either  STOL  or  CTOL  operation,  and  its  FCS  is  a step  advance- 
ment with  its  inclusion  of  triplex,  digital  flight  control  electronics.  Digital  electron- 
ics, which  are  making  rapid  technical  and  cost  advances,  will  likely  become  standard  in 
future  flight  control  systems.  The  ability  of  digital  computers  to  perform  logic  functions 
enables  comprehensive  system  tests,  failure  monitoring  and  identification.  These  advan- 
tages will  be  reflected  in  favorable  maintenance  coats.  Finally,  the  flexibility  of 
digital  computers,  if  used  with  discipline,  enables  development  of  superb  flying  qualities 
with  fine  tuning  as  required  during  flight  teat,  and  permits  aircraft  FCS  growth  as  needed 
for  new  missions. 
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RtSUMt 

Le  contrdle  actif  du  flottement  a ete  experimente  en  soufflerie  sur  une  maquette  d'alle  munie 
d’un  reservoir  exterieur.  Les  forces  aerodynamiques  du  controle  etaient  engendrees  par  un 
aileron  classique.  pilote  par  une  servo-commande  miniature,  & partir  d'un  signal  provenant  du 
mouvement  de  I'aile.  Une  loi  de  contrdle  unique  a dte  utilisde  dans  tout  le  domaine  de  vitesse. 
Un  gain  de  plus  de  1 5 a ete  obtenu  sur  la  vitesse  critique  de  flottement. 


Descripteurs  (thesaurus  CEDOCAR)  : Vibration  structure  — Flottement  aeroelostique  — Oscil- 
lation aile  — Amortissement  vibration  — Stabilisation  — Seivocommandes. 


WIND  TII>(NEL  STUDY  OF  AN  ACTIVE  FLUTTER  SUPPRESSION  SYSTEM 


SUMMARY 

Active  flutter  control  has  been  experimented  in  a wind  tunnel  on  a model  of  wing  carrying 
an  external  tank.  The  aerodynamic  forces  of  the  control  system  were  generated  by  a classical 
aileron,  piloted  by  a miniaturized  servo-control  from  a signal  issued  by  an  accelerometer  detect- 
ing the  wing  movement.  A single  control  law  was  used  in  the  whole  velocity  range.  A gain 
of  more  than  1 5 has  been  obtained  on  the  flutter  critical  velocity. 


Descriptors  (NASA  thesaurus):  Structural  vibration  — Flutter  — Wing  oscillations  — Damping 
— Suppressors  — Stabilizers  (fluid  dynamics)  — Automatic  control  — Dynamic  control  — Servo- 
control. 


(*}  Ch«f  de  Division  de  Recherche  ii  I'ONERA 


i. 


I.  — INTRODUOTON 

II  y a d6j&  longtemps  que  Ton  a imaging  de  stabiliser  un 
avion  qui  aurait  naturellement  tendance  au  flottement  en 
utilisant  une  gouveme  a^rodynamique  asservie  au  mouve- 
ment  de  I'aile.  Mais  la  mise  au  point  de  ce  syst^me  de 
controle  actif  ezige  une  technologie  trds  d^velopp^e  dans 
les  domaines  de  I’^lectronique  et  des  servo-commandes 
^lectrohydrauliques  et  une  connaissance  precise  des  forces 
a^rodynamiques  qui  sont  & I'origine  des  flottements  ^ven- 
tuels  et  des  forces  a^rodynamiques  stabilisatrices  engen- 
dr6es  par  la  gouveme  de  controle. 

Dans  r^tat  actuel  de  la  technologie  et  des  connaissances 
en  a^rodynantique  instationnaire,  le  controle  actif  ne  peut 
pas  &tie  consid6r6  comme  une  solution  applicable  ct  tous  les 
cas  de  flottement.  II  constitue,  par  centre,  un  moyen  valable 
pour  ^liminer  des  flottements  apparaissant,  sur  une  struc- 
ture initiidement  saine,  4 la  suite  de  changements  de  confi- 
gurations dMavorables;  les  conditions  cons^cutives  k ces 
changements  de  configuration  sont  souvent  marginales 
pour  le  flottement,  et  celui-ci  peut  etre  assez  facilement 
maitris^. 

La  presence  de  charges  sous  la  voilure  des  avions  mili- 
taires  provoque  des  flottements  qui  entrent  dans  cette 
catdgorie.  En  fixant  une  charge  militaire  ou  un  reservoir 
ezt^rieur  de  grande  inertie  avec  un  support  dont  la  rigi- 
dity est  forcyment  limitye,  on  introduit  de  nouveaux  degrys 
de  liberty  qui  modifient  les  modes  propres  de  vibration 
de  la  voilure.  Certains  modes  se  dydoublent  sans  que  la 
dyformye  de  I'aile  change  fondamentalement  de  nature. 
Par  example,  le  mode  de  torsion  se  dydouble  en  gynyral 
en  une  torsion  k fryquence  relativement  yievye,  dans 
laquelle  la  charge  vibre  en  opposition  avec  I'aile,  et  en  un 
mode  y basse  fryquence  oit  la  dyformye  de  I'aile  prysente 
toujours  le  caractyre  d'une  torsion,  metis  oil  le  mouvement 
prydominant  est  I'oscillation  de  la  charge  en  phase  avec 
I'aile.  Ce  mode  se  couple  avec  la  flexion,  comme  une  tor- 
sion classique,  mais  le  fait  que  sa  fryquence  propre  soit 
plus  faible  que  la  fryquence  de  torsion  de  I'aile  lisse  contri- 
bue  k yiargir  le  domaine  de  flottement.  Toutefois,  I'oscilla- 
tion  de  la  charge,  qui  est  prydominante  dans  ce  mode,  lui 
donne  une  grtmde  inertie  (ou  masse  gynyralisye)  ; on  dit 
que  le  mode  est  lourd.  Cette  particularity  contribue  y 
diminuer  « I'explosivity  » du  flottement. 

Bien  sfir  ce  flottement,  dont  le  mycanisme  reste  celui  du 
flottement  de  flexion-torsion  classique,  n'est  qu'un  ezemple 
montrant  comment  la  prysence  de  charges  peut  aggraver 
les  couplages  ayroyiastiques. 

Mais  la  difficulty  du  probiyme  posy  par  la  prysence  des 
charges  rysulte  surtout  de  la  multiplicity  des  configurations 
de  charge  qu'un  avion  peut  emporter.  Le  flottement  dypend 
de  la  mcisse  et  des  inerties  de  ces  charges,  de  leur  rypar- 
tition  sous  I'aile  et  de  la  rigidity  des  supports,  n est  forte- 
ment  influency  aussi  par  I'interaction  ayrodynamique 
voilure-charge,  e'est-y-dire  qu'il  dypend  de  la  forme  exty- 
rieure  des  charges.  Or,  les  configurations  de  charges 
emportyes  par  un  avion  donny  sont  successivement  variyes  : 
la  rypartition  de  masse  est  souvent  asymytrique  et  elle 
yvolue  au  cours  d'une  myme  mission  par  suite  de  la  consom- 
mation  du  carburant  et  des  largages  successifs. 

Comme  il  est  pratiquement  impossible  de  ryaliser  une 
structure  qui  assure  I'absence  de  flottement  quelle  que  soit 
la  rypartition  des  charges  avec  un  devis  de  masse  raison- 
nable,  il  faut  admettre  que  I'avion  peut  se  trouver  y tut 
moment  donny  dans  une  configuration  de  charge  amenant 
le  flottement,  et  il  faut  disposer  de  mythodes  et  de  moyens 
souples  permettant  de  le  faire  disparaitre  sans  pynaliser 
la  structure  du  point  de  vue  massique. 

Le  contrdle  actif  apparait  comme  un  moyen  relativement 
simple  et  particuliyrement  bien  adapty  k cette  application. 
C'est  la  raison  pour  laquelle  I'ONERA  en  a entrepris  I'ytude 


en  se  donnant  comme  support  expyrimental  une  maquette 
d'aile  dans  laquelle  le  flottement  est  provoquy  par  la 
prysence  d'une  charge  reprysentant  un  ryservoir  extyrieur. 

Les  forces  ayrodynamiques  de  contrdle  sont  engendryes 
par  une  gouveme  de  bord  de  fuite  classique  actionnye  par 
une  servo-commande  miniature  pilotye  par  un  signal 
yiectrique  yiabory  suivant  une  certaine  loi  de  contrdle, 
y partir  du  signal  d'un  accdiyromytre  convenablement 
disposy  dans  I'aile.  La  difficulty  principale  du  contrdle 
actif  rdside  dans  la  dytermination  d'une  loi  de  contrdle 
efficace.  Celle-ci  a yty  ajustde  en  soufflerie.  Une  loi  de 
contrdle  unique  a permis  d'assurer  la  stability  dans  tout 
le  domaine  de  vitesse  explory. 


II.  — DESCRIP'nON  DE  LA  MAQUETTE 

L'aile  rectangulaire,  d'allongement  S,3,  est  fixye  y la 
paroi  (fig.  1).  L'ypaisseur  du  profil  est  de  12  %.  Un  ryser- 
voir de  grande  dimension  est  fixy  sous  la  voilure,  y 45  % 
de  I'envergure,  par  une  liaison  yiastique. 


1 

I 


I 


Fig.  1 . — Contrdle  de  flottement  sur  demi-aile  rectang ufaire  {/.  5,3). 

Soufflerie  Sj  Chalais. 


La  gouveme,  yquilibrye  autour  de  son  axe  d'oscillation, 
peut  dtre  actionnye  par  une  servo-commande  yiectro- 
hydraulique  miniature  donnant  un  couple  pur  et  pilotye 
par  le  capteur  d'accyiyration  En  I'absence  de  signal 
de  contrdle,  la  gouveme  est  verrouiliye  par  la  servo- 
commande  qui  lui  assure  une  frdquence  propre  supdrieure 
y SO  Hz. 

Aprds  ajustement  des  caractyristiques  stmcturales,  le 
calcul  pryvoyait  un  flottement  aux  environs  de  75  m/s 
en  I'absence  de  contrdle.  Les  trois  modes  retenus  dans 
rytude  du  flottement  dtaient  la  flexion  fondamentale  y une 
fryquence  de  9,30  Hz,  le  tangage  du  bidon  y 13,55  Hz  et  la 
torsion  de  la  voilure  y 18,67  Hz.  Le  calcul  montre  que  les 
autres  modes  de  la  structure  sont  suffisamment  yioignys 
en  fryquence  pour  pouvoir  6tre  ndgligys. 

Les  caractyristiques  de  ces  trois  modes  ont  dtd  dytermi- 
ndes  pfu:  un  essai  de  vibration  effectud  sans  loi  de  contrdle 
(gouveme  verrouillde  sur  la  servo-commande).  Les  ddfor- 
mdes  et  les  masses  gyndratisdes  sont  prycisdes  sur  la 
figure  2.  On  remarque  que  les  dyformdes  des  modes  2 et  3 
sont  trds  semblables  sur  I'aile  : seule  la  participation  du 
bidon  diSdre  en  phase  et  en  amplitude. 

Toutes  les  ddfoimdes  ayant  dtd  normdes  au  mdme  point 
de  I'aile,  on  constate  que  la  -nasse  gyndralisde  du  mode 
de  tangage  est  nettement  plus  grande  que  celle  des  autres 
modes.  L'aile  dtait  dqulpde  d'un  grand  nombre  de  capteurs 
accdldromytriques,  tous  rdglds  pour  avoir  la  mdme  sensi- 
bility. 


Fig.  2.  — Modes  propres  de  I’aile.  Ddformdes  normalts^es  di  1 m^Cre  au  point  A.  frequence  f,  masse  gdndralis^e  p. 


m.  — ESSAIS  EN  SOUFFLERIE 
SANS  LOI  DE  CONTROLE 

Dans  une  premiere  phase  la  vitesse  critique  de  flottement 
a d^tennin^e  exp4rimentalement  sans  loi  de  contrdle, 
et  compar^e  aux  r^sultats  des  calculs  eSectuds  par  la 
m^thode  de  doublets  et  par  une  m^thode  de  tranches 
(coefficients  bidimensionnels  coziigka  d’un  efiet  d'allonge- 
ment).  Les  deux  m6thodes  de  ctdcul  foumissent  des  r^sul- 


tats  sensiblement  identiques  et  en  assez  bon  accord  avec 
les  essais.  Les  conditions  subcritiques  ont  M6  ditenainies 
en  soutflerie  k partir  de  I'analyse  des  densit^s  spectrales 
de  la  r^ponse  de  la  maquette  k la  turbulence  naturelle  de 
la  soufflerie. 

Dans  ces  essais  la  gouveme  est  verrouill^e  sur  I'aile  par 
la  servo-coimnande.  La  comparaison  des  frequences  et  des 
amortissements  theoriques  et  experimentaux  est  donnee 
dans  la  figure  3.  On  remarque  que  la  cassure  de  la  courbe 


Fig.  3.  — Compsrsison  chiorie-experitnc*  sans  loi  de  contrMe. 


Calcul  sans  loi  Essais  sans  loi 


A 

flexion 

*0-0 

A 

flexion 

*0*3  **/oO 

X 

tanfafe  bidon 

«0  0 

X 

canfage  bidon 

*0-3  ”/oo 

o 

torsion 

*0—0 

• 

tension 

*0-5  “/oo 

7-4 


T 


Fjg.  4.  — Depart  naturel  en  flottement.  76,5  ? V ^ 77,1  m/s. 


1.  Aileron,  2,  Ampliftcateur,  3.  Incegrateur,  4,  Filtre  passe  haut,  5.  Rejecteur  basse  frequence.  6.  Rejecteur  servo. 
7.  Phase  globale.  8.  Servo-commande,  9.  Servo-valve.  10.  Amplificaceur.  11.  Sommaceur.  12.  Filtre. 


d'amoTtissement  du  mode  n°  2,  au  voisinage  de  la  vitesse 
critique,  se  retrouve  avec  des  pentes  comparables.  L’6vo- 
lution  de  I'amortissement  est  lente  jusqu'i  I'entr^e  en 
flottement,  mais  I'ezplosivit^  du  ph^nom^ne  croit  au  del^ 
de  la  vitesse  critique. 

Sur  Ja  figure  4 on  voit  un  depart  en  flottement  enregistr^ 
par  difi^rents  capteurs  de  la  maquette.  Durant  cet  enregis- 
trement  la  vitesse  4volue  tr6s  lentement  de  76,S  m /s  jusqu'i 
77,1  m/s  en  15  secondes  environ.  Lorsque  I’amplitude  priss 
par  un  point  de  r^terence  de  la  maquette  atteint  une  valeur 
donn^e,  un  dispositif  automatique  de  s6curit6  arrete  le 
flottement  en  verrouillant  le  degr4  de  liberty  de  flexion 
de  la  maquette  li  I’aide  de  cables.  Dans  I'exemple  pr^sent^, 
le  point  au  1 /4  avant  de  la  corde  d'extr^mit^  avait  au  maxi- 
mum une  amplitude  de  10  mm  4 une  frequence  de  12  Hz: 
dans  le  m6me  temps,  la  pointe  avant  du  bidon  atteignait 
80  mm. 


IV.  — ESSAIS  AVEC  LOI  DE  CONTROLE 

Le  capteur  i partir  duquel  la  loi  de  controle  est  61abor6e 
est  plac6  dans  I'aile  pr^s  de  la  ligne  des  nceuds  des  modes 
de  tangage  et  de  torsion,  et  ddtecte  done  la  r6ponse  du 
mode  de  flexion.  La  reinjection  de  ce  terme  se  fait  4 travers 


une  loi  qui  tient  compte  de  la  fonction  de  transfer!  de  la 
servo-commande,  selon  le  schema  de  la  figure  S. 

Dans  les  essais  en  sotrfflerie,  le  reglage  des  parametres 
de  phase  et  de  gain  global  de  la  chaine  a ete  recherche 
manuellement  4 une  vitesse  de  75  m Is.  II  est  interessant 
de  remarquer  que  ce  reglage  a pu  etre  maintenu  sans 
changement  dans  toute  la  gamme  des  vitesses  explorees  : 
il  assure  la  stabilite,  du  moins  jusqu'4  88  mis  qui  est  la 
vitesse  maximum  realisee,  ce  qui  represente  un  gain  de 
plus  de  15  sur  la  vitesse  critique. 

Les  essais  ont  donne  les  r4sultats  de  la  figure  6.  On  voit 
que  les  frequences  de  tangage  et  de  torsion  sont  peu 
aSectees  par  le  controle,  tandis  que  la  frequence  de  flexion 
est  considerablemem  reduite. 

Les  reactions  du  contrdle  4 la  turbulence  en  soufflerie 
provoquent  des  debattements  de  gouverne  2"  alors 
qu'une  amplitude  de  10"  etait  possible  sans  saturation. 

Pour  montrer  I'efficacite  du  systeme  de  contrOle,  on  a 
procede  au-del4  de  la  vitesse  critique  4 des  coupures 
suivies  de  remises  en  marche  du  contrdle.  La  figure  7 
montre  des  enregistrements  effectues  4 une  vitesse  de 
78,5  m/s.  Aprds  coupure  du  contrdle  et  ddpart  en  flotte- 
ment, la  loi  est  rdtablie  alors  que  l'acc414ration  4 la  pointe 
avant  du  bidon  atteint  10  g;  le  braquage  maximum  de  la 
gouverne  de  contrdle  passe  par  une  valeur  de  + 8".  Les 
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Fig.  7.  — Coupure  et  remise  en  marche  de  la  loi  de  concrSle. 


departs  en  flonement  permettent  de  determiner  la  valeur 
de  I’amortissement  negatif  pour  des  vitesses  au-del4  du 
flottement.  Sur  le  m6me  enregistrement  on  voit  un  deuxieme 
depart  en  flottement  e la  coupure  du  contrdle  et  un  arret 
par  le  systeme  de  blocage  automatique. 

A la  meme  vitesse  de  78, S m /s,  une  remise  en  marche  du 
contr&le  a ete  efiectuee  pendant  le  flottement  alors  que 
I'amplitude  de  la  gouveme  avait  deptuse  la  limite  de  satu- 
ration. Meme  dans  ces  conditions,  le  contrdle  a amorce  le 
retour  e la  stabilite,  mais  il  I'a  fait  avec  un  temps  de  reponse 
plus  long. 

II  faut  noter  que  le  choiz  de  la  phase  de  contrdle  est  trds 
critique.  Dans  le  rdglage  optimum  trouvd  e 7S  m/s,  le 
dephasage  de  1' ensemble  de  la  chaine  est  de  168"  e 12  Hz. 
La  tolerance  par  rapport  4 ce  reglage  est  faible,  puisqu'un 


ecart  de  _ 7"  amdne  des  flottements  difierents  du  flotte- 
ment initial.  Cette  difficulte  provient  beaucoup  du  fait  que 
le  dephasage  de  la  chaine  de  contrdle  n'est  pas  constant 
dans  la  bande  de  frequence  utile,  mais  varie  de  140"  4 
198"  entre  7,8  et  17  Hz.  Avec  de  telles  variations  de  phase, 
le  contrdle  peut  avoir  un  effet  stabilisateur  sur  le  mode 
dormant  le  flottement  initial,  au  voisinage  de  12  Hz,  et  intro- 
duire  de  I'energie  dans  un  mode  de  frequence  supdrieure 
ou  inferieure.  Cette  difficulte  s'ezplique  par  la  complezite 
des  fonctions  de  transferts  des  servo-commandes  eiectro- 
hydrauliques;  elle  n'est  pas  particuUere  aux  maquettes 
en  soufflerie  et  doit  se  retrouver  4 I'echelle  avion.  Elle 
empdche  de  mettre  au  point  un  systeme  fondd  sur  une  loi 
de  contrdle  simple.  Dans  le  cas  present,  on  remarquera 
que  le  dephasage  de  la  chaine  de  contrdle  dvolue  avec  la 
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{riquence  de  part  et  d'autre  de  180°.  La  portance  induite 
par  I’oaciUation  de  la  gouveme  comporte  done  une  compo- 
■ante  pr6dominante  en  opposition  avee  l'acc616ration 
ditectie  par  le  capteur  de  contrdle  Yg,  e'est-i-dire  en  phase 
avec  la  diformation  de  flexion.  Cette  force  agit  done  sur 
le  mode  de  flexion  conune  une  rigidity  negative;  elle  fait 
baisser  sa  frequence  et  I'^carte  ainsi  de  la  frequence  du 
mode  n°  2,  conune  on  le  voit  sur  la  figure  6 et  favorise  ainsi 
la  stability. 


V.  — TENTATIVE  D’IDENTinCATION 
DES  FORCES  AfRODYNAMIQUES 

L'dquation  du  flottement  s*6crit  : 

[P^[ir]  + P[P1  +[A+  Yl]  = [F^]G 
oi  P = /u. 

(X,  p,  A et  Y 8ont  respectivement  les  matrices  de  masse 
g6n6ralis6e,  de  partie  imaginaire  et  r6elle  des  forces 
a6rodynamiques  instationnaires,  et  y la  matrice  des  rigi- 


dit6s;  P^O  est  la  colonne  des  forces  induites  par  la  gouveme 
pour  un  braquage  6,  U6  par  I'inteimidiaire  de  la  loi  de 
contrdle  aux  coordonndes  gdniralisdes  q. 

Sans  loi  de  contrdle  le  second  membre  est  nul.  La  loi 
de  contrdle  est  de  la  forme  : 

+>c;g]jCe'» 

1 — tP 

oil  9 = un  ddphasage  pur  donnd  par  

1 -I-  tP 

K = gtun  pur. 

+ >C^  dtant  les  colonnes  des  forces  adrodynamiques 
induites  par  la  rotation  de  la  gouveme  6.  Le  rdglage 
est  manuel. 

A partir  de  ce  schdma,  plusieurs  difficultds  apparaissent. 
Tout  d’abord  la  fonction  de  transfert  de  I’ensemble  servo- 
valve-servo-commande  et  dlectronique  associde  n'est 
constante  ni  en  module,  ni  en  phase  dans  la  bande  des 
frdquences  considdrdes.  II  en  rdsulte  que  des  risques 
d'instabilitd  peuvent  apparaitre  pour  d'autres  modes  que 
celui  qui  est  contrdld,  ce  qui  rend  le  critdre  de  la  phase  de 
contrdle  trds  critique. 


Fi(.  8.  — Influcnc*  ds  I'idtntiflcation  >ur  la  plan  fist. 

_ , I e ealcul  «ani  corraction  _ . 1 + c»le«l  »»"•  eorraction  ru  **"•  corraction 

Toralon  j ^ ealcul  avac  eorraction  ' »"f»fa  j ^ ealcul  avae  corraction  riaxion  j ^ ealcul  avac  corraction. 
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Pour  pallier  rinauffisuice  des  theories  qui  aont  i la  baae 
des  calculs  dea  (orcea  airodynamiquea  inatationnairea, 
inauffisancea  qui  afiectent  aurtout  lea  eSeta  de  gouveme, 
U a paru  niceaaaire  de  proc^der  & une  identification  dea 
forcea  airodynamiquea  en  aouffierie.  Cette  identification 
a pour  but  de  determiner  lea  matricea  A,  B,  C et  C". 

Le  probieme  a ete  acinde  en  deux  partiea. 

a)  Determination  dea  matricea  /I  et  fi  relativea  au  plan 
fixe.  — Une  methode  de  tranchea,  dana  laquelle  on  neglige 
lea  forcea  aerodynamiquea  agiaaant  aur  le  reaervoir  et  lea 
efieta  de  couplage  aerodynamique  entre  reaervoir  et  aile, 
a ete  utiliaee  pour  calculer  lea  forcea  aerodyiumiquea 
generatricea  du  ayateme  e trois  mode#  (matricea  A et  B). 

La  reaolution  dea  equationa  montre  que  lea  coefficienta 
de  couplage  entre  modes  (coefficienta  non  diagonaux 
dea  matricea  A et  B)  ont  un  effet  negligeable  juaqu'i  tine 
viteaae  de  40  m/a,  c'eat-A-dire  que  pour  dea  viteaaea  infe- 
rieurea  e 40  m/a  on  p>eut  conaiderer  que  la  frequence  et 
I'amortiaaement  de  chacun  dea  modea  eat  determine  par 
lea  coefficients  diagonaux  correapondanta  dea  matricea  A 
et  B.  Dana  ce  caa,  lea  coefficienta  aerodynamiquea  de  tran- 
che aont  lies  par  des  relations  lineaires  aux  variations  de 
frequence  et  d’amortissement  et  peuvent  done  etre  deter- 
mines experimentalement  ail ‘on  saitmeauret  lea  frequencea 
propres  et  lea  amortiasementa  dana  le  vent  avec  une  pred- 
aion  suffiaante. 

Lea  valeurs  experimentales  ainsi  detsrminees  different 
senaiblement  des  valeurs  theoriques,  comme  le  montre  le 
tableau  I dans  lequel  : 


It,  — Jr ' -I-  jWgk'J'  eat  le  coefficient  de  portance  dfi  i la 
deflexion  verticale, 

kf,  — k'l,  + jw/fk'^'  eat  le  coefficient  de  portance  dfl  e la 
deOexion  de  torsion, 

=■  + /e^nt'"  eat  le  coefficient  de  moment  dO  k 

la  deflexion  verticale, 

nij,  = -t-  eat  le  coefficient  de  momen'  dfi  4 

la  deflexion  de  torsion, 
cig  eat  la  pulsation  reduite  du  moment, 
kj  = k'j  + Ji^gk'J'  eat  le  coefficient  de  portance  induit 
par  la  rotation  de  la  gouveme, 

irij  = iti'j  + Jbtgm'j'  eat  le  coefficient  de  moment  induit 
par  la  rotation  de  la  gouveme, 

eat  le  coefficient  de  moment  de  la 

gouveme. 

TABLEAU  I 

Cotfficianu  ralatift  i I'ailt 


fc. 

• 1 

*t'»  + 

m, 

Ujl  V 1 

m* 

k-. 

k * + 

m\ 

Th4ori« 
lr«.  1| 
Exp4rttnce 

0,054 

0,050 

0,943 

0,880 

m rs 
oo 
1 1 

0,88 

0,94 

1,58 

1.33 

1.08 

0,79 

Dans  le  calcul  des  valeurs  ezp^rimentales  on  a suppose 
que  les  coefficients  (')  et  ("')  ^talent  ind^pendants  de  la 
frequence  r^duite.  L'iniluence  des  coefficients  corrigte 
SUT  le  Oottement  est  montr^e  figure  8. 

b)  La  seconde  identification  porte  sur  les  terines  c'  et  c" 
de  la  colonne  de  contrdle.  EUe  a 6t6  r^alis^e  par  la  mesure 
des  champs  de  pressions  instationnaires  induits  par  une 
rotation  de  la  gouverne  ft  diff^rentes  frequences  et  difie- 
rentes  vitesses.  La  comparaison  avec  la  distribution  de 
instationnaire  calcuiee  par  une  methode  de  doublets  est 
donnee  sur  la  figure  9.  La  corde  choisie  est  la  corde  mediane 
de  la  gouverne.  Les  pressions  etaient  mesurees  ft  I'eztrados 
en  10  points. 

Les  differences  entre  theorie  et  experience  sont  impor- 
tantes,  sur  le  module  atissi  bien  que  sur  la  phase. 

L'integrale  des  pressions  ezperimentales  donne  des 
coefficients  de  force  et  de  moment  qui  peuvent  §tre 
compares  avec  les  valeurs  theoriques. 


TABLEAU  II 

Coefficients  relatifs  i la  gouverne 


k ,/ 

m ,/ 

n rt 

k .1 

m ,1 

r ./ 

Theorie 
(ref.  1| 
Essais 

0.51 

0.36 

0.45 

0.32 

0,022 

0.015 

- 0,020 
+ 0,058 

0,20 

0,17 

0,025 

0,14 

Le  recoupement  entre  theorie  et  experience  porte 
egalement  sur  les  fonctions  de  transfert  de  la  maquette  e 
differentes  vitesses.  Ces  fonctions  de  transfert  expriment 
le  rapport  entre  la  reponse  de  I'aile  et  le  braquage  de  la 
gouverne.  Elies  ont  ete  determinees  experimentalement 
avec  une  excitation  par  bruit  blanc  sur  la  gouverne. 

Le  calcul  a 6t^  effectu^  avec  les  coefficients  exp4rimen- 
taux.  La  figure  10  montre,  k une  vitesse  de  6S  m /s,  la  bonne 
concordance  entre  calcul  et  essai  pour  le  point  de  I'aile 
ou  est  plac^  le  capteur  qui  foumit  les  informations  i la  loi 
de  contrdle  (fonctior  de  transfert  en  boucle  ouverte). 
On  peut  observer  qu’en  ce  point  la  rdponse  du  mode  de 
flexion  / = 9,8  Hz  est  prddominante  devant  les  autres  modes. 

L'introduction  de  la  loi  de  contrdle  expdrimentale  dans  le 
calcul  corrigd  donne  finalement,  en  fonction  de  la  vitesse, 
une  Evolution  satisfaisante  des  paramdtres  de  frequence  et 
d'amortissement.  En  particulier  I'effet  de  rigiditd  adrody- 
namique  negative  qui  dcarte  le  mode  1 du  mode  2 est  bien 
reprdsentd  figure  11.  Le  calcul  prdvoit  que  la  maquette 
reste  stable  jusqu’4  la  vitesse  maximum  expdrimentde, 
ce  qui  est  conforme  ft  I’expdrience.  Afin  de  verifier  la 
sensibilitd  de  la  loi  de  contrdle  k une  variation  de  phase, 
le  calcul  a dtd  effectud  en  introduisant  un  ddphasage  suppld- 
mentaire  de  + 5"  dans  la  loi. 


L'amortissement  du  mode  de  tangage  k 88  m/s  varie 
dans  ces  conditions  de  + 8 (valeur  avec  la  loi  nominale) 
d 2 De  mdme  un  ddcalage  de  phase  de  — 8"  conduit  d 
un  amortissement  de  11  Cette  trds  grande  sensibilitd 
aux  paramdtres  de  contrdle,  qui  est  conforme  d I’expdrience 
est  aggravde  par  la  similitude  des  ddformdes  d'aile  dans 
les  modes  de  tangage  et  de  torsion  qui  rend  la  ma trice 
modale  presque  singulidre. 


VI.  — CONCLUSION 

On  a montrd  qu’il  dtait  possible  de  contrdler  un  flottement 
provoqud  par  la  presence  d’un  rdservoir  de  grande  dimen- 
sion en  utilisant  une  gouverne  de  bord  de  fuite  classique 
pour  engendrer  des  forces  adrodynamiques  stabilisatrices. 
La  stabilitd  a dtd  obtenue  avec  une  loi  de  contrdle  unique 
dans  toute  la  gamme  de  vitesse  explorde  (de  0 d 88  m /s) 
alors  que  le  flottement  se  produisait  d partir  de  73  m /s 
en  I'absence  de  contrdle. 

Du  fait  de  la  complexitd  des  fonctions  de  transfert  de  la 
servo-valve  et  de  la  servo-commande,  il  est  difficile  de 
rdaliser  un  systdme  de  contrdle  dont  le  fonctionnement  se 
prdte  d une  interprdtation  simple,  n est  dgalement  difficile 
de  prdddtermiper  la  loi  de  contrdle  d cause  des  incertitudes 
qui  existent  dans  le  calcul  des  forces  adrodynamiques 
instationnaires  et,  surtout,  des  forces  adrodynamiques 
induites  par  le  braquage  des  gouvernes.  Dans  le  cas  prd- 
sent,  la  meilleure  loi  de  contrdle  a dtd  obtenue  par  un 
rdglage  manuel  effectud  en  soufflerie;  mais  cette  solution 
n'a  dtd  possible  que  parce  qu'il  n'existait  qu'un  rdglage 
de  phase  et  un  rdglage  de  niveau. 

Par  centre,  il  est  assez  encourageant  de  constater  que 
I'identification  des  forces  adrodynamiques  effectude  en 
soufflerie  a permis  de  corriger  les  forces  adrodynamiques 
thdoriques  et  de  retrouver,  par  un  calcul  a posteriori,  des 
dvolutions  de  frdquence  et  d'amortissement  et  des  fonctions 
de  transfert  en  assez  bon  accord  avec  I'expdrience. 

Manuscrit  remis  le  17  Juin  1976. 


RfirtRENCES 

PI  CHOPIN  S.  et  SALAUN  P.  — Coefficients  oirodynamiques 
instationnaires  thEoriques  en  regime  subsonique  pour  une  voilure 
de  faible  allongement.  — Document  ONERA  non  publid  (1959). 

|2|  HAIDL  G.  — Active  flutter  suppression  on  wings  with  external 
stores,  denstt  Active  control  systems  for  loads  alleviation,  flutter 
suppression  and  ride  control  ».  AGARDograph  n"  175  (1974), 
p.  57-76. 

|3|  SENSBURG  O.  and  HONLINGER  H.  — Dynamic  testing  in 
wind  tunnels,  dans  « Flight  ground  testing  facilities  correlation. 
AGARD  Conf.  Proc.  n"  187  (1976),  mdmoire  n"  5,  part.  I. 

|4|  NIS5IM  E.  — Active  flutter  suppression  using  trailing  edge  and 
tab  control  surfaces.  — AIAA  paper  n"  75-822, 


REPORT  DOCUMENTATION  PAGE 


1.  Recipient’s  Reference 

2.  Originator’s  Reference 

3.  Further  Reference 

4.  Security  Gassification 
of  Document 

AGARDCP-228 

ISBN  92-835-0200-0 

UNCLASSIFIED 

S.Origiiuitor  Advisory  Group  for  Aerospace  Research  and  Development 


North  Atlantic  Treaty  Organization 
7 rue  Ancelle,  92200  Neuilly  sur  Seine,  France 

6.Title 


STRUCTURAL  ASPECTS  OF  ACTIVE  CONTROLS 


7.  Presented  at  44th  Meeting  of  the  AGARD  Structures  and  Materials  Panel  held  in  Lisbon, 
Portugal  on  21  April  1977. 


8.Author(s) 


9.  Date 


Various 


August  1977 


10.  Autfior’s  Address 


11. Pages 


Various 


12.  Distribution  Statement  This  document  is  distributed  in  accordance  with  AGARD 

policies  and  regulations,  which  are  outlined  on  the 
Outside  Back  Covers  of  all  AGARD  publications. 

1 3.  Keywords/Descriptors 

Structures  Design  criteria  Feedback  control 

Servomechanisms  Transfer  functions 

Adaptive  systems  Tests 


106 


14.  Abstract 

^^.^he  Meeting  dealt  with  the  philosophy  and  approach  on  the  use  of  active  controls  to 

realize  structural  improvements.  The  question  of  what  constitutes  a good  balance  of  effort 
to  achieve  a successful  active  control  system  was  posed.  Specifics  dealt  with  the  techniques 
for  evaluating  the  system  transfer  function,  with  the  relative  roles  of  ground  vibration 
testing,  bench  testing  of  component  parts  and  the  merits  of  open  and  closed  loop 
testing.  The  question  of  what  is  an  appropriate  index  of  performance  is  of  central 
significance.  /A 


AG ARD  Conference  Proceedings  No.  228  AGARD-CP-228  AGARD  Conference  Proceedings  No.  228  AGARD-CP-228 


S'"  C 75 

EE  2 2 

s a -g  5 = 

I § s 

S u a.  C ^ M 

Sa  44  ^ O O 

« S g-.sp  c i3  ^ 

c C JS  CQ  ^ 44 

w 44  ^ ^ U 44  44 

C/3  C/)  < O H H tu 


52  C *o 

E e .2  2 

.2  W ea  - 

C*^  'r*  O 5 

5 >»  « 5 u 

s "S  j;;  *c  ^ 

U 44  If  CJ  u C4 

S E -S  c -2  2 

a g g-.2f  c -o 

t ^3  Si  ?L  & 

;/j  Vi  ^ Q H H U. 


■ S B ;=  '5 

e ^ 
44  O s -2 

.a  « o 2 
*e3  3 

* 44  W «4J 

*-  j=  2 

'o  » 

^ > « 
® .2  Cu 
‘O  e -C 

• u 5 o 
tS  2 


c -3 

.leg 

u ^ c 
c . 54 

3 M Oi 

C o 1 


:§ 

? B « <2 

i 

^ Urn  ^ 

4>  U.  o 

•o  o . ® g 

an  «2  44  44 

•S  I 5 = ^ 

« = E-2  - 

^lirs 

« „ a o c 
£ E E o o 
H o .a  BO  o 


c ^ , 

««  g *r  ! 

- -2  « i 

2 2 f I 

« .0  « ' 
w ■>  S j 
y-g  7>  •. 

« § g 

- 2 <2 

c *2  2 J 

-00.^ 

3 44  C 
*e3  O 2 I 

S u O' 

g ai 

u .s  e 
o ^ c , 

«£  ««  o 
«i  "5  o 

44  ^ 

S 44  O 
O'JS  ^ 

S - a’ 

^ ^ W 

^ .tS  ' 

W J V - 


cO  -' 

cS  I 

Q b o. 

OS  g.2 


! <Q  44  44 
i -C 

«/5 

'•3 

: *1  3 ’i  * 


j ♦-  3 

4 ea  c4) 

* *5  I-  o 

> ^ S 


5^0*^ 

> 44  44 

3 ^ 

- a- 

► O 

J iC  ^ rt 

^ p « ^ ■ 


a:  < 

5 ^ «/5 
^ *t3  44 
P 44  ttO 
O J=  SO 

D .2  3- 

^ .0  sO 

H 3 O 

C/3  Ou  ^ 


E « £ b w-  -g 

3 E -B  c <2  « 

u o o.  so  2 VI  ^ 

g E:  j9  '5!  M t;  u 

3 S "o  u g u 5; 

VI  </3  <;  Q H H U. 


2 VI  c 

E E .0 

.i2  44  cQ  M 

g to  ‘C  2 
2 >*  ii  § 

U5  »3  .ts  t2 

2i  44  £ b 

2 E B c >52 

C4  O D.  00  ^ (O 

E E:  .g  8 « »; 

*3  44  ^ w u,  4# 

C/3  C/3  ^ Q H H I 


>p?s « 


’ 25  e 

- .c  E (j 
; O c 

> _ '-S  44 

) X C4  .0 

.t:  c 
^ .3  00 

i ^ ^ c 

z: 

, — .44  VJ 
u 2 »5 
! "O  S c 

‘ VI  2 o 
1.2  " 
S E 2 

> O 44 

• 44  »•  .m 

I Q.  > 
: c/5 

> ..  c 


ja-s's 

> CO  to 

1^3^ 

' 2 « p 

5 > V. 

’ E 

> 44 ,0  t; 

» W CO 

i S>s  K 

I 3 44 
) -or  3*4C 
O -3  ••-» 
l te  3 — 


*2  jL  V 
S o Si 
i:  c 

e CL 

I 0.  " is 
o c c 

CO  00 

0 .j2 

M -M 

•^S  ^ 2 
n £ ^ 

S ? C 
E tZ  44 

44  O 
.C  r-  ^ 

* § o 

•O  r t/i 
c ^ • 

3 44 

ir  °'c 

Sum 

-fSi 

1 Ob'S 

o .s  a 

Q,  *ii 

E S-s 
O 

a X 

2 iJ 
o ^ ^ 

00  w ‘E 
3 Ps 

•—  £ 44 
nl  .2  W 

2 u c 


f » »• 

<3  1 

o er  o. 
2 ^ 

< '?  £ 
O *0  V 

<<a 


^ E 

oc  < 

^ t/5 

b«  ^ 44 

n 44  00 

y 2 

3 <2  CL 


so  44  44 

> >3 
v>  *.S 

44  o 
*2  CO  .C 

i a * 

g ^ ii 

S 8 3 

S 

^ ^ t/j 

■g  o 

» « i- 

44  *3 
O 44  CL 
^ S 5/5 


t/5  O 
44  ^ V 

3 2 

«yt-  2 


• ® E 

44  44 

; F c ^ 

' E — ^ 

ii:i 

CL  O 3 
: e O O 
.S  00  C4 


techniques  for  evaluating  the  system  transfer  function,  techniques  for  evaluating  the  system  transfer  function, 

with  the  relative  roles  of  ground  vibration  testing,  bench  with  the  relative  roles  of  ground  vibration  testing,  bench 

testing  of  component  parts  and  the  merits  of  open  and  testing  of  component  parts  and  the  merits  of  open  and 

closed  loop  testing.  The  question  of  what  is  an  approp-  closed  loop  testing.  TTie  question  of  what  is  an  approp- 
riate index  of  performance  is  of  central  significance.  riate  index  of  performance  is  of  central  significance. 


NATO  ^ OTAN 

7 RUE  ANCELLE  • 92200  NEUILLY-SUR-SEINE 
FRANCE 


DISTRIBUTION  OF  UNCLASSIFIED 
AGARD  PUBLICATIONS 


Telephone  745.08.10  ■ Telex  610176 


AGARD  does  NOT  hold  stocks  of  AGARD  publications  at  the  above  address  for  general  distribution.  Initial  distribution  of  AGARD 
publications  is  made  to  AGARD  Member  Nations  through  the  following  National  Distribution  Centres.  Further  copies  are  sometimes 
available  from  these  Centres,  but  if  not  may  be  purchased  in  Microfiche  or  Photocopy  form  from  the  Purchase  Agencies  listed  below. 


BELGIUM 

Coordonnateur  AGARD  - VSL 
Etat-Major  de  la  Force  Adrienne 
Caserne  Prince  Baudouin 
Race  Dailly,  1030  Bruxelles 

CANADA 

Defence  Scientific  Information  Service 
Department  of  National  Defence 
Ottawa,  Ontario  KIA  OZ2 

DENMARK 

Danish  Defence  Research  Board 
Osterbrogades  Kaseme 
Copenhagen  Q 

FRANCE 

O.N.E.R.A.  (Direction) 

29  Avenue  de  la  Division  Leclerc 
92  Chatillon  sous  Bagneux 

GERMANY 

Zentralstelle  fiir  Luft-  und  Raumfahrt- 
dokumentation  und  -information 
Postfach  86088L 
D-8  Miinchen  86 

GREECE 

Hellenic  Armed  Forces  Command 
D Branch,  Athens 

ICELAND 

Director  of  Aviation 
c/o  Flugrad 
Reykjavik 


NATIONAL  DISTRIBUTION  CENTRES 
ITALY 

Aeronautics  Militate 

Ufficio  del  Delegato  Nazionale  all'AGARD 
3,  Piazzale  Adenauer 
Roma/EUR 
LUXEMBOURG 
See  Belgium 
NETHERLANDS 

Netherlands  Delegation  to  AGARD 
National  Aerospace  Laboratory,  NLR 
P.O.  Box  126 
Delft 
NORWAY 

Norwegian  Defence  Research  Establishment 
Main  Library 
P.O.  Box  25 
N-2007  Kjeller 
PORTUGAL 

Direccao  do  Servico  de  Material 

da  Forca  Aerea 

Rua  de  Escola  Politecnica  42 

Lisboa 

Attn:  AGARD  National  Delegate 
TURKEY 

Department  of  Research  and  Development  (ARGE) 
Ministry  of  National  Defence,  Ankara 
UNITED  KINGDOM 

Defence  Research  Information  Centre 
Station  Square  House 
St.  Maty  Cray 
Orpington.  Kent  BR5  3RE 


UNITED  STATES 

National  Aeronautics  and  Space  Administration  (NASA), 

Langley  Field,  Virginia  23365 

Attn:  Report  Distribution  and  Storage  Unit 


THE  UNITED  STATES  NATIONAL  DISTRIBUTION  CENTRE  (NASA)  DOES  NOT  HOLD 
STOCKS  OF  AGARD  PUBLICATIONS,  AND  APPUCATIONS  FOR  COPIES  SHOULD  BE  MADE 
DIRECT  TO  THE  NATIONAL  TECHNICAL  INFORMATION  SERVICE  (NTIS)  AT  THE  ADDRESS  BELOW. 


PURCHASE  AGENCIES 

Microfiche 
Technology  Reports 
Centre  (DTI) 

Station  Square  House 
St.  Maty  Cray 
Orpington,  Kent  BR5  3RF 
England 

Requests  for  microfiche  or  photocopies  of  AGARD  documents  should  include  the  AGARD  serial  number,  title,  author  or  editor,  and 
publication  date.  Requests  to  NTIS  should  include  the  NASA  accession  report  number.  Full  bibliographical  references  and  abstracts 

of  AGARD  publications  are  given  in  the  following  journals: 

Scientific  and  Technical  Aerospace  Reports  (STAR),  Government  Reports  Announcements  (GRA), 

publiidied  by  NASA  Scientific  and  Technical  published  by  the  National  Technical 

Information  Facility  Information  Services,  Springfield 

Post  Office  Box  8757  Virginia  221 51 . USA 

Balthnoie/WaAington  International  Aitport 
Maryland  21240,  USA 


Microfiche  or  Photocopy 
National  Technical 
Information  Service  (NTIS) 
5285  Port  Royal  Road 
Springfield 
Virginia  22151,  USA 


Microfiche 

Space  Documentation  Service 
European  Space  Agency 
10,  rue  Mario  Nikis 
75015  Paris,  France 


tinted  by  Technical  Editing  and  Reproduction  Ltd 
Harford  Houte,  7-9  Charlotte  St,  London  WIP  ]HD 


ISBN  92-83S-0200-0 


